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v oltage stability has been defined [1] as the ability of the system to maintain voltage so that 
when load admittance is increased, load power also increases and both power and voltage re- 
main controllable Voltage collapse is the process by which voltage instability leads to very 
low voltage profile in a significant part of the system resulting into partial or total blackout. 

oltage stability studies have been carried out for variation (increase) in system loading or for 
c ontingency conditions 

Most of the works have considered voltage stability as static phenomena due to slow variation 
of voltage over a long time until it reaches near to maximum loadability or collapse point Load 
flow equations have been used as model and to explain the phenomena, singularity of Newton- 
Raphson load flow Jacobian has been observed [2-4]. Minimum singular value and condition 
number of the load flow Jacobian [5-7] have been popularly used as indicators to predict static 
voltage stability in A C systems Several other techniques such as optimisation method [8], 
multiple load flow solutions [9] etc have been used to predict extreme loading conditions and 
static voltage stability margin 
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From the literature survey, it appears that most of the works [10,11] on voltage stability 
analysis of AC-DC systems have represented AC systems by their Thevenm equivalent The 
concept of minimum singular value and condition number has not been applied to combined 
AC-DC systems considering the detailed representation of AC systems 

One of the recent concerns of maximising voltage stability margin has led the researchers 
to find out new objectives for rescheduling real and reactive output of the sources An optimal 
generation scheduling scheme with respect to maximisation of minimum singular value of the 
Jacobian has been presented in ref [12] Traditionally, optimal power flow (OPF) is run for 
generation rescheduling m view of minimisation of the total fuel cost of generation or minimi- 

s 

sation of the system transmission loss From the literature survey it appears that the effect of < 
the various optimal power flow schemes on the voltage stability margin has not been explored 
One important outcome of the studies reported in the literature is the general consensus 
that bifurcations in underlying mathematical models of a power system are closely linked with 
various types of instabilities This is especially true for voltage collapse precipitated by the slow 
variation in a system parameter such as load etc A bifurcation is a qualitative change in the 
phase portrait of a dynamical system that occurs as a system bifurcation parameter is varied 
Bifurcations and chaos have been demonstrated by several researchers [13] m simple models 
of power system Hopf bifurcation phenomena has been studied in ref [14-16]. A fourth order 
generator model along with excitation system have been considered in ref [15] and it has been 
shown that a complex pair of eigenvalue associated with excitation system results m Hopf 
bifurcation Most of the works [17, 18] showing chaotic invariant set have considered only the 
load dynamics along with swing equations and have observed only the period doubling route to 
chaos 

The literature survey reveals that the detailed model of induction motor with speed depen- 
dent generalised loading and speed governing control loop has not been considered Also the 
dynamics of various other dynamic elements such as OLTC and SVS have not been investigated j 
The bifurcation and chaos are unwanted phenomena in the power system Hopf bifurcation 
reduces the stability margin, especially the subcritical Hopf bifurcation in which the operating 
point is stable but its region of attraction is reduced by the surrounding unstable periodic 
orbit A power system though being dynamically stable can exhibit a bounded behaviour when 
the stable operating print is perturbed to attracting region of chaos The chaotic systems 
have continuous broad band spectrum and can introduce harmonics especially low frequency 
harmonics m the system, even in the absence of any harmonic source 
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One of the simplest ways to avoid bifurcation and chaos is to move the system out of the 
chaotic region by restricting the parameter to enter this region However, it reduces the system 
loadability Recently new devices for the control of voltage and power have been developed [19] 
The prominent amongst these components are static var compensators (SVC), static phase angle 
regulator (PAR) and controllable services capacitors (CSC) The effectiveness of these devices 
for controlling bifurcation and chaos has not been explored 

Therefore, the motivation behind the work reported in this thesis are 

( I ) To extend the concept of minimum singular value and condition number of load flow Ja- 
cobian to combined AC-DC systems considering the detailed representation of AC system 
and also to AC systems having voltage dependent loads 

( II ) To explore the effect of different generation rescheduling schemes on the static voltage 
stability margin considering voltage dependent capabilities of synchronous machine 

(m) To observe various types of local bifurcations such as static and Hopf bifurcations and their 
effect on stability considering the dynamic models of generator and loads, to determine 
critical loading and control parameters’ values and to ascertain the predominant states 
causing Hopf bifurcation using participation analysis 

(iv) To observe various types of global bifurcations and chaos considering the OLTC, SVS, 
load and generator dynamics and possibility of voltage collapse due to these phenomena 

(v) To study the effectiveness of FACTS devices in eliminating dynamic bifurcation and chaos 

A brief description of the work reported m the thesis is given below 

Chapter 1 introduces various aspects of voltage stability and presents the state-of-art survey 
on the subject 

In Chapter 2 the minimum singular value and condition number of load flow Jacobian has 
been used to predict static voltage stability of integrated AC systems considering voltage de- 
pendent characteristics of loads and extended, for the first time, to combined AC-DC systems 
considering detailed representation of AC systems. The effect of control mode switching in the 
DC link and bus type switching in handling the reactive power limits of generators on the static 
voltage stability has also been demonstrated 

In Chapter 3 four different schemes of generation rescheduling such as minimisation of total 
cost of generation, system transmission loss, maximisation of minimum singular value and a 
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new scheme have been used to maximise the static voltage stability margin Voltage dependent 
capabilities of synchronous machines have been considered m the model The new proposed 
scheme of generation rescheduling has consistently resulted into maximum voltage stability 
margin for various loading conditions 

Chapter 4 investigates the local bifurcations namely static and Hopf bifurcation in various 
sample power systems Effect of static bifurcation on power system stability has been demon- 
strated by studying the associated unstable manifold of type-1 unstable equilibrium points on 
the stability boundary Center manifold theory has been used for studying the stability of Hopf 
bifurcation and determining critical values of parameters The effect of local bifurcations on the 
stability in parameter space has also been studied 

Chapter 5 presents the sample system study results on observance of various types of global 
bifurcations such as period doubling bifurcation, bifurcation of periodic orbit to torus, chaos via 
period doubling bifurcations & via torus breakdown and boundary crisis Lyapunov exponents 
and dimension of the attractor have been determined Broad band frequency spectrum has 
been observed under chaotic oscillations Occurrence of voltage collapse due to subcritical Hopf 
bifurcation and boundary crisis has been demonstrated 

Chapter 6 demonstrates the application of FACTS devices such as PAR, SVC and CSC 
for control of bifurcation and chaos A first order delay model of these devices have been 
considered It has been found that the dynamic bifurcations and chaos can be eliminated with 
proper selection of control signals and the controller gam parameters 

Chapter 7 concludes the main findings m this thesis and gives suggestions for further work 
in this area 
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Chapter 1 


Introduction 


1.1 General 

A modern interconnected power system network can be described as highly nonlinear and 
stressed system. Bemg one of the most complex man made dynamical system, its stability 
and control has always been challenging The power system stability may be broadly defined as 
the property of the system that enables it to remain in a state of operating equilibrium under 
normal operating conditions and to regain an acceptable state of equilibrium after bemg sub- 
jected to a disturbance [176,189] Instability in a power system may be manifested m many 
different ways depending upon the system configuration and operating mode Traditionally, the 
problem of stability has been one of maintaining synchronous operation — rotor angle stabil- 
ity However, due to economic and environmental pressures the power system of today is bemg 
operated close to the maximum loadability limit, and is facing the problem of maintaining the 
required bus voltage [92, 130,146] — voltage stability 

1.1.1 Voltage Stability, Voltage Collapse and Voltage Security 

Some of the terms related to voltage stability have been formally defined m IEEE report [73] 
and by Concordia [82] which are given below 

• Voltage stability is the ability of a system to maintain voltage so that when load admittance 
is increased, load power will increase, and so that both power and voltage are controllable 

• Voltage collapse is the process by which voltage instability leads to very low voltage profile 
in a significant part of the system 


1 
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• Voltage security is the ability of a system, not only to operate stably, but also to remain 
stable (as far as the maintenance of system voltage is concerned) following any reasonable 
credible contingency or adverse system change. 

Voltage stability encompasses a wide range of phenomena. It has been treated as fast phenomena 
by engineers involved with induction motor loads or HVDC links It has also been viewed as 
slow phenomena by those involved with on load tap changer dynamics of transformers 

1.1.2 (Voltage) Angle Vs Voltage (Magnitude) Stability 

The problem of rotor angle instability is well understood and documented [7, 14,25, 109, 176] 
However, the voltage instability has recently attracted the attention of researchers world wide 
and is an active area of research 

Just as real power is the key variable to rotor angle instability analysis, reactive power is 
central to voltage stability analysis Deficit of reactive power either locally or globally leads to 
poor voltage profile and with increase in loadings it may lead to voltage collapse 

Imbalances in active power input to the generators and their electrical power outputs are 
balanced by the changes in angular momentum of the rotor resulting into change m the average 
system frequency, whereas the imbalances in the reactive power result into change m the bus 
voltage profile affecting the energy stored in shunt capacitors and reactors Usually the R/X 
ratio of the transmission system is very small Because of low resistances of the network the 
flow of active energy between the generators and the loads is lightly damped However, because 
of large reactance of the transmission system, the flow of reactive power is highly damped, 
causing the voltage instability to start as a local phenomenon The angle instability leads to 
loss of synchronism and operation of protective relays causing disconnection of loads and/or 
system breaking into separate subsystems The problem of angle instability has now been 
overcome due to improvements in the design philosophy of machines and control strategy such 
as power system stabilisers, fast excitation system etc The loss of voltage (magnitude) stability 
follows upon a change in either the reactive power supply or demand beyond the ability of 
the system to accommodate it by changes in bus voltage magnitude resulting into stalling of 
induction motors and/or disconnection of equipments due to protective relay operation leading 
to cascade tripping The power factors of the generators affect the voltage stability Under 
voltage instability conditions, some of the existing controls may prove counter productive 
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1.2 State-of-the-Art 

Voltage collapse was identified and analysed for simple power systems even in the first half 
of this century as given in text books [1,3,4] However, only recently the problem of voltage 
stability has become serious operating concern The voltage stability studies have been carried 
out on static as well as on dynamic considerations A large number of papers exist in the area 
of voltage stability and it is difficult to present a detailed review of all of them Hence, only 
a representative survey of literature relevant to the work carried out in the present thesis is 
presented in this section 

Mercede et al [48] and Chow et al [67] presented a framework for analysis and prediction 
of voltage stability problems Depending upon the nature of disturbance and the subsystem 
dynamics, the voltage instability may be regarded as slow or fast phenomena [121] A special 
issue on the subject [129] contains some of the state-of-the-art papers 

The static or steady state voltage stability has often been regarded as viability problem 
suitable for static power flow 1 The ability to transfer reactive power from production sources 
to consumption nodes during steady state operating conditions is a significant aspect of voltage 
stability A 1987 CIGRE report [42] recommends analysis methods and power system planning 
approaches based on static models The singularity of power flow Jacobian has been linked to 
steady state stability [6,36,78] and its positive definiteness has been proved to be the sufficient 
condition for an operating point to be small disturbance stable [41,56] Various indices for 
prediction of steady state voltage stability margin has been presented m CIGRE report [169] 

It has been observed that when the operating point reaches the nose point of P— V or 
Q—V curve, the sensitivity of bus voltage to the load power becomes alarmingly high and at 
the nose of P—V or Q—V curve, jp or ^ becomes infinity [8] The concept of load flow 
Jacobian sensitivity has been used by Flatabo et al [69, 134] and has been applied on the 
Norwegian system for maintaining required security level The sensitivity of total reactive power 
generation to variation in real and reactive load has also been proposed as security index by 
Carpentier [28] and Begovic et al [102] Carpentier has termed this index as Voltage Collapse 
Proximity Indicator (VCPI) Another sensitivity based indirect criteria which relates the change 
in generator voltage E with corresponding change in load bus voltage Vj(^) has been proposed 
by Venikov [8] and used by Borremans et al [27] for voltage stability studies At the point of 
voltage instability, ^ becomes zero, meaning thereby the load voltage can not be controlled by 


1 Load flow and power flow have been interchangeably used in this thesis 
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generators, a fact which has been mathematically proved by Hiskens et al. [56,84]. The deficit of 
reactive power has been cited as one of the important reasons responsible for voltage instability 
by Lachs [31]. Tamura et al. [24] observed that voltage instability is closely related to multiple 
load flow solutions which are caused by the nonlinear nature of load flow equations. A heavily 
loaded system has two very close power flow solutions. One is the higher voltage power flow 
solution which is usually obtained by the conventional power flow calculations and the other is 
the lower voltage solution which is obtained by the methods proposed in ref. [10]. When these 
two solutions coalesce, a saddle node bifurcation occurs. The proximity of these two solutions 
have been used by many researchers as an index for voltage stability margin. Use of energy based 
measures has been suggested to determine the closeness of multiple solution [68,90,91,155]. Use 
of P—V and Q-V curves have also been made to assess the vulnerability of an operating point to 
voltage collapse [49,53,75,94,97,103,116,126,181]. Using a conventional power flow programme, 
the maximum power transfer can be found by a binary search without computing an entire P—V 
curve. Near the nose point, the Newton-Raphson algorithm may diverge. To overcome this 
difficulty, the continuation power flow method has been introduced by Ajjarapu et al. [100], 
Canizares et al. [166] and Iba et al. [86]. 

Tire divergence of power flow solution has been used as one of the indicators for assessing 
voltage instability [37,71]. However, the divergence of power flow may also be due to the load flow 
Jacobian being ill-conditioned. The proximity of load flow Jacobian to singularity can be better 
measured in terms of its minimum singular value [51,52,58,113,141-143,179]. A fast algorithm 
based on parallel processing for computation of minimum singular value has been suggested 
by Tiranuchit et al. [52]. A minimum singular value based security index for scheduling real 
and reactive power output of sources has also been suggested by Tiranuchit et al. [51]. Lof et 
al. [113,141] presented a fast algorithm based on inverse iteration to compute minimum singular 
value and tested it on a large system. Prediction of static voltage stability has also been made 
by Lof et al. [142, 179] using minimum singular value of power flow Jacobian and considering 
voltage dependent reactive capability of generators. Pai et al. [58] considered an exponential 
type voltage dependent load model and considered a new index called condition number for 
static voltage stability prediction. Begovic et al. [81] demonstrated that by minimising the 
weighted sum of the absolute values of control actions subject to equality constraints (real and 
reactive power balance equations) and inequality constraints (limits on the state variables, real 
and reactive power generations) voltage stability could be improved. 

Eigenvalue analysis has been used to assess static voltage stability and sensitivity analysis has 
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been earned out to determine critical group of buses responsible for voltage collapse [64,93, 110, 
148,175] Participation analysis has also been carried out to determine participation of various 
states m different modes of power system model by Hamdan [135] and Rajagopalan et al [61] 
The concept of coherency among load buses with respect to voltage dynamics was introduced 
by Begovic et al [65] and Vargas et al [159] An algorithm for approximate assessment of 
minimum singular value [65] and eigenvalue [159] of Jacobian was presented based on partial 
state information about network obtained from coherent clusters of load buses Ajjarapu et 
al [162] have presented a scheme to improve voltage stability margin by use of fixed capacitors 
at appropriate locations However, the use of fixed capacitors may prove detrimental to voltage 
stability in some operating conditions (example 11 2 of ref [8] and example 14 1 of ref [176]) 

The prevailing approach to real time control of voltage collapse (as opposed to system plan- 
ning) consists of security monitoring coupled with algorithms for allocation of resources for 
minimising risk of voltage collapse [102, 140, 152] A related problem is that of maintaining max- 
imum distance (in parameter space) to the set of parameter values at which bifurcation from 
the nominal operating point occurs The development of voltage stability indices is essential 
to studies of this kind [67, 107, 132, 133, 153, 163, 178] Dobson [106] has studied the geome- 
try of saddle node bifurcation and presented a method to avoid saddle node bifurcation [107] 
Overbye [182, 183] has also presented a method to restore power flow solvability 

On dynamic considerations, the problem of voltage instability has been studied by solving 
set of nonlinear differential algebraic equations [63,112,131,143] It has been shown that when 
the Jacobian associated with the algebraic system becomes singular, the system model breaks 
down and the load bus voltages can no longer be controlled [63,84,118,173] The eigenvalue 
analysis has also been carried out using dynamic models of power systems for analysis of voltage 
instability [61,87, 115, 125, 144, 160, 186, 190] 

Voltage stability analysis of an AC-DC system has been studied both as a steady state as 
well as dynamic phenomena In steady state analysis, the converters are assumed to act very 
fast Dynamic voltage stability analysis of a two terminal HVDC link and a AC system are 
presented in ref [5,12] The influence of the short-circuit ratio (SCR), DC line resistance and 
control parameters on the voltage stability characteristics of the systems is highlighted In 
the analysis, the authors have assumed constant current and constant extinction angle (CEA) 
control at the rectifier and inverter terminals respectively The converter control is assumed to 
be of proportional type Hammad et al [30,35] and Franken et al [70] have used the concept of 
voltage stability factor (VSF) for steady state voltage stability analyses The VSF was defined 
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as the ratio of the incremental change in converter bus voltage to the incremental change m 
capacitive susceptance at the same bus Effect of short circuit ratio (SCR) of both rectifier and 
inverter AC system on voltage stability of HVDC link has been investigated by Srivastava et 
al [96] 

With developments and increased understanding of the nonlinear dynamic system theory 
along with availability of inexpensive computing facility, it is being increasingly used to study 
the power system behaviour [161] Varaiya et al [117] presented a survey of nonlinear dynamical 
phenomena covering bifurcation and chaos They demonstrated that even the simplest represen- 
tation of a power system can exhibit both local and global bifurcations Saddle node bifurcation 
was hypothesised by some researchers as the only cause of voltage collapse [36, 54, 66, 101] 
However, later work demonstrated the possible role of other bifurcations such as Hopf bifur- 
cation, chaos and boundary crisis in power system voltage instability leading to voltage col- 
lapse [60, 76, 88, 122, 127, 139, 154, 191] Abed et al [26] demonstrated the effect of variable 
damping, frequency dependence of electrical torque, transmission line resistance and excita- 
tion control on the classical swing equation dynamics and observed the emergence of periodic 
orbit through Hopf bifurcation Alexander [33] demonstrated the existence of sub critical and 
supercritical Hopf bifurcations for different values of resistance to reactance (R/X) ratio of trans- 
mission lines A fourth order generator model along with excitation system have been considered 
by Rajagopalan et al [60] and it has been shown that a complex pair of eigenvalues associated 
with the excitation system results in Hopf bifurcation Ajjarapu et al [101] have considered 
a dynamic load model as suggested by Waive [40] along with generator swing dynamics and 
demonstrated multiple Hopf bifurcations and existence of a series of period doubling bifurcation 
leading to chaotic oscillations m a sample 3-bus system Hopf bifurcation has been demonstrated 
m synchronous motor dynamics [83]. Chen et al [47] demonstrated degenerate Hopf bifurcation 
m power system Frequency domain approach has been adopted by Kwatny et al [74] to study 
Hopf bifurcation in power system Kopell et al [17] have demonstrated the presence of chaotic 
invariant set m a 3 generator power system network described by swing equations using Mel- 
nikov’s method Nayfeh et al [76,88] investigated a single machine (quasi) infinite bus system 
in which periodic orbits became unstable through a series of period doubling bifurcations that 
led to chaos A single generator under excitation control may exhibit chaos as demonstrated by 
Salam [22] Salam et al [23, 29] demonstrated that in multimachine system even the classical 
swing equations could exhibit complex chaotic motion called Arnold diffusion Using the same 
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model of power system as in [101], chaotic dynamics have been studied using Lyapunov expo- 
nents by Chiang et al [127] and voltage collapse through boundary crisis was demonstrated by 
Abed et al [122] and by Wang et al. [191] Lee et al. [139] investigated m detail the period dou- 
bling route to chaos using the same model of power system Liu et al [177] proposed a technique 
to detect chaotic behaviour by computing largest Lyapunov exponents through measured values 
of state variables A new bifurcation term, singularity induced bifurcation (SIB) was introduced 
and applied to power system voltage stability by Venkatasubramaman [118] 

Voltage stability has also been called to be load stability [8] as the load characteristics have 
profound effect on the system dynamics [185] Dynamics of induction motor especially large 
induction motors have been reported to be crucial for voltage stability studies [99, 185] Effect of 
various load characteristics including synchronous and induction motor loads on voltage stability 
have been studied [55,56,58,62,79,87,114,125,150,160,170,172,184,188] Camzares [167] have 
demonstrated the effect of load model on the bifurcation characteristic of the system model 
The classical stability simulation of Swedish blackout on 27th December 1983 could not 
explain the system behaviour that led to cascade tripping This prompted the researchers to 
develop new load models for voltage stability studies such as by Hill [136], Waive [40], Xu et 
al [192] and Hill et al. [174] Probabilistic load models have also been considered m [43,149,164] 
The effect of on load tap changer (OLTC) is to restore the load However it has been found 
to aggravate the voltage instability [15,32,57,89, 119] Instability due to OLTC controls may 
lead to the possibility of either hunting or voltage collapse [44]. The behaviour of SVC has been 
studied by Hiskens et al [111] under voltage collapse 

Recently, the system stability has also been linked to its structural stability [185] The 
structurally unstable set form the bifurcation surfaces In order to eliminate bifurcation and 
chaos, Chiang et al [127] suggested that the system parameters should be restricted to enter 
the range where these phenomena would be present This, however, limits system loadabihty 
It has been demonstrated by Tan et al [154] that by varying damping the bifurcation and chaos 
can be controlled Wang et al [157] employed a feedback law using generator speed deviation 
as feedback signal to control stability and amplitude of bifurcated solution 


1.3 Objective of the Thesis 

From the survey of literature presented in section 1 2, it appears that most of the works on 
voltage stability analysis of AC-DC system [5, 12,35,70,96] have represented AC system by their 
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Thevenm equivalent. Though the concept of minimum singular value and condition number of 
power flow Jacobian has been successfully used m literature to assess static voltage stability of 
AC systems, it has yet to be applied to AC-DC systems considering the detailed representation 
of AC systems 

One of the recent operating concerns in the power system network is to maximise the voltage 
stability margin A limited research efforts have been made to achieve this objective through 
rescheduling of real and reactive power output of the sources Tiranuchit et al [51] have used 
maximisation of minimum singular value as the objective for optimal rescheduling Modern en- 
ergy management system (EMS) run the traditional optimal power flow (OPF) for rescheduling 
of sources in view of minimisation of the total fuel cost of generation or minimisation of the 
system transmission loss In all these formulations, fixed reactive power limits of generators have 
been considered Also, the effect of these optimal power flow schemes on the voltage stability 
margin has not been explored in the literature 

One important outcome of the studies reported in the literature is the general consensus 
that bifurcations in underlying mathematical models of a power system are closely linked with 
various types of instabilities A power system is usually designed to operate at a stable fixed point 
and a well designed power system should not exhibit any bifurcation under normal operating 
range of parameters Regions of the state space, where these bifurcations and the associated 
instabilities in behaviour occur, are reached only under disturbances or unanticipated loading 
conditions Nevertheless, the study of these bifurcations are important because the region of 
stable operating point is determined by the location of bifurcation The detailed model of 
induction motors with speed dependent generalized loading and speed governing loop along 
with excitation system dynamics has not been considered Also the dynamics of OLTC and 
SVS have not been investigated for bifurcation studies 

Recently new devices for the control of voltage and power have been developed [108,137,138] 
The prominent amongst these components are static var compensators (SVC), static phase angle 
regulator (PAR) and controllable series capacitors (CSC) The effectiveness of these FACTS 
(Flexible AC Transmission System) devices for controlling bifurcation and chaos has not yet 
been explored 

Therefore the motivation behind the work reported in this thesis was 

(l) To extend the concept of minimum singular value and condition number of load flow Ja- 
cobian to combined AC-DC systems considering the detailed repiesentation of AC system 
and also to AC systems having voltage dependent loads 
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(n) To study the effect of different generation rescheduling schemes on increasing the static 
voltage stability margin considering the voltage dependent capability of the synchronous 
machines and to explore a new computationally attractive scheme for generation reschedul- 
ing that maximises voltage stability margin especially under stressed condition 

(in) To observe various types of local bifurcations such as static and Hopf bifurcations and their 
effect on stability considering the dynamic models of generators and loads, to determine 
critical loading and control parameter values and to ascertain predominant states using 
participation analysis causing Hopf bifurcation 

(iv) To observe various types of global bifurcations and chaos considering OLTC, SVS, load and 
generator dynamics and possibility of voltage collapse due to these phenomena and also 
to explore other routes to chaos than the period doubling route m power system models 

(v) To study the effectiveness of FACTS devices m eliminating dynamic bifurcation and chaos 

1.4 Outline of the Thesis 

The present chapter introduces the problem of voltage stability, presents a brief survey of liter- 
ature in the area of voltage stability & collapse and lays down the motivation of the research 
work presented in this thesis 

In Chapter 2 the minimum singular value and condition number of the load flow Jacobian 
has been used to predict static voltage stability of integrated AC system considering voltage 
dependent characteristic of loads and extended them, for the first time, to combined AC-DC 
systems considering detailed representation of AC system The effect of control mode switching 
m the DC link and bus type switching in handling the reactive power limits of generators on 
static voltage stability has also been demonstrated 

In Chapter 3 four different schemes of generation rescheduling using objective functions as 
minimisation of total cost of generation, system transmission loss, maximisation of minimum 
singular value and a new scheme in which slack bus reactive power injection is minimised, have 
been used to maximise the static voltage stability margin Voltage dependent capabilities of 
synchronous machines have been considered in the model 

Chapter 4 investigates the local bifurcation namely static and Hopf bifurcations in sample 
power systems Effect of static bifurcation on power system stability has been demonstrated by 
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studying the associated unstable manifold of type-1 unstable equilibrium points on the stability 
boundary Center manifold theory has been used for studying the stability of Hopf bifurcation 
and determining critical values of parameters The effect of local bifurcations on the stability 
in parameter space has also been studied The hard limiters m excitation and speed governing 
control loop has been represented by a new approach using tanhyperbolic function 

Chapter 5 presents the system study results on observance of various types of global bifur- 
cations such as period doubling bifurcation, bifurcation of periodic orbit to torus, chaos via 
period doubling bifurcations and via torus breakdown and boundary crisis Lyapunov expo- 
nents and dimension of the attractor have been determined Broad band frequency spectrum 
has been observed under chaotic oscillations Occurrence of voltage collapse due to subcntical 
Hopf bifurcation and boundary crisis has been explored 

Chapter 6 demonstrates the application of FACTS devices such as CSC, PAR and SVC for 
control of bifurcation and chaos A first order delay model of these devices have been considered 
The control signals have been selected through participation analysis and the effect of varying 
controller gam parameters on elimination of dynamic bifurcations have been observed 

Chapter 7 concludes the main findings of the thesis and gives a few suggestions for further 
work in this area 



Chapter 2 

Static Voltage Stability Prediction of 
Integrated AC and AC/DC Systems 


2.1 Introduction 

Voltage stability, in most of the works, have been studied as static phenomena due to the reason 
that in many incidents of voltage collapse, it has been characterised as slow variation of voltage 
over a long time until it reaches close to the collapse point Moreover, the static consideration 
provides several conceptual observations about the phenomena 

Several researchers such as Venikov et al [6], Kwatny et al [36] and Sauer et al [78] have 
linked the singularity of power flow Jacobian to the steady state stability limit Wu et al [41] 
and Hiskens et al [56] demonstrated that the sufficient condition for an operating point to 
be small disturbance stable is the positive definiteness of the power flow Jacobian Kwatny et 
al [36] demonstrated that the static bifurcation of load flow equations are associated with either 
divergence type instability or loss of causality A generic assumption in considering the power 
flow equations as the static model for voltage stability studies is that there is some underlying 
dynamics, the fixed point of which is given by the power flow equations [106] In fact, the power 
flow solution is a subset of the fixed point of any generic dynamic model of the power system 
However, the power flow Jacobian is topologically diagonally dominant and almost symmetrical 
matrix This implies that the possibility of instability is associated with a real eigenvalue or 
singularity of the power flow Jacobian Assessment of singularity of power flow Jacobians by 
observing the divergence of power flow solution may sometimes be misleading as the solution 
may also diverge due to many other reasons In view of this, other indices such as minimum 


11 



CHAPTER 2 STATIC VOLTAGE STABILITY SYSTEMS 


12 


singular value and condition number have been used to determine the singularity of the power 
flow Jacobian 

Tiranuchit et al [51,52] presented an approach to predict the occurrence of voltage collapse 
using singular value decomposition of the load flow Jacobian and an efficient algorithm to com- 
pute the singular value using parallel processing Lof et al [113, 141] presented a fast method 
based on inverse iteration scheme to compute minimum singular value of the Jacobian suitable 
for large size systems Pai et al [58] extended the concept of singular value decomposition for 
detecting voltage collapse, considering an exponential type voltage dependent loads and Q-limits 
of the generators, and also proposed the condition number as a new indicator 

Voltage stability analysis of HVDC links has also been studied as static [35, 70, 96] and 
dynamic [5, 12] phenomena But m all these studies, the AC system has been represented by 
Thevenm equivalent 

From the literature survey, it appears that the minimum singular value and condition number 
of the load flow Jacobian have been effectively used for the prediction of point of static voltage 
instability in AC system and are yet to be applied to combined AC-DC systems Moreover, m 
all the voltage stability studies of AC-DC systems, the AC systems have been represented as 
Thevenm equivalent sources and detailed representation of combined AC-DC systems has not 
been considered Hence, the motivation behind the work presented in this Chapter is to extend 
the concept based on minimum singular value and condition number of the load flow Jacobian 
to combined AC-DC systems considering the detailed representation of the AC system and also 
to AC systems having polynomial type voltage dependent load characteristics 

The effect of considering reactive power limit of synchronous machines and control mode 
switching in DC link on occurrence of the voltage instability has also been investigated Studies 
have been conducted on IEEE- 14 bus AC system, and its three modified configurations to include 
a DC link, IEEE-30 bus AC system and a 156-bus Indian system having a DC link 


2.2 Singularity of Jacobian as a Condition of Voltage 
Collapse 

The singularity of the power flow Jacobian has been utilised to determine the point of voltage 
collapse To illustrate the validity of this concept, consider a 2-bus system as shown m Fig 2 1 
with an infinite source connected to bus-1 The real and reactive power loads (P&zQ) at bus-2 
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Figure 2 1 Sample AC system 


are given by the following equations 

P = ^smd (2 1) 

A 

V V 2 

q = _ cos *__ (22) 

The contours of P and Q in the (V, S) plane will look like those shown in Fig 2 2 The 
intersection of the contours of P and Q in the (V, 6 ) plane yields the steady state operating 
point, that is, the load flow solutions It is seen from Fig 2 2 that for a load Q = Qi and 

P = Pi, there are two solutions of voltages V\ and V \ Keeping Q = Q\ constant, as P is 

increased from P x to Pi, V\ decreases to V 3 and V 2 increases to V 4 If P is further increased to 
P 3 both the solutions coalesce at V 5 This point corresponds to static bifurcation (saddle-node) 
point At this point, the gradient vectors of P and Q will be colhnear, that is, 

VP = aVQ (2 3) 

where a is a scalar quantity The power flow Jacobian for the sample system in Fig 2 1 , 
assuming bus -1 as slack, will be 


ap dP 

as av 
aq ag_ 

d6 dV 


VP r 


Equations (2 3) and (2 4) imply that the load flow Jacobian will be rank deficient at P = P 3 
and Q = Qi which describe the voltage collapse point Thus, the singularity of Jacobian implies 
the point of static voltage instability 

This concept can be similarly extended to larger systems Assume a system having (n + 1) 
buses, (n + l ) tl1 being the slack bus The load flow equations can be written as, 


F(A,X,N) = 0 
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Figure 2 2 Loci of P and Q 


where A = [Pi, 

,Pn,Q 1, 

,Q n ]eR 2n 

£ 

II 

X 

) KiA? 

A] 6 R 2n 

II 

fa 

,F 2n ]eR 2n 


N = Network parameter such as [Y^us] 

The intersection of the contours of F in X € R 2n will yield a load flow solution for a given N 

If there exist at least two elements of F, say F % and F } such that either they do not intersect or 

they intersect at one point only, the gradient vector associated with them will be collmear, that 
•* 

is 


VP, = aVFj (2 6) 

Equation (2 6) denotes the linear dependence of two rows of the Jacobian, which implies that 
the load flow Jacobian is rank deficient or singular One approach to test the singularity of 
the Jacobian for increased loading or different operating conditions in the system is to check 
the divergence of the solution This approach, however, sometimes gives erroneous results as 
the load flow may also diverge due to numerical solution difficulties A more effective index to 
check the singularity of the Jacobian is through calculation of the minimum singular value or 
condition number as described below 
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2.2.1 Minimum Singular Value as a Measure of Proximity to Sin- 
gularity 

Consider a matrix A £ R nXn Let A A be the perturbation matrix such that (A + AA) is 
singular The notion of singularity is intimately related to the ability to compute its inverse If 
A is non singular, the following holds good 

(A + AA) -1 = (I + A^AA) -1 A -1 (2 7) 

The term (I + A -1 A A) can be shown to have an inverse if ||A -1 AA|| < 1, which is guaranteed 
if || AA|| < 1 /|| A -1 1| since ||A _1 AA|| < |j A - 1 1| || AA|| for some suitable matrix norm || || Hence, 
a measure of the nearness of matrix A to singularity is the number ||A _1 || _1 

Define ||A|| 2 = X max (A T A) = cr^ ax (A) where A max (A T A) maximum eigenvalue of A T A In 
terms of singular values, ||A _1 || = cr mai (A _1 ) The singular value decomposition is given by A = 
UEV T where U and V are orthonormal matrices (UU T = /, VV T = I) and S = diag[a x , , a n ] 
with cr max = a x ><r 2 > > <j n = a mtn Therefore, A" 1 = VE~ l U T and HA- 1 !!" 1 = a mtn {A) 

Hence, minimum singular value is a measure of proximity of a matrix A to singularity 

The minimum singular value of a power flow Jacobian matrix J can be computed by applying 
inverse iteration scheme [113, 141] The following are the steps applied for computation of 
minimum singular value 

(l) Perform LU factorisation of J 

(n) To start the iteration scheme, choose the right singular vector v such that ||u|| ^ 0 To 
perform computations on the power flow Jacobian matrix J for minimum singular value, 
the components of v corresponding to bus angles were set as zero, and those for bus 
voltages were taken as 1 This choice results in better convergence [141] 

(m) For k th iteration, solve J T u k = v k , for the left singular vector u k 

(iv) Compute minimum singular value <r k nxn = 

( V ) ^ Ikmm -*m7nll < to1 ’ st °P 

(vi) If convergence is not achieved in step (v), the right singular vector v k is updated for 
( k -f I)'' 1 iteration by solving Jv fc+1 = u k 

(vn) Compute the minimum singular value cr k im = 
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(vill) If ||<7*tn - °Ln\\ < io1 S to P ? else g° step (ill) 

This method converges rapidly when the ratio of cr n _i/cr n is much larger than 1 


2.2.2 Condition Number to Test Singularity of Jacobian 


The condition number of a matrix J is defined [19] as the ratio of its maximum singular value 
to minimum singular value 

r ma x( J) 


k(J) 


.( J ) 


(2 8 ) 


A small value of k(J) refers to a well conditioned matrix which is not very sensitive to per- 
turbations On the other hand, if a matrix is highly sensitive to perturbation, it is said to be 
ill-conditioned and is characterised by a very large value of k( J) The condition number k(J) 
can be used to measure the nearness to singularity of the power flow Jacobian When Jacobian 
is singular, cr mm will be zero and k(J) becomes infinite 

For computing fc(J), cr max is also needed m addition to cr mtn The cr max can be computed 
by applying inverse iteration scheme for finding the largest eigenvalue of the matrix using the 
following steps 


(0 

00 

(ill) 

(iv) 

(v) 

(Vl) 

(vn) 


Choose left singular vector u such that ||u|| 2 ^ 0 In this work, each component of u is 
taken as 1 

For k tfl iteration, find the right singular vector v k such that v k = J T u k 
Compute maximum singular value cr^ ar = jj^jjj 

If Ikmax - ^mai 1 1 < to1 , s t°P 

If convergence is not achieved in step (iv), the left singular vector u is updated for (k + l) th 
iteration by solving u fc+1 = Jv k 

Compute the maximum singular value cr^ax — 

If bLx - Vmal II < t0 1 sto P. else S° t0 Ste P ( u ) 


This method will converge rapidly when the ratio of <7i/<T2 is much larger than 1 
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2.3 Problem Formulation 


2.3.1 AC System with Voltage Dependent Load Model 

Voltage stability of a system is greatly affected by its load characteristics In power flow sim- 
ulations loads are often modelled as function of the bus voltage magnitude Usually frequency 
is assumed to be constant at its rated value Steady state composite load characteristics of 
equivalent consumption nodes are often approximated by the polynomials [168] 

Pl — a o + a \V + diV 2 (2 9) 

Ql — bo A b\V -f- (2 10) 

where Pl, Ql and V are expressed in per unit and ao, a\, a 2 , 6 0 , b\, 62 are the load coefficients By 
appropriately choosing coefficients in equations (2 9) and (2 10), a variety of load characteristics 
can be modelled The composition of loads may be different at each of the buses m a system 
However, a typical example of complex consumption is taken from ref [168] for the purpose of 
the present study as given below 

P L = 083 -0 3V + 0 47V 2 (2 11 ) 

£x, = 6 7— 15 3V + 9 6 V 2 (2 12) 

The above model has been utilised to represent loads at each of the buses in the AC system 
At nominal voltage (V = 1 pu), Pl = Ql — 1 pu Since the real and reactive complex 
consumption of a load node is the function of voltage at that bus, the load flow Jacobian 
elements corresponding to derivatives of P and Q with respect to V will be modified If H, L, 
M and N are the submatrices of the AC system load flow Jacobian corresponding to |^V, 
and |^V terms, the Newton- Raphson (N-R) equations with voltage dependent loads will 
be of the form 


' A P ' 


‘ H 

L' ' 


‘ A 6 " 

. A Q . 


M 

N' _ 


AV 

. V . 


(2 13 ) 

where only diagonal elements of the submatrices L' and N' are modified as follows (say corre- 
sponding to bus-p) 


V„ = ipp + (oWp + 2a,v‘) 

N'=N„ + (b l V t +2h, v !) 


(2 14 ) 
(2 15 ) 


The off-diagonal elements remain unaltered 
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T 


HVDC LINK 


1 flr 


a, 1 


PDp -hlQDp 


L — — j 


J 


PDq +jQDq 


PRp +jQRp 


-PIq +jQIq 


Figure 2 3 Sample HVDC link 

2.3.2 AC/DC System 

For steady state voltage stability analysis of an AC-DC system, the steady state performance 
equations of DC link can be integrated in the AC power flow equations 

Consider a monopolar HVDC link connected between two AC buses p and q, as shown m 
Fig 2 3 The performance equations can be written as 


VdR = Vdi + Rdh 

(2 16) 

V d R = k 2 a R V v cos a R - k 2 X c I d 

(2 17) 

V d i = k x aiV q cos 7 / - k 2 X c I d 

(2 18) 

Pr v - V dR Id 

(2 19) 

Piq = Vdi Id 

(2 20 ) 

Sr p = l-P/tp d” jQr p | = kk\Q,RVpld 

(2 21 ) 

Slq = \-Plq + jQlql = kk X aiVqI d 

(2 22 ) 

Qb, = {SI - 

(2 23) 

Qh = (Si, - Pi,)'' 2 

(2 24) 


where k , k x and k 2 are constants and their value are taken [18] as k = 0 995, k x = -hp, k 2 = | 
In the above equations, V p and V q are rectifier and inverter side AC bus voltages, V d R and 
V d i are the DC voltages on rectifier and inverter sides, a R & cn the rectifier and inveiter side 
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converting transformer tap ratio, an the firing angle of rectifier, 7 / the extinction angle of 
inverter, Pr p k Qr p the real and reactive power drawn by the rectifier from AC bus, Pj q k Qi q 
the real power injection into AC bus and reactive power drawn by the inverter from AC bus 
respectively and Id the DC line current Sr p = \Pr p + jQr p \ and Sj q = \—Pi q + jQi q \ are the 
apparent powers drawn by the converters Rd is the DC line resistance and X c the commutating 
reactance The following choice of base quantities have been made to define p u system for 
interfacing equations of DC link with equations for AC system 

For AC system Sr three-phase power base 
Vr line to-line voltage base 
Ib (current base) = 

For DC system Sg = S B , Vg c = V B , 7| c = \Z3I B 

The above performance equations (2 16) to (2 24) involve seven DC variables 

Xdc = tydR, Vdh h , ur, ai, <xr, 7 /] (2 25) 

Only three equations out of equations (2 16) to (2 24) are independent Hence, m order to have 
a unique solution, four variables should be prespecified In the normal operation of an HVDC 
link, one terminal controls the voltage and the other the current in the DC link Usually under 
normal operation, the inverter controls the voltage and the rectifier controls the current m the 
DC link When the AC voltage at the rectifier goes down, the rectifier may hit minimum ocr In 
such situation, the voltage is controlled by the rectifier and the current by the inverter When 
the inverter is controlling DC voltage, the possible combinations of variables to be prespecified, 
result into four possible control modes Ci to C 4 as given in Table 2 1 When the rectifier is 
controlling the voltage three control modes C 5 to C-j are possible depending on the combination 
of variables given in Table 2 1 In this formulation the tap ratios cir and aj are treated as 
continuous variables Equations (2 14) to (2 24), with any combination of the four prespecified 
variables as in Table 2 1 , will be a polynomial in V p and V q 

A DC link in an integrated AC-DC system may exist in any of the following configurations 

Case-A a DC link forming a part of an AC system with no AC line in parallel 
Case-B a DC link forming a part of an AC system with an AC line in parallel 
Case-C a DC link acting as an asynchronous tie between two AC systems 
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Table 2 1 Control Modes of HVDC Link 


Control Mode 

Specified Variables 

Ci 

Vdl, Pdl 

c 2 

or, 7/, Vdii Pdl 

C 3 

ocR,yi,ai, Pdi 

C 4 

a Ri a I, lli Pdl 

c 5 

0-R, &R, Gil Pdl 

C 6 

a R, lli a Ri Pdl 

C 7 

a Ri ah Vdl, Pdl 


The AC-DC power flow model used in the present study is due to Smed et al [95] Since the 
converters have been considered as loads on AC buses, the elements of the mismatch vector at 


p and q will be modified as follows 


AP p = AP"-P R „(K,,V„A-,fe) 

(2 26) 

AQ p = AQ;'-Q Kr (v r ,v,,x ic ) 

(2 27) 

AP, = AP“ C + P/,(V p , Vg, Xjc) 

(2 28) 

AQ, = AQ“ - Q,,(V P , V„ X dc ) 

(2 29) 


where A P ac and A Q ac are the real and reactive power mismatches without considering the DC 
link 

The modifications in the Newton- Raphson (N-R) load flow equations to include the DC link 
for the three possible configurations is described below 


Cases- A and B: Since the real and reactive power consumed by the converters are functions 
of only AC bus voltage magnitudes, the load flow Jacobian elements corresponding to the deriva- 
tives of P and Q with respect to V will be modified If H , L, M and N are the submatrices of 
the AC system load flow Jacobian corresponding to the |j, |£V, and terms, the N-R 
equations for the AC-DC system will be of the form 


' A P ' 


' H 

V ' 


i 

<1 

I 

.A Q 


. M 

N ' . 


L v J 


(2 30) 


where the modified elements of the submatrix L' will be 
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T ' — T | T/ ^ PRp 

^pp — v pp + v v Qy 

(2 31) 

V — T \ V ^ Prp 

L pq - L Pd + v g gy g 

(2 32) 

r - r v ^ Plq 

u qp ~ v t> qy v 

(2 33) 

V -L V dPlq 

U qq ~ V q QY g 

(2 34) 


The submatrix N' is also modified on similar lines The attractive feature of this formulation is 
that only four mismatch equations and upto eight elements of the Jacobian have to be modified 
and no new variables or equations are added to the solution vector when a DC link is considered 
m the power system network 


Case— C: Consider two AC systems (systems A and B) connected by a DC link The methodol- 
ogy described for Cases- A and B for modifying the N-R load flow equations can also be extended 
to this case Both the systems have to be handled simultaneously, that is, two slack buses, one 
in each system, are required to be considered If Ha, La , Ma , and Na represent the Jacobian 
submatnces of AC system A and Hg, Lg , Mg , and Ng for the system B, the modified N-R 
load flow equations for the two AC systems and the DC link combined will be 


AP a 

A P p 
AQ a 
AQ P 
A Pg 
AP q 
A Qg 

AQ, 


where N’ Aq = V a 



L'm = V, 


1 dV q 
9 Pr r 


9 dV a 


L' n „ = -V„ 


dPfo 


W' 


K 


3Qio 


p av a 


Ma 

A6 P 

A Va/V a 
A V p /V p 

ASg 

A S q 

A VgjVg 

Avyv, 


(2 35) 
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The elements of N' A , N' B , L\ and L' g will be similar to those defined for Cases-A and B For 
predictmg the voltage collapse, the singularity of the Jacobian, as defined m equation (2 30) 
for Cases-A and B and in equation (2 35) for Case-C of the combined AC-DC systems, can be 
studied 


2.4 Computational Procedure 

2.4.1 AC System 

Voltage collapse study for the real and reactive loads at all the buses were increased simultane- 
ously by a loading factor a (a > 1) till the load flow Jacobian approached to singularity The 
load flow was solved for each loading condition The minimum singular value and condition 
number were computed for each load flow solution In each iteration of the load flow, Q-limits 
of the generators were checked When a Q-limit of the generator bus violated, it was converted 
to a PQ type bus 

2.4.2 AC/DC System 

For voltage collapse study, the load flow of AC-DC systems was solved at different loading values 
increased gradually only at the inverter side AC bus and also m another case, simultaneously at 
all the buses The minimum singular value and the condition number were computed for each 
load flow solution To solve the load flow using the N-R method initial estimates of the DC 
variables are required The initial estimates for the DC variables Xj c were obtained from the 
specified value of Pi q corresponding to the rated DC voltage of the HVDC link and assuming a 
converter power factor of 0 9 and an inverter operating under minimum extinction angle 

For each of the AC-DC systems, the maximum and minimum limits of transformer taps on 
both inverter side (a/) and rectifier side (a B ) were taken as 0 9 pu and 1 1 pu respectively The 
firing angle was assumed to vary from a mtn = 5 deg to a max = 15 deg and extinction angle from 
7mm = 10 deg to 7 mai = 20 deg Resistance of DC line was considered as 6 006 x 10 -3 pu and 
commutating reactance as 0 0061 pu All AC bus voltages were set to 1 0 pu except at PV 
buses for which specified values were considered In each iteration of the load flow, Q-hnuts 
of the generators and operating limits of the DC variables were checked When Q-limit of a 
generator bus violated, it was converted to PQ type bus and m the case of limit violation of DC 
variables, they were fixed to their respective limiting values, which changes the control mode 
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of the operation accordingly. A scheme for checking the limits of the DC variables and control 
mode switching is given in the flow chart presented in Fig 2 4 The main steps for estimating 
the voltage collapse point are given in Fig 2 5 

2.5 Results and Discussions 

The system studies were conducted on HP-9000/735 computer system to test the potential of 
the proposed method to predict the point of voltage collapse of AC systems as well as AC-DC 
systems 

2.5.1 AC System 

The studies were conducted on three sample systems namely IEEE- 14 bus (Appendix-G), IEEE- 
30 bus (Appendix-I) and an Indian 156-bus systems ( Appendix- J) For each of the three 
systems, the base case loads were assumed to have the same composition at all the buses as 
given by equations (2 11) and (2 12) 156-bus UPSEB system has a DC link between 26 and 

113 buses The link has been modelled as a constant power source at bus-26 ( Pg — 10 0 pu, 

Qg = 0 0 pu) and a constant PQ load at bus-113 (P L = 1021 pu, Q L = 0 0) The following 
studies were carried out on all the three systems by simultaneously varying loads at all the buses 
by a loading factor a (a > 1) 

Case-1 Considering constant P, Q load model and without checking Q-limits of generators 
Case-2 Considering constant P, Q load model and the generators’ Q-limits 

Case-3 Considering voltage dependent P, Q load model and without checking Q-limits of generat 
Case-4 Considering voltage dependent P, Q load model and the generators’ Q-limits 

The variation of minimum singular value and condition number of load flow Jacobian with 
increase in loading have been observed for all the three systems The variation for the above 
four case studies of cr mtn for IEEE- 14 bus, 30-bus and 156-bus systems have been plotted in 
Figs 2 6, 2 7 and 2 8 respectively and the variation of the condition numbers are shown in Figs 
2 9, 2 10 and 2 11 respectively From these figures it is found that the cr mm reduces close to zero 
while the condition number k(J) increases sharply and goes to infinity near the point of voltage 
instability A sharp change in the value of cr mm and condition number has been observed at a 
loading when any of the generators in the system hits its reactive power limit This also reduces 
the maximum loadability of the system This can be clearly observed by comparing Figs 2 6 to 
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Figure 2 4 Flow chart for checking the limits of the DC variables 



Figure 2 5 


Flow chart for voltage collapse prediction 
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Figure 2 6 Plot of cr mln for IEEE- 14 bus system 

2 8 for Case-1 and Case-2 or for Case-3 and Case-4 m each system Further it has been found 
that with voltage dependent loads, the system maximum loadability increases (compare Case-1 
with Case-3 or Case-2 with Case-4) This is due to the fact that the voltage dependent loads 
reduce as the bus voltage decreases 

2.5.2 AC/DC System 

To test the potential of the proposed method to predict the voltage collapse of an AC-DC 
network, four sample AC-DC systems were considered 

System A an IEEE 14-bus system (Appendix-G) modified to replace line 20 (Fig G 1) by 
a DC link 

System B an IEEE 14-bus system with a DC link in parallel with line 20 (Fig G 1) 

System C an IEEE 14-bus system with the rectifier end connected to bus 14 of the IEEE- 14 
bus system and inverter end connected to bus 1 of a sample 2-bus system (Fig 
2 12 ) 

System D a 156-bus Indian system representing UP State power system network with a DC 
link (Fig J 1) 

The study was carried out on all the above four systems by keeping the DC power from inverter 
fixed in one case and by varying it in the other case. In addition to this, the following studies 
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Figure 2 7 Plot of cr mtn for IEEE-30 bus system 


Minimum SV 



Figure 2 8 Plot of <r mm for 156-bus Indian system 
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Figure 2 9 Plot of &(t/) for IEEE- 14 bus system 
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Figure 2 10 Plot of k(J) for IEEE-30 bus system 
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Figure 2 12 IEEE-14 bus extended system (System C) 
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0 135 0 572 1-010 1 6.6.8 1 886 2 323 

LOAD 

Figure 2 13 System A effect of Q-limit check and load variation on a min 
also performed on System A 

Variation of the load at the inverter bus without checking Q-limit 

Variation of the load at the inverter bus with generator Q-limit checked at each iteration of 

the load flow L 1A , . 

Simultaneous variation of the load at all the buses, considering the generators Q-limits 

he results of the studies conducted for system A are given m Figs 2 13 to 2 21, and those for 
stem B m Figs 2 22 to 2 24, for system C in Figs 2 25 to 2 27 and for system D in Figs 2 28 

30 repectively 

In system A, the load was increased at the inverter bus 14 from its base value of 0 135 pu 
. the minimum singular value became quite close to zero It is seen from Fig 2 13 that when 
generator Q-limit is checked, the maximum loading corresponding to the voltage collapse at 
inverter bus reduces Moreover, the generators are found to be constrained at their upper 
power limits Thus, the reduction in the maximum loading is due to a deficit in the 
power generations At a load of 1 99 pu when Q-limit of the generators are not checked 
1 01 pu when it is checked, a sudden discontinuity is observed in Fig 2 13 At this 
however, it is noted that the voltage control shifts to the rectifier from the inverter and 
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Figure 2 14 System A effect of Q-hrrut check and load variation on k(J) 



Figure 2 15 System A effect of Q-lnrut check and load variation on the bus voltage 
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Figure 2 16 System A effect of DC power and load variation on a mtn 
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Figure 2 17 System A effect of DC power and load variation on k(J ) 
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Figure 2 20 System A: effect of simultaneous load variation, on k(J) 
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Figure 2 21 System A effect of simultaneous load variation on the bus voltages 
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Figure 2 22 System B effect of DC power and load variation on cr mm 



Figure 2 23 System B effect of DC power and load variation on k( J) 
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Figure 2 28 System D effect of DC power and load variation on a min 



Figure 2 29 System D effect of DC power and load variation on h(J) 
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Figure 2 30 System D effect of DC power and load variation on the bus voltages 

the current control shifts to the inverter from the rectifier Discontinuities are also observed in 
Fig 2 13 when the reactive power outputs of the generators reach their limiting values During 
control mode shift (control mode C 3 to mode C 6 ), the elements of the Jacobian change and 
when a generator hits its Q-limit, one row and one column are added to the Jacobian due to 
switching of a PV type bus to a PQ type bus, and the minimum singular value decreases more 
rapidly From Fig 2 14, it is seen that the condition number is more sensitive and approaches 
infinity at the maximum loading point Fig 2 15 shows the variation of AC bus voltages at the 
rectifier and inverter ends The lowest bus voltage is observed at the inverter end in both the 
cases 

The study on system A was also performed by increasing the load at the inverter AC bus 
while the inverter end DC power was kept constant in one case and in the other case, was varied 
by the same amount as the load change It is seen from Fig 2 16 that a load of 0 944 pu can 
be supported when the DC power is kept fixed and 0 800 pu when it is also varied The reason 
for the higher maximum loading in the first case is due to the supply of the additional powei 
through line 13 which is connected to bus 7 having Q-support (shunt capacitors) Also in this 
case, the lowest voltage as observed from Fig 2 18 is at bus 14 when the inverter end DC powei 
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is varied by the same amount as the load The increase in the load is being met by the DC power 
transmitted through the DC line from bus 13. Since there is no Q-support at either bus 13 or 
at bus 12, which is connected to bus 13 through line 19, the voltage decrease at the rectifier end 
(bus 13) is faster than in the previous case. Fig 2 17 shows that the condition number increases 
to a very large value at the collapse point. 

In the next set of studies on system A, the loads at the load buses were increased by a factor 
a (o: > 1 ) a was varied from 1.0 and the collapse was observed at a = 1 755 (Fig 2 19) 
There are exactly four discontinuities corresponding to each of the four generators hitting their 
Q-limit Fig 2 20 shows a large value of the condition number at the collapse point and Fig 
2 21 indicates that all the bus voltages (except the slack) approach the nose of the P-V curve 
simultaneously 

The study was conducted on system B with the inverter end DC power fixed m one case and 
varying (by the same amount as the load at the inverter bus) m the other case Fig 2 22 shows 
that the minimum singular value assumes a sufficiently low value for the maximum loading of 

1 065 pu and 1 028 pu when the DC power is kept fixed and varying respectively It can be 
noted that the maximum loading has increased in system B as compared to system A due to 
additional power support provided by the AC link in parallel with the DC link m system B The 
condition number approaches infinity at the maximum loading, as seen from Fig 2 23 The 
lowest voltage is seen at the inverter bus when the DC power is varied, for the same reasons 
mentioned for system A (Fig 2 24 In a similar set of studies on system C, it is seen from Fig 

2 27 that if the DC power is kept fixed the voltage collapse takes place m the 2 -bus AC system 
and the 14-bus AC system remains almost unaffected However, the discontinuity m the curve 
is due to control mode change (from mode C 3 to mode C&) But when the DC power is also 
varied by the same amount as the load at the inverter bus, that is, assuming that the increase 
in load is met by the power flow in the DC link, the collapse of bus voltages takes place in the 
14-bus AC system as the power flow in the bus system remains almost unaffected. Fig 2 25 
shows that the minim um singular value goes sufficiently close to zero at the maximum loading 
and Fig 2 26 indicates that the condition number is sufficiently large at this point 

Studies were carried out on the 156-bus Indian system (system D) for the conditions studied 
in systems B and C In this system a DC link (see Fig J 1 ) exists between buses RHND(S) 
and NCR for transportation of 1000 MW of power from RHND(S) to NCR (Dadn) bus NCR 
has a real power load of 735 MW (7.35 pu) The load flow for the 156-bus Indian system was 
run without considering the DC link and it was found that the load flow did not converge The 
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convergence of the load flow was possible only when the DC link was considered From Fig 
2 28, it is seen that the maximum loading at the inverter bus is 9 987 pu when the increase in 
the load is met by the DC power, and when the DC power is kept fixed, the maximum load 
at which the collapse occurs is 9 565 pu From Fig 2 29, it is seen that the condition number 
is very large at the collapse point Fig 2 30 shows that the lowest voltage is observed at the 
rectifier end when the DC power is varied and at the inverter end when DC power is kept fixed 
The variation of voltage magnitude at the rectifier end is only about 5% when the DC power 
is varied This confirms the observation in ref [40] that the drop in bus voltage magnitude can 
not be taken as a voltage collapse indicator 


2.6 Conclusion 

From the works presented in this chapter, the following conclusions are drawn 

(l) The concept of minimum singular value and condition number of load flow Jacobian can 
be used to predict the static voltage stability of an integrated AC-DC network as well as 
the AC system having a polynomial type voltage dependent loads 

(n) At maximum loadabihty point, the condition number approaches to infinity and the mini- 
mum singular value approaches to zero A sharp change in the values of condition number 
and minimum singular value of Jacobian is observed when Q-lmnts of generators are vio- 
lated or control mode switching in a DC link takes place 

(m) Considering the static model, it has been found that if the load is varied at only one bus 
in an integrated AC-DC system, the collapse of voltage takes at the bus at which the load 
is being increased However, with simultaneous changes in the loadings at all the buses, 
all the bus voltages approach the nose point of the P-V curve simultaneously. 

(iv) Consideration of reactive power output limits of synchronous machines reduces the loading 
at which the system becomes voltage unstable Thus it reduces the maximum loadabihty 
of the system 

(v) The voltage dependent characteristic of loads enhances the voltage stability of the system 



Chapter 3 


Effect of Generation Rescheduling on 
Voltage Stability Margin 


3.1 Introduction 

Maintaining the required voltage stability margin has become one of the recent concerns of 
power utilities Extensive works exist in the area of steady state voltage stability and margin 
predictions [168, 169, 189] Some of the indices suggested for the static voltage stability are 
reactive power sensitivity to bus voltage [69,134], stability margin in terms of real or reactive 
power [116,181], VCPI (Voltage Collapse Proximity Indicator) derived from the sensitivity of 
total reactive power generation to real or reactive load increase [28,102], VIPI (Voltage Instability 
Proximity Indicator) based on the closeness of multiple load flow solutions [24,68], minimum 
singular value, eigenvalue or condition number of power flow Jacobian [51,52,58, 113, 141, 142, 
179] 

Tiranuchit et al [51] proposed an optimal generation scheduling scheme to maximise the 
minimum singular value of the Jacobian using a linear programming technique Generators 
and condensers were modelled by their fixed reactive power limits and their voltage dependent 
capabilities were not considered Begovic et al [81] demonstrated that by using optimal power 
flow strategy to minimise the weighted sum of the absolute values of the control actions, the 
voltage stability can be improved 

Traditionally generation scheduling is performed in view of the minimisation of total fuel 
cost of generation or system transmission loss These schemes form a part of modern energy 
management system Ref [85] presents a review of various optimal power flow strategies being 
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used m practice However, effect of these conventional optimal power flow schemes on the voltage 
stability margin has not been reported in the literature Hence, in this Chapter, the following 
0 P F schemes for generation rescheduling have been considered and their relative impact on 
voltage stability margin has been studied 

(1) Minimisation of the total fuel cost of the generation 

(2) Minimisation of the total system transmission loss 

(3) Maximisation of the minimum singular value of the Jacobian 

(4) A new scheme which minimises the reactive power injection at the slack bus 

The voltage stability margin has been compared for the above four formulations m terms of 
minimum singular value of power flow Jacobian as well as the total system reactive power 
margin in each case for IEEE- 14 bus, IEEE-30 bus and 156-bus Indian System The reactive 
power margin has been computed by estimating the closest distance of an operating point to 
the static bifurcation set in reactive parameter space [132] 

In addition, the voltage dependent real and reactive power capability of the synchronous 
machines (generators and condensers) have been considered 


3.2 Generator and Synchronous Condenser Model 


One of the major factors responsible for the long term voltage instability is the heating of 
armature and field circuit of synchronous machines [189] This underlines the need for their 
appropjiate modelling The synchronous machines for the present study have been modelled as 
constant voltage source behind the synchronous reactance The real and reactive power outputs 
of a salient pole generator ignoring armature resistance is given as follows [2] 


o V ‘ E * li 1 ? 

P 0 = _ slni + T 

Qa = j^-cos6-V? 


x„ 


1 

X d 


cos 2 6 sin 2 5 


X„ 


(3 1) 


(3 2) 


where E q is the induced emf, Vj the terminal voltage, 6 is the phase angle between V t and E q 
Xd-, X q are the d-axis and q-axis synchronous reactances 

The corresponding equations for power angle characteristic of a cylindrical rotor machine 
having synchronous reactance X a can be written as 
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V<E 

Pa = (3 3) 

V*E V 2 

Qg = \r~^~ cos 8 (3 4) 

The machine real and reactive power outputs (Pa and Qq) and the other variables such as E q 
and 8 are physically constrained because of limits on maximum electrical and mechanical stress 
and temperature rise to which the machine can be subjected to 

The maximum real power output PG max °f the generator is usually constrained by turbine 
rating Also, it is not desirable to reduce the real power output of generator m certain power 
plants below a minimum value PG mxn unless it is taken off-line For example, m a fossil fuel 
plant, constant boiler temperature must be maintained to prevent slagging These operating 
limits on the real power output of generators can be expressed as 


p G m>n < Pg < Pg. 


(3 5) 


The synchronous condensers are connected in a network to meet the lagging and leading VAR 
requirement Neglecting all losses, their real power output and the load angle 8 can be assumed 
to be zero The reactive power output of a synchronous condenser is, thus, given by 


Qg 


V t E q V* 


(3 6) 


3.2.1 Rotor Heating Limit 

If the effect of the saturation is ignored, the field current of a synchronous machine (both 
generator and condenser) is proportional to the induced emf E q Therefore, the maximum 
allowable rotor current can be mapped to a maximum allowable induced emf E qmax Thus, the 
rotor heating limit can be expressed as, 

E q < E qmax (3 7) 


3.2.2 Underexcitation Limit 

Stability considerations impose a limit on the reactive power that a generator can absorb while 
operating in leading power factor mode Under AVR action, the leading power factor operation 
causes increase in the load angle Hence, it is possible to map the underexcitation limit of a 
generator to a suitable maximum limit on the load angle 
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6 < S m ax (3 8) 

For stable synchronous operation of the condenser with the rest of the network, synchronising 
power must be greater than or equal to zero Under these conditions, if the excitation voltage 
is allowed to be negative, the maximum reactive power that can be absorbed by the condenser 
is given as follows [2] 

V? 

Q mtn = -T7- (3 9) 

A q 


3.2.3 Stator Heating Limit 


The real and reactive power output of a generator is also limited by maximum allowable armature 
current „ 

'-“max 


Ia = 


\J{Pg + QV) 

v t 


< L 


d-max 


(3 10) 


The reactive power output of a synchronous condenser is also limited by maximum allowable 
armature current 


Qa < (3 11) 

The equations (3 1) and (3 2) for the salient pole or equations (3 3) and (3 4) for the cylindrical 
rotor together with equations (3 5), (3 7), (3 8) and (3 10) represent the generator operating 
capability under steady state Though it is possible to overload the generators [9], the overload 
capability of the generators has not been considered in this thesis 

Since the terminal voltages of synchronous machines are not usually allowed to exceed a 
specified maximum value V™ ax , the following additional constraints will apply 

V t , < V™ ax (3 12) 


3.3 Generation Rescheduling Schemes 

The voltage dependent model of synchronous generators and condensers, as discussed in section 
3 2, can be integrated in the power flow programme The nonlinear real and reactive power 
balance equations at any bus i can be written as follows 

(3 13) 


P,(»,V)=Pg,(^E,)-Pl, 
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Q t (0,V) = Q Gl (6,E q )-Q Li (3 14) 

where P, and Q, refer to the real and reactive power injections at bus-i, Pi, and Qi, the loads, V 
and 6 the bus voltage magnitude and angle vectors Since in this formulation, a large number of 
variables are constrained, the load flow for determining the base case operating point has been 
formulated as an optimisation problem The following objective function has been minimised 
to obtain the load flow solution subject to equality constraints (power balance) equations (3 13) 
and (3 14) and inequality constraints as given in equations (3 1) to (3 12) 

Minimise £ Wv,(V,,-V~^y (3 15) 

where Wv, is a weighing factor and N pv , the set of PV buses in the network 

Four different formulations of the generation rescheduling have been attempted in the present 
work considering voltage magnitude and angle at all the system buses and the internal nodes 
of the generators and the synchronous condensers as optimisation variables The objective 
functions considered for the generation rescheduling schemes are described below 

3.3.1 Minimisation of the Total Fuel Cost (OBJ-I) 

The cost characteristics of a generator-i is a function of its real power output and can be 
approximated by a quadratic function as below 

C(P a ,) = ^,Pg. 2 + b t P G , + c, (3 16) 

The total fuel cost of the generation for a system can be given as 

Ct = £ C(P c .) (3 17) 

teNg 

where N g is the set of thermal generating plants in the system In the present formulation, 
scheduling of hydrogenerators have not been considered 

3.3.2 Minimisation of the Total Transmission Loss (OBJ-II) 

The total real power transmission loss of the system is taken as the sum of the real power 
injections P, n] , at all the buses Thus, for a N bus system, the real power system transmission 
loss P/ 034 can be given as, 

^ = ( 318 ) 

t=i 
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P tnJi can be expressed m terms of bus voltages and angles as follows 

N 

Ptnj, = VtjVtV] cos (0, -6j- 7y) (3 19) 

J=1 

where y tJ L is the ij th element of U, U3 


3.3.3 Maximisation of Minimum Singular Value (OBJ-III) 

The power balance equations (3 13) and (3 14) can be linearised around an operating point to 
yield the following Newton-Raphson power flow relationship 


" A P ' 


' Pe 

Pv ' 

‘ A 9 ~ 

A Q 


_ Qe 

Qv 

AV 

- V . 


s. v ✓ 

= J 


The Jacobian matrix J is composed of the four submatrices Pg, Pv, Qe and Qv The system 
voltages are infinitely sensitive to bus injections if the Jacobian is singular l e when det { J} = 0 
Using Schur complement [141] and assuming Pg to be non-singular, it can be shown that 


det { J} = det {P e } det {Q v - Q g P g 1 P V } 

' v ' 

-Jr 


(3 21) 


If both the active power injections, A P and the changes m the bus voltage angle A# m equation 
(3 20) are set to zero, the relationship between changes m reactive power injections and voltage 
magnitudes can be described by the submatrix Qv and written as 


[A<2] = [Qv] 


rAUi 

. V . 


(3 22) 


However, there also exists considerable coupling between change in reactive power injections 
and change m bus voltage angles especially during stressed system conditions [189] One way of 
considering it is to consider the full Jacobian matrix J From equation (3 21), it is seen that in 
order for the Jacobian J to be singular, the reduced matrix Jr must be singular if the submatrix 
Pg is non-smgular 

Upon setting A P = 0 in the equation (3 20), the change in reactive power injections can be 
related to change in bus voltage magnitudes as follows 


[A Q]=[Qv-Q e Pf'Pv 


rAU] 

. V . 


(3 23) 
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Minimum SV 



1 1 002 1 004 1 006 1 008 1 01 1 012 1 014 1 016 1 018 

Load Parameter 


Figure 3 1 Plot of cr mtn of J, Jr, Pg and Qy for increase in system loading 

In order to observe whether the matrices J, Jr , Qy and Pg become singular simultaneously 
for increase m loading, their minimum singular values were plotted as a function of loading 
parameter a (a > 1) for 156-bus Indian System as shown m Fig 3 1 The real and reactive 
loads at all the buses were increased simultaneously by the factor a, retaining the base case 
power factor of loads at all the buses It was observed that out of the four matrices, the 
minimum singular value of J and Jr approached zero corresponding to a = 1 0175 Only these 
two matrices become singular when the loading m the system was increased, while the matrices 
Pg and Qy did not become singular at the loading corresponding to the static bifurcation set 
Similar observation has also been noted by Lof et al [141] In the present formulation, minimum 
singular value of matrix Jr has been maximised 

Though the matrix J is highly sparse, the matrix Jr is dense The following algorithm [1 13] 
can be used for sparse formulation to compute a mtn of Jr 

(i) Find LU factorisation of J 

(n) To start the iteration scheme, choose the right singular vector v such that ||u|| /0 To 
perform computations, the components of v corresponding to bus angles were set as zero, 
and those for bus voltages were taken as 1 This choice result in better convergence [141] 
Set k = 1 
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(m) For k th iteration, solve J T u k = v k , for the left singular vector u k 

(iv) Set u t = 0 for all i corresponding to the active power injections 

(v) Compute a mtn (J R ) = jjjjjj£ 

( V1 ) If \Wmin - ^mmll < tol. Stop 

(vn) If convergence is not achieved in step (vi), the right singular vector v k is updated for 
( k + l) </l iteration by solving Jv k+1 = u k 

(vm) Set v t = 0 for all i corresponding to the node angles 

(ix) Compute a min (J R ) = 

( X ) If Ikmtn ~ a Ln\\ < to1 sto P, else g° to ste P 0») 

3.3.4 Minimisation of Slack Bus Reactive Injection (OBJ-IV) 

Tiranucb.it et al [51] observed that the effect of maximising the minimum singular value of power 
flow Jacobian is to reduce the slack bus reactive power injection and call upon the other reactive 
sources to inject more reactive power This can be achieved by a simpler formulation utilising 
a new objective function as given below in that the slack reactive injection is also minimised 
along with the load flow objective 

Minimise Y W v , (v t , - + W Q Q tn , 3 (3 24) 

*GAT pt , 

where Wq is weighing factor associated with the slack bus reactive power injection Q t n Js If 
subscript s stands for slack bus, can be expressed as 

Qtn U = Yu y»3 V » V J S1I1 (^ “ ~ 7 ( 3 2 5 ) 

J = 1 

3.4 Estimation of System Reactive Power Margin 

For a given reactive loading vector Q° L , assume that the resulting operating state is stable 
When Qi is varied, the operating state may change For some Ql = Q* L , let the load flow 
Jacobian is singular In reactive parameter space set of all Ql, such that the Jacobian is 
singular, define a boundary of the feasible region say Cl for stable operation Cl typically consists 
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of hypersurfaces and their intersections The saddle node bifurcation instability, when the 
reactive loads encounters Cl, can result in catastrophic collapse of the system voltages and 
blackouts [54, 66] One of the ways to monitor the proximity of an operating point to voltage 
collapse is to measure the closest distance of the operating point Q° L to Cl Although the geometry 
of 0 is not known, fl has been assumed to be convex m the present study In particular, it is 
useful to calculate a critical reactive load power Q* L m fi for which ||Q£ — lb 1S a local 
minimum [13,97, 107, 132,133,163] The line segment Q° L to Q* L represents a worst case reacti\e 
load variation and || Q* L — II2 measures the proximity of Q° L to Cl This problem has been 
formulated as an optimisation problem whereas it was solved by using Newton-Raphson method 
in ref [132, 163] A sequential quadratic programming algorithm has been used to find Q* L such 
that ||Q* — Q ° l || 2 is a minimum while satisfying the power balance equations (3 13) and (3 14), 
implying that the load flow Jacobian J in equation (3 20) is singular 


3.5 Results and Discussions 

The studies were conducted on the following three sample AC system using HP 9000/735 com- 
puter system 

System-A IEEE- 14 bus system given in Appendix-G 
System-B IEEE-30 bus system given in Appendix-I 
System-C 156-bus Indian system given in Appendix-J 

Sequential quadratic programming algorithm was used for the optimisation in each of the four 
generation scheduling formulations 

The numerical values of Xd, X q and E qmax for the generators and synchronous condensers 
were taken from [7] in accordance with their capacity The maximum allowable load angle 6 max 
was taken to be 70 deg for all the generators m the three systems The studies were carried 
out by varying real and reactive loads at all the buses simultaneously by a factor ex (a > 1) 
The weighing factors WV, m equations (3 15) and (3 24) were taken as 1 0 in each of the system 
studies and Wq in equation (3 24) was selected as 0 1 The value of Wq was decided through 
experimentation and was taken as the one resulting into the maximum voltage stability maigin 
(VSM) The study results of the three systems are given below 
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Table 3 1 Voltage Stability Margin for System-A 


Loading Factor 

VSM 

LFLOW 

OBJ-I 

OBJ-II 

OBJ-III 

OBJ-IV 

a= 1 0 

Q-margm (Zi) 

1 081325 

1 074358 

0 998987 

1 080723 

1 078052 

Q-margm (I 2 ) 

0 371041 

0 368844 

6 341533 

0 370492 

0 370179 

&mtn 

0 6008 

6.5984 

0 5914 

0 6035 

0 6014 

O 

II 

s 

Q-margin (Lx) 

0 332611 

0 349178 

0 454444 

0 539591 

0 538171 

Q-margm (I 2 ) 

0 105002 

0 112257 

6 146889 

0 173682 

0 173885 

&mtn 

0 2810 

0.2478 

0 2952 

0 3399 

0 3345 


3.5.1 System— A 

There are 3 generators and 2 condensers in this system. Generator # 1 at bus-1 was taken as 
the slack Two sets of results were obtained corresponding to the base case loadings (loading 
factor a = 1 0) and the stressed conditions at higher loading (a = 2 0) Loads at all the buses 
were raised to two times their base values in the stressed condition The reactive power margins 
(both L\ and L 2 norm) using the method described in section 3 4 and the minimum singular 
value a mtn of final Jacobian submatrix Jr obtained m each of the cases were computed and are 
given in Table 3 1 It is seen from the Table 3 1 that at base case loading amongst the four 
objectives, OBJ-III yields a maximum of a mtn and reactive power margin (both L\ and L 2 norm) 
which are very close to the results obtained from OBJ-IV and OBJ-I However, none of them 
could improve the reactive power margin available at base case condition though the minimum 
singular value has increased slightly with OBJ-III and OBJ-IV At higher loading, though OBJ- 
III resulted in maximum cr mm , but provides (slightly) less reactive power margin (both L\ and 
L 2 norm) in comparison to OBJ-IV The voltage stability margin, however, has substantially 
increased by using either of the objectives from the unoptimised case The optimum real and 
reactive power outputs were obtained using the four objectives and their values are given in 
Table 3 2 along with load flow values at a = 2 0 

It is found that OBJ-III which maximises the minimum singular value does not yield a 
maximum of reactive power margin, for example at a loading (a = 1 0) much lower than critical 
loading The correlation between minimum singular value and reactive power margin is lost 
because at lower loading, a singular value which is minimum may not remain so when loading 
is increased 
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Table 3 2 Real and Reactive Power Outputs (in p u ) for System-A (a = 2 0) 


Pg&Qg 

LFLOW 

OBJ-I 

OBJ-II 

OBJ-III 

OBJ-IV 


5 3130 

4 1762 

4 3112 

4 8080 

4 7510 

Pg 2 

0 4000 

0 8466 

0 7143 

0 4636 

0 4168 

Pg 3 

0 2000 

0 7287 

0 7096 

0 5124 

0 6027 

Qg i 

1 4288 

1 4332 

1 0497 

0 8511 

0 8093 

QG 2 

0 8767 

0 4658 

0 6740 

0 8729 

0 8983 

Qg 3 

0 9826 

0 6402 

0 7121 

0 9095 

0 8451 

Qg 4 

0 4188 

0 4224 

0 4346 

0 4400 

0 4412 

Qgs 

0 3982 

0 3950 

0 4071 

0 4163 

0 4159 


3.5.2 System-B 

There are 3 generators and 3 condensers m this system Generator # 1 at bus-1 was taken as 
the slack Two sets of studies were performed One corresponding to normal (base case) loading 
at a = 1 0 and the other for a higher loading corresponding to a = 1 6 Value of a higher than 
1 6 could not be considered as it resulted into divergence of the load flow Table 3 3 presents the 
values of voltage stability margin (Q-margin and <r mtn ) obtained through load flow and with the 
different objectives described in the earlier section At normal loading, though OBJ-III results 
in a maximum of cr mtn , the proposed new OBJ-IV yields the maximum reactive margin Also 
at higher loading condition, the OBJ-IV turns out to yield a reactive margin which is almost 
same as with OBJ-III At normal (low) loading condition, the generation rescheduling does not 
help in increasing the voltage stability margin However, distinct improvement in the margin is 
observed during stressed conditions The generator outputs obtained with the four objectives 
along with load flow results at a = 1 6 is presented in Table 3 4 

3.5.3 System-C 

The data for this system has been slightly changed from one given in Appendix- J At bus- 
24 a synchronous condenser has been assumed to be present, and generators are placed at 
bus-1 to 23 The loading in this system already corresponds to the stressed condition Since 
the generator cost coefficients were not available for this system, the studies were confined to 
generation rescheduling using only OBJ-II, OBJ-III and OBJ-IV From the studies conducted 
on System-A and System-B (Table 3 1 and Table 3 3), it was found that minimum singular 

CENTRAL LIBRARY 
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Table 3 3 Voltage Stability Margin for System-B 


Loading Factor 

VSM 

LFLOW 

OBJ-I 

OBJ-II 

OBJ-III 

OBJ-IV 

<2=10 

Q-margin L\ 

0 710141 

0 723411 

0 725283 

0 727966 

0 716998 

Q- margin Li 

0 194651 

0 184339 

0 186435 

0 188566 

0 193180 


0 2171 

0 2081 

0 2130 

0 2191 

0 2175 

a = 1 6 

. 

Q-margin Li 

0 314804 

0 338730 

0 336078 

0 339490 

0 339442 

Q-margin Li 

0 071728 

0“ 077342 

0 076710 

0 077504 

0 077496 


0.1242 

0 1277 

0 1275 

0 1285 

0 1285 


Table 3 4 


Real and Reactive Power Outputs (m p u ) for System-B (a = 2 0) 


PakQa 

LFLOW 

OBJ-I 

OBJ-II 

OBJ-III 

OBJ-IV 

Pqi 

4 4310 

2 2250 

2 0087 

2 0648 

2 1868 

Pg 2 

0 4000 

1 3618 

1 4230 

1 3414 

1 2253 

Pg 3 

0 2000“ 

1 2061 

1 3473 

1 3739 

1 3736 

Qg i 

-0 1896 

0 0441 

0 1218 

0 0297 

0 0018 

Qg 2 

1 5745 

0 9642 

0.8660 

0 9981 

1 0433 

Qg 3 

1 3142 

0 8807 

0 8552 

08119 

0 8123 

Qg< 

0 2545 

0 2565 

0 2564 

0 2568 

0 2568 

Qg s 

0 3792 

0 3817 

0 3812 

0 3825 

0 3824 

Qg 6 

0 2519 

0 2542 

0 2541 

0 2546 

0 2546 
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Table 3.5: Real and Reactive Power Outputs (in p.u.) for Systera-C (a = 1) 


IH8S9 

Real Power Output 

Reactive Power Output 

LFLOW 

OBJ-II 

OBJ-III 

OBJ-IV 

LFLOW 

OBJ-II 



i 

7.6781 

■iartai 

4.2247 

3.3033 


0.8752 

0.7374 

0.6205 

2 



0.5250 

0.5250 



0.0493 

0.0729 

3 



3.2493 


1.1045 

WMMMM 

■BEWiFl 

0.4343 

4 


2.5666 

2.2618 

i»fi| 

0.7143 

0.5152 

■mm 

0.6575 

5 


0.2500 

wmM 

M%?*S1 

■O 

0.0321 

0.0393 

0.0656 

6 

1.0000 

2.5000 

■HI 


MBS 

0.2644 

0.2431 

Hill 

7 

2.5000 

0.5407 

0.7555 

■TOfgEl 

TiIMil 

BSSfl 


WEMSM 

Q|[| 

2.3000 

2.4673 

3.3752 

3.4925 

0.6544 

— 1 


BilBi 

9 

1.4000 

2.1623 

Mtiifitol 


ISES! 

0.0935 

0.3395 

wmmwm 

10 

1.8000 

iiliMttil 

■ 

I 

■Ell 

MgM$| 

0.0487 

0.0357 

— ■ 

11 

HIa’il'8'111 

mom 

— BMI 

8hhhmii 

0.1152 

0.1018 

0.0442 


12 

HU 

aBalMifcj 
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value cr m i n can act as a reliable indicator of voltage stability margin during stressed condition 
of the network and it gives similar information as the Q-margin. However, at light loading 
condition (away from saddle node bifurcation points), the correlation between Q-margin and 
cr min is lost. Since the base case loading of the System-C corresponds to a stressed condition, 
effect on voltage stability margin for this system has been analysed in terms of <r m , n only. The 
source real and reactive power outputs (in p.u.) obtained from the three rescheduling schemes 
and their base values are presented in Table 3.5. 

Fig. 3.2 presents the variation of cr mjn in each case as the loading is increased in the system. 
Loading has been increased upto a = 1.38 beyond which the load flow did not converge. From 
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Figure 3 2 Effect of load variation on cr min for 156-bus Indian System 

this figure, it can be observed that at all the loadings, the voltage stability margin is enhanced 

using either of the three (OBJ-II, OBJ-III or OBJ-IV) schemes from the unoptimised case 

3.6 Conclusion 

The studies presented in this Chapter reveal the following 

(i) The generation rescheduling using any of the four objectives in stressed condition of system 
increases the voltage stability margin which can be measured either in terms of system 
Q- margin or the minimum singular value of the power flow Jacobian 

(n) At light load conditions, the rescheduling of generation using any of the schemes is unable 
to increase the Q-margm cr mm , however, is slightly increased with OBJ-III or OBJ-IV In 
such cases, any correlation between cr mtn of the power flow Jacobian and Q-margin is lost 

(in) Amongst the four objectives tried out, the new OBJ-IV has consistently resulted into 
better results than other objectives which is close to that obtained with OBJ-III However, 
computationally it is less expensive in comparison to OBJ-III Hence the new objective 
(OBJ-IV) of minimising the slack bus reactive injection along with load flow objective can 
be effectively used for voltage stability margin enhancement 




Chapter 4 


Observance of Local Bifurcations 


4.1 Introduction 

An electric power system is a complex nonlinear dynamical system and its behaviour can be 
described by nonlinear algebraic and differential equations The behaviour of a power system 
may undergo qualitative changes under system disturbances and parameter variations In dy- 
namical systems’ lingo, this qualitative change m the phase portrait of a system, when a system 
parameter is quasistatically varied is known as bifurcation The parameter under variation is 
called bifurcation •parameter It is found in the literature [117, 161] that bifurcations m underly- 
ing dynamical model of a power system are closely linked with its instability This is especially 
true for voltage collapse precipitated by slow variation in a system parameter such as load. This 
chapter deals with the local bifurcations namely static and Hopf bifurcation A local bifurcation 
means bifurcations occurring in the neighbourhood of the fixed operating point 1 

The region of attraction of stable operating point (s e p ) is decided by unstable operating 
points. (u ep ) and closed orbits (especially the type-1) lying on the stability boundary of the 
stable operating point Following the instability, the unstable trajectory is determined by the 
unstable manifold of the controlling unstable operating point [128] 

In this Chapter, the study of the unstable manifold of the unstable operating point has 
been conducted to differentiate between angle and voltage stability The studies have been 
conducted on the two sample systems, one being a 3-bus system and the other a 9-bus system 
having multiple generators Further, Hopf bifurcation theory has been used to determine the 

1 Fixed operating point, fixed point, equilibrium point and singular point have been used interchangeably in 
this thesis It simply means a point in the state space where the velocity vector vanishes l e x = 0 
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critical values of control system and load parameters A detailed nonlinear dynamical model 
of generator along with its excitation control system and turbine governing system loop and 
induction motor load have been considered The hard limiters m the excitation control loop 
and turbine speed governing loop have been modelled using a new approach in an analytic form 
and the load torque as nonlinear function of speed The analysis have been carried out on two 
sample systems, one consisting of a single generator (2-bus system) and the other containing 
multiple generators (10-bus system) Finally, the stability region in the parameter space has 
been characterised considering the saddle- node bifurcation and Hopf bifurcation 

Since the phenomena of bifurcation and chaos are relatively new to the power system field, 
some of the basic definitions and concepts of these phenomena have been briefly explained in 
this Chapter and also in Chapter 5 


4.2 Local Bifurcation 

4.2.1 General 

Consider an autonomous dynamical system described by the following parameterised vector field 

y = g(y, A), yeR", \eRF (4 1) 

where g is a C T function on some open set m RJ 1 x R p The degree of differentiability will be 
deter min ed by the need for the number of terms to be retained in Taylor expansion of equation 
(4 1) Usually C 5 will be sufficient Let equation (4 1) has a fixed point at (y, A) = (yo, A 0 ), i e 

5 r (yo,A o ) = 0 (4 2) 

Two questions immediately arise’ 

1 Is the fixed point stable or unstable ? 

2 How is stability or instability affected if A is varied ? 

To check the stability of a fixed point, the first step is to examine the linear vector field obtained 
by linearising equation (4 1) about the fixed point (y, A) = (yo,Ao) The linear vector field is 
given by 

(4 3) 


t = D v g(y 0 ,\ 0 )£, f € R n 
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where £ = y — yo, D y g == j- is Jacobian matrix 

If the fixed point is hyperbolic (1 e none of the eigenvalues of D y g(yo, Ao) lie on the imaginary 
axis), the stability of ( y 0 , A 0 ) of equation (4 1) is determined by the linear equation (4 3) 

If the fixed point is hyperbolic, then varying A slightly does not alter the nature of the 
stability of the fixed point If D y g(yo , Ao) is hyperbolic and is invertible, then by implicit 
function theorem, there exists a unique C T function y(A), such that 

9(y( A),A) = 0 (4 4) 

for A sufficiently close to Ao with y(A 0 ) = yo 

Now by continuity of the eigenvalues with respect to parameters, for A sufficiently close to A 0 , 
D y g(y(X), A) has no eigenvalue on the imaginary axis Therefore, for A sufficiently close to A 0 , the 
hyperbolic fixed point (yo, Ao) of equation (4 1) persists and its stability type remains unchanged 
To summarise, m a neighbourhood of Ao an isolated fixed point of equation (4 1) persists and 
always has the same stability type However, if the fixed point (y 0 , A 0 ) of equation (4 1) is not 
hyperbolic i e when D y g(y 0 , A 0 ) has some eigenvalues on the imaginary axis, radically new 
dynamical behaviour can occur for A very close to Ao , for example fixed points can be created 
or destroyed and time-dependent behaviour such as periodic, quasiperiodic or even chaotic 
dynamics can be created In certain sense, the more are the eigenvalues on the imaginary axis, 
the more exotic the dynamics will be 

4.2.2 A Zero Eigenvalue 

If D y g(yo, Ao) has a single zero eigenvalue with the remaining eigenvalues having non zero real 
parts, the following three types of bifurcations may occur. 

4 2.2.1 Pitchfork Bifurcation 
Consider a vector field 

x = f(x , g) = gx — x 3 , x € R 1 , g € R 1 (4 5) 

The fixed points of equation (4 5) are x = 0 and i 2 = g as shown in Fig 4 1 At g = 0, the 
associated Jacobian matrix is singular 

For g < 0, there is one fixed point l e x = 0 which is stable For g > 0, x = 0 is still 
a fixed point, but two new fixed points have been created at g = 0 and are given by j 2 = g 
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Figure 4 2 Transcntical Bifurcation 

In the process, x = 0 has become unstable for p > 0, while other two fixed points are stable 
This type of bifurcation in which the original branch changes its stability character and a pair 
of fixed points are born when the Jacobian matrix is singular, is called pitchfork bifurcation 

4.2. 2. 2 Transcritical Bifurcation 
Consider a vector field 

x = f(x, p) = px — x 2 , xGfi 1 , p G R 1 (46) 

The fixed points of equation (4 6) are x = 0 and x = p as shown in Fig 4 2 

At p = 0, the associated Jacobian matrix is singular For p < 0, there are two fixed points 
x = 0 is stable and x = p is unstable The two fixed points coalesce at p = 0 and for p > 0, 
x = 0 is unstable and x — p is stable Thus, an exchange of stability has occurred at p = 0 
This type of bifurcation is called transcritical bifurcation 
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Figure 4 3 Saddle-node Bifurcation 

4. 2. 2.3 Saddle— node Bifurcation 

Consider a vector field 

* = /(*» P) = P ~ x e R 1 , fie R 1 (4 7) 

The fixed points of equation (4 7) are /x = x 2 This represents a parabola m the fi-x plane as 
shown in Fig 4 3 For fi < 0, equation (4 7) has no fixed point For fi > 0, there exists two 
fixed points, one of them being stable and the other unstable At fi = 0, the associated Jacobian 
matrix is singular This particular type of bifurcation where on one side of a parameter value 
there is no fixed point and on the other side there exists two fixed points is referred to as a 
saddle-node bifurcation 


4.3 A Pair of Pure Imaginary Eigenvalues 

When the matrix D y g(y 0 , Ao) has a pair of purely imaginary eigenvalues, a Hopf bifurcation (HB) 
arises and is characterised by emergence of a periodic orbit (limit cycle) around an equilibrium 
point At Hopf bifurcation, the number of fixed points are preserved The Hopf bifurcation is 
said to be subcntical when unstable periodic orbit emerges around stable fixed point as shown 
in Fig 4 4 The Hopf bifurcation is said to be supercritical when stable periodic orbit emerges 
around unstable fixed points as shown in Fig 4 5 The nature of Hopf bifurcation (subcriti- 
cal/supercritical) can be determined by solving the set of nonlinear differential equations 
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Figure 4 4 Subcntical Hopf Bifurcation Figure 4 5 Supercritical Hopf Bifurcation 

4.3.1 Determination of Static Bifurcation 

Occurrence of transcntical bifurcation in power systems have been shown m very few literatures 
[105, 165] Transcntical and pitch fork bifurcations are non generic whereas the saddle-node 
bifurcation (SNB) is generic 

The saddle-node bifurcation in a power system model occurs when the number of solutions 
to the load flow equations changes under parameter variation The singularity of D y g(yo , A 0 ) 
has been popularly used to identify SNB At SNB, a stable operating point and a saddle coalesce 
to form an unstable equilibrium at which the Jacobian is singular and has a zero eigenvalue A 
zero eigenvalue of D y g(y 0 ,\ 0 ), also arises when a type-2 and a type-1 equilibrium coalesce to 
form a fold bifurcation (FB) or limit point (LP) 

Several researchers such as Alvarado et al [53] and Camzares et al [103,126] have presented 
a method called Point of Collapse (POC) method to determine SNB In POC method, the 
parameter A is regarded as a variable and the condition of zero eigenvalue of the Jacobian is 
appended to the set of nonlinear equations Thus, if v is right eigenvector corresponding to 
the zero eigenvalue of the Jacobian, the SNB or LP can be obtained by solving the following 
equations considering A € R 1 , 


g(yA) = o 

(4 8) 

D y g(y, A)t> = 0 

(4 9) 

IMI 

(4 10) 


The POC method can also be formulated in terms of the left eigenvector One of the major 
limitations of the POC method is that it yields only (y, A) corresponding to a SNB or LP In 
the neighbourhood of the SNB or LP, the equation (4 8) can not be solved by Newton Raphson 
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method as the Jacobian becomes ill-conditioned To overcome this difficulty, the continuation 
method [50] which is a modification of Newton-Raphson method has been used by Ajjarapu 
et al [100] and Iba et al [86] In the present work AUT086 software [34] has been used 
for numerical analysis which also utilises the continuation method The continuation method 
involves a predictor step and a corrector step which have been described below 


Predictor Step: Let a solution t/o corresponding to Ao is known Prediction of the next 
solution can be made by taking an appropriately sized step in a direction tangent to the solution 
path Thus, the first task in the predictor process is to calculate the tangent vector This tangent 
calculation is derived by first taking the derivative of both sides of the equation (4 8) 

= 0 (4 11) 


<%(y,A)] = \9y,g\] 


dy 
d A 


On the left side of this equation is a matrix of partial derivatives multiplied by a vector of 
differential The above equations involve (n + 1) variables 

Thus one more equation is needed to have a unique solution for the tangent vector This 
problem can be overcome by choosing a non-zero magnitude (say unity) for one of the compo- 
nents of the tangent vector (Let d\ = 1) This results in 


9y 

9 \ 


dy 


' 0 ' 

0 

1 


. dx 


1 


The AUT086 software usually chooses the largest component of the tangent vector to have a 
non-zero magnitude The corresponding parameter is called continuation parameter 

Once the tangent vector has been found by solving equation (4 12) the predicted solution 
(y*> A*) can be obtained as follows 


V* 

X* 


y o 
Ao 


4 - a 


dy 

dX 


(4 13) 


The software has the capability to choose variable magnitude of step size a to improve the 
convergence 


Corrector Step: Once the prediction has been made, a method of correcting the approximate 
solution is needed This is achieved by specifying the value of the parameter A at its piedicted 
value The new set of equations would be, 



CHAPTER 4 OBSERVANCE OF LOCAL BIFURCATIONS 


62 


aivA) 

A- A* 


(4 14) 


The equation (4 14) is solved using Newton-Raphson method with predicted solution as initial 
guess The convergence of this equation yields a point say (yi,Aj) on the branch and repeating 
the above process the solution curve in (n + 1) dimensional space R n x R 1 is traced Along the 
solution branch, the eigenvalues of matrix y v (y 0 , A 0 ) is monitored to detect the local bifurcations 
such as saddle-node and Hopf bifurcations 


4.3.2 Determination of Hopf Bifurcation 

For determination of Hopf bifurcation, BIF0R2 program [11] has been used It makes use of 
center manifold theory [11,21,80] The structure of orbit near (y, A) = (yo, A 0 ) is determined 
by the vector field m equation (4 1) restricted to the centre manifold This restriction results m 
a p-parameter family of vector fields on a two-dimensional center manifold However BIFOR2 
program utilised in the present work, caters to one parameter only 

On the center manifold, the vector field in equation (4 1) has the following form [11,80] 


Z i 


ReX(n) —Im\(fi) 

' Zx ' 

4- 


Z 2 . 


Re\(fi) 

Z 2 

i 

*). 


(Z x ,Z 2 ,/x) G R 1 x R 1 x R} 

where, / 1 and / 2 are nonlinear functions in Z\ & Z 2 and A(/z) = a(n)+ju>(fi) are the eigenvalues 
of the vector field linearised about the fixed point If a(/x) = 0 and ^ 0, the fixed point is 
not hyperbolic The dynamical system undergoes a Hopf bifurcation [11,21,38,80] and a periodic 
orbit emerges In order to analyse the dynamics, the equation (4 15) needs to be simplified and 
cast into its Poincare-normal form [11,80] as follows, 

Z x — a(n)Zi — uj{h)Z 2 4- ( a(n)Zi — b{n)Z 2 )(Z 2 4- Z 2 ) 4- 0(\Z x f , \Z 2 \ 5 ) (4 16) 

and 

Z 2 = 4- a(n)Z 2 4- (a(/i)Zj 4- b(n)Z 2 )(Zi 2 + Z 2 2 ) 4- 0{\Z X \ , \Z 2 \ ) (4 17) 

For the sake of simplicity, it is desirable to work in polar coordinates ( vL6 — Z\ 4- jZ 2 ) and the 
equations (4 16) and (4 17) will assume the following form, 
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r = a(fi) + a(/z)r 3 + 0(r 5 ) 

(4 18) 

9 = uj(fi) + b(/j,)r 2 + 0(r 4 ) 

(4 19) 

Upon Taylor series expansion of coefficients around pt = 0, the equations (4 18) and (4 19) yield, 

r = a'(0 )/J,r + a(0)r 3 + 0(^i 2 r, fir 3 , r 5 ) 

(4 20) 

9 = u>(0) + u/(0)/i + 6(0)r 2 + 0(fi 2 , fir 2 , r 4 ) 

(4 21) 

Neglecting higher order terms, the dynamics of equation (4 1 ) m 
given as follows 

Poincare-normal form can be 

r = a'(0 )fir + a(0)r 3 

(4 22) 

9 = w(0) + u'(0)fi -f- 6(0)r 2 

(4 23) 


The value of r > 0 and fi for which r = 0 and 9 — 0 correspond to periodic orbits of equations 
(4 22) and (4 22) Therefore, for — oo < fia'(0)/a(0 ) < 0 and fi sufficiently small, the periodic 
orbit for equations (4 22) and (4.23) is given as follows 

(r(*M(<)) = 

where 9 0 is the initial phase angle The periodic orbit is asymptotically stable for <z(0) < 0 
(supercritical) and unstable for a(0) > 0 (subcntical) 

4.3.3 Determination of Unstable Manifold 

Consider a type-1 unstable equilibrium point (u e p) say x Let be the only eigenvalue with 
a positive real part and its eigenvector be r The one dimensional unstable manifold W u (x) of 
x can be broken by removing x from the manifold, into two half manifolds W+(x) and VU“(i) 
Each part in itself is an invariant manifold consisting of a single trajectory as shown in Fig 4 6 
Calculation of each part involves following two major steps [45,59,128] 


-/i<y(0) 

a(0) 


w(0) + M0)- 


6(0)a'(0) 


t + 9 o 


(4 24) 
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Figure 4 6 Half-manifolds at an equilibrium point 

Locating a point on W+(x): Near the controlling (u e p ) W+(x) can be approximated to 
the first order, by the unstable eigenvector 7 / ;i Therefore, m order to derive IT “ (x), a point 
close to x that lies on can be chosen More specifically, for some small a > 0. choose the 
point 

x a =x + arjn (4 25) 

Care should be taken when choosing a If a is too large, x a might not lie close enough to W+ ( x ) 
in which case the trajectory through x a is not a good approximation to W+(x) a is chosen to 
locate a point x a on W+(x), as follows 

(1) Given a small enough e > 0, create a ball around x such that | x — x |< e holds 

(2) Give the initial a (a > e) 

(3) Obtain x a = x + arj^ 

(4) Integrate the system from x a in reverse time by a specified number of steps 

(5) If the system trajectory remains in the ball, stop Otherwise, let cx — \a. and go to 
Step 3 

Integrating on W+(x): Removal of x from W u (x), splits W u (x) into two segments The 
larger segment of VF£(x) is found by integration from x a in forward time The smaller segment 
is the line segment connecting x and x a 

Similarly can be computed by considering the reverse direction of the eigenvector In 

order for an (u e p ) to be on stability boundary, one of the two half manifolds should conveige 
to the stable equilibrium point (s e p ) 
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4.4 Power System Model for Bifurcation Studies 

4.4.1 Generator Model 

Different models of the generator [7,25,176] used in the studies are given as follows, 

4.4 1.1 Second Order Generator Model 

In this model, the generator is represented by a constant emf source behind the transient reac- 
tance and modelled by the swing equation 

dS 


dt 


= U 


du 


M g — = Pm — P ge — Du 


(4 26) 
(4 27) 


where <5 is generator’s internal bus voltage angle, u the speed deviation, M g the moment of 
inertia of generator, P ge the electrical power output of the generator, D the damping constant 
and Pm is the mechanical power input to the generator 

4 4 1.2 Third Order Generator Model 

This model ignores the saliency and accounts for the changes in the flux linkage of the field 
winding although the changes in the flux linkages of other windings are ignored Therefore, 
there is one more state equation for the third order generator model, in addition to the swing 
equations (4 26) and (4 27) 


T' d0 




X d -X' d dt 


Ejd - E' Q 


X d -X' d 




(4 28) 


where E' q is the voltage behind the d-axis transient reactance X d , T' do is the d-axis transient 
open circuit field time constant, X d is the d-axis steady state reactance, Ej d the voltage applied 
to the field and Q ge is the reactive power output at the internal node of the generator 


4.4 1 3 Fourth Order Generator Model 

If the effect of saliency is considered, the changes in flux linkages of the field windings have to 
be accounted for along the d- and q-axes Therefore, two more additional state equations (4 29) 
and (4 30) along with the swing equations (4 26) and (4 27) have to be considered 
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dE' 

= Ef d + - X' d )I d - E' q 

T'J-§ = ~(X, - X',)I, - E' d 


(4 29) 
(4 30) 


along with the voltage equation, 

' v q -v q 

. E' d -Vd 

where, T' qo is the q-axis open circuit transient time constant, R a the armature resistance, X g the 
steady state q-axis reactance, X’ q the q-axis transient reactance, E' d the voltage behind X q , I d 
& Iq the d- and q- axes armature current and V d & V q are the generator terminal voltage along 
d- and q-axes 


, 

p 


i — 

-r 

i 

i 

R a 

1 

p? 


(4 31) 


4.4.2 Load Model 

The load characteristics are known to have profound effect on the system dynamics [184, 185] 
The classical models of load such as constant P and Q , constant current or constant impedance 
have been reported to be insufficient and inappropriate to capture voltage dynamics [40, 114] 
Different load models used for the bifurcation studies are as follows 

4.4 2.1 Load Model— I 

This load model considers a composite load consisting of a dynamic induction motor (Pjm + 
jQim ) in parallel with a static constant P, Q load (P, + jQ s ) The model of induction motor 
has been taken from ref [40] and expressed as a function of rate of change of its terminal voltage 
magnitude and angle The model is given as follows 


< 

ii 

(4 32) 

Q, = Qi 

(4 33) 

P,m = P«+K^9+K„(V + TV) 

(4 34) 

Qim = $0 + Kqu6 + K qv V + K qV2 V 2 

(4 35) 


where P 0 + jQo is the constant part of induction motor load, K puJ , K pv , T, K quJ , K qv and K qV2 
are load coefficients 
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Figure 4 7 Model of induction motor 


4. 4. 2.2 Load Model-II 

This model also represents a composite load consisting of static and dynamic components The 
equations describing the load (P d + jQd) using this model [104, 150] are as follows 


Pd = 8 pf + DO + aV 

(4 36) 

Q d = i^/l-pf 2 + b6 + kV 

(4 37) 


where i is nominal MVA demand at the load bus, pf the power factor and D, a , b and k are 
the dynamic load parameters 

The model-I and model-II tacitly ignore the stator/ network transients and if stator tran- 
sients are included, these models are not be valid [150] 


4. 4. 2.3 Load Model-Ill 


This model represents a third order induction motor model that ignores stator transients If 
V m L8 m be the voltage behind the stator resistance R\ and transient reactance X' of the motor 
as shown in Fig 4 7, ur its rotor speed, the equation governing its dynamics can be written 
as [99] 


r 0 dv m 

X 3 -X' dt 


Qm Vjji 

.k:~xt^Y'. 


ds m , P m {X, - X') 

—— = U)R — UJ S + TTJTr' 


(4 38) 
(4 39) 




duR 

dt 


= Pm —Tl OJr 


(4 40) 


where M m = Moment of inertia of induction motor 
T' a = Open circuit transient time constant 
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Figure 4 8 IEEE type-1 exciter 
X s k X' = Steady state and transient reactance 

P m k Q m = Real and reactive power drawn at the internal node of the motor 
u>r = Rotor speed 

Tl = Ao + + CoUJft, speed dependent load torque 

w, = co 0 + where u> 0 is the synchronous speed and <f> 3 is the phase angle 

of the terminal voltage 

4.4.3 Exciter Model 


IEEE type-1 exciter [7,25] as shown m Fig 4 8 has been considered in this study and equations 
governing the dynamics is given as follows 


T e ~ = v a -(K e + S E )E fd 

(4 41) 

T 'ii = K,B > d ~ v ' 

(4 42) 

T a ^ = K a (V ref -V t -V f )-V a 

(4 41) 

and V a m,n <V a < V™ ax 


where S E = A ex exp (B ex Ef d ) is a saturation function 


V a = Regulator amplifier output voltage 
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P max 



Figure 4 9 General purpose speed governing system 

K a = Regulator amplifier output gam 

T a = Regulator amplifier output time constant 

Vf = Stabilising circuit output voltage 

Kj = Stabilising circuit output gam 

T f = Stabilising circuit output time constant 

K e = Exciter gam 

T e = Exciter time constant 

V re f = Reference voltage setting 

V t = Terminal voltage 


4.4.4 Speed Governor Model 


A general purpose governor model [7] as shown m Fig 4 9 has been used and its model can be 
described by the following equations, 


Tl dt 


u>0 -0J - T 2 j t f 


LjJqR, 


X l 


eg 


(4 44) 


Xl — X\ + Pmo 


(4 45) 


(4 46) 


X$ = X 2 
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(4 47) 
(4 48) 

(4 49) 

(4 50) 

where R eg = Turbine steady state regulator setting or droop 
Pmax — Maximum turbine output 
Ti = Governor response time constant 
T 2 = Hydro reset or pilot valve time constant 
Tz = Servo or hydro gate time constant 

T 4 = Steam valve bowl time constant (zero for hydrogovernor) 

T 5 = Steam reheat time constant or half hydro water starting time constant 
F = pu shaft output ahead of reheater (-2 0 for hydro units) 


Pmax ^ ^3 ^ 0 

rp dX 4 _ v v 

l3 ~dT - x *- x < 
dX 5 


and 


rp dPfi4 y | prn dX$ 

T5 ~dF ~ x * + FT5 ~dT ~ 


4.4.5 Hard Limiter Model 


The hard limiters used in the excitation control as well as in turbine speed governing control 
loop for providing ceiling to amplifier output voltage and generator real power output The 
presence of hard limiters make the nonlinear equations non-analytic Thus the model can not be 
straightaway analysed for Hopf bifurcation using BIFOR2 [8] which is based on the assumption 
that the function is C T (r > 5) l e analytic m nature In order to overcome this difficulty, a 
new model of the limiters is proposed representing them by tanhyperbohc function 

A typical hard limiter is shown in Fig 4 10 having input variable as X and output as Y It 
is seen that if X is between —Xlimit and Xlimit , V = X i e the slope of the curve is unity 
and when X > X limit or X < -X limit, V = ±X limit i e the slope is zero In order to 
approximate the hard limiter, a function of the following form has been selected 


Y = X limit tan h 


(-*-) 

\X LIMIT' 


(4 51) 


This function has unity slope at X = 0 and saturates when X > Xlimit or X < — Xlimit Io 
order to establish the validity of this function, a sample study was conducted on the excitation 
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Figure 4 10’ A typical hard limiter 

system for Xlimit = 3 5 pu and the signal YvsX has been plotted m Fig 4 11 and 4 12 for 
the sake of comparison It can be observed that the new function as given by equation (4 51) 
closely approximate the hard limiter Hence, tanh function has been used to the model hard 
limiter action which is smooth and analytic 

4.5 Case Studies 

4.5.1 Static Bifurcations 

Static bifurcations arise when number of fixed points changes if a system parameter is varied 
In order to illustrate the occurrence of static bifurcations and its implications, studies were 
conducted on a 3-bus and a 9-bus sample systems derived from ref [66] and [7], respectively 
and also presented in Appendices-B and E respectively The studies were assisted by AUT086 
software [34] on HP-9000/735 computer systems and are presented below 

4. 5. 1.1 3-Bus System 

The topology and data for 3-bus system is shown in Fig B 1 in Appendix-B It consists of an 
infinite bus, a generator and a load bus The generator is modelled by using its second order 
model equations (4 26) and (4 27) and load by load model-I given by equations (4 32) to (4 35) 
With this formulation, the resulting equations become purely a system of ordinary differential 



CHAPTER 4 OBSERVANCE OF LOCAL BIFURCATIONS 


72 



x 


Figure 4 11 Hard limiter output 
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Figure 4 12 Tanh limiter output 
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Voltage 



Bifurcatiopn Parameter Q static 


Figure 4 13 Variation of voltage magnitude V vs bifurcation parameter Q\ for 3-bus system 

equations The static reactive load component Q x was taken as bifurcation parameter The 
bifurcation diagram has been plotted showing variation of load bus voltage V versus the bifur- 
cation parameter Q\ is given in Fig 4 13 The following bifurcations are observed from Fig 
4 13 

• HB at Q x = 6 272709 and 6 542023 

• SNB at Qj = 6 543211 

• LP at Qr = 0 7841441 and 1 978656 

From Figs 4 13, it is seen that there are different numbers and types of fixed points in different 
regions of the parameter They are listed as below. 

Region A (0 < Q\ < 0 7841441) one stable and one type-1 unstable fixed points 
Region B (0 7841441 < Q x < 1 978658) one stable, two type-1 unstable and a type-2 
unstable fixed points 

Region C (1 978658 <Q\< 6.272709) one stable, one type-1 unstable fixed points 
Region D (6 272709 < Q\ < 6 542023) one type-2 unstable and one type-1 unstable 
fixed points 

Region E (6 542023 < Q x < 6 543211) one stable and one type-1 unstable fixed points 
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All along the bifurcation parameter, it is seen from the bifurcation diagrams of Fig 4 13 that 
unstable equilibrium points exist along with stable equilibrium point In view of this, the stable 
equilibrium point is not globally asymptotically stable, and its boundary of region of attraction 
is given by the union of the stable manifolds of type-1 equilibrium points 

Concentrating in the region B, it is found that there exist two type-1 unstable equilibrium 
points One type-1 u e p is characterised by high voltage and the other one by low voltage It 
will be of interest to know about the unstable manifold because for a type-1 u e p , the critical 
unstable trajectory can be approximated, after a finite time, by the unstable manifold of the 
type-1 controlling u e p [128] Since the dimension of unstable manifold of a type-1 u e p is one, 
the unstable mode relative to critical unstable trajectory whose corresponding controlling u e p 
is type-1, is uniquely defined by the one-dimensional unstable manifold For a type-2 (or higher 
type) controlling uep, the critical unstable trajectory can be approximated, after a finite time, 
by the unstable manifold of the controlling uep But the dimension of the unstable manifold of 
a type-k uep is k This causes the unstable mode relative to the critical unstable trajectory, 
whose corresponding controlling uep is type-2 (or higher), not to be uniquely defined by the 
corresponding unstable manifold [128] 

Figs 4 14 and 4 15 show the half unstable manifold, W^x) computed m the reverse direction 
of the eigenvector associated with the real positive eigenvalue for the high voltage uep at 
Q\ = 1 0 pu From Figs 4 14 and 4 15 it is evident that the half manifold for this uep 
converged to the stable equilibrium point (sep) Thus, this type-1 uep will be on the 
boundary of the region of s e p 

Fig 4 16 and 4 17 show the plot of the other half manifold, W+{x) m the forward direction 
of the eigenvector It is seen that the state variables 8 — (x[l]) and u = (x[2]) in Fig 4 16 
associated with the generator becomes unbounded while the variables 6 = (x[3]) and V = (x[4]) 
m Fig 4 17 remain bounded Therefore, a post-fault trajectory for which this type-1 uep is 
the controlling one, will be accompanied by angle instability (or synchronous instability) only 
The voltage instability would not be observed m this case 

Figs 4 18 and 4 19 show the plot of half manifold in the reverse direction of the eigenvector 
associated with the positive real eigenvalue of low voltage type-1 uep corresponding to the 
same parameter value Q x = 1 0 pu , as studied earlier for the case of high voltage uep It is 
established from the Figs 4 18 and 4 19 that the uep will be on the stability boundary of 
sep as it converged to the sep Figs 4 20 and 4 21 show the plot of the other half manifold 
of this uep determined m the direction of the eigenvector It is seen from Fig 4 20 that the 
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Figure 4 14- W*{x) plot of 8 and u for high voltage uep of 3-bus system 
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Figure 4 15 W^(x) plot of 9 and V for high voltage u e p of 3-bus system 
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S = (x[l|) and u) = (z[2]) variables remain bounded while m Fig 4.21 the variables 9 — (x[3]) 
and V = (x[4]) become unbounded If this type-1 u e p is the controlling uep, the fault-on 
trajectory will be accompanied by collapse of the voltage 

However for parameter value Q\ = 1 0 pu. both the type-1 uep would exist Therefore, 
the system will be vulnerable to angle as well as voltage stability in the region B 

In region A, the system will be prone to angle instability For Q i > 1.978658, the system 
will be threatened by only voltage stability However, for Qi > 1 978658, several types of global 
bifurcations would exist which have been explored in the next Chapter (Chapter 5) 

4. 5. 1.2 9— Bus System 

The single line diagram and data of 9-bus system is presented m Fig. E and Appendix-E 
respectively The network was slightly modified to take the line 5-4 out The existing static 
load of (125 + j50) MVA was reduced to (100 + j30) MVA This change m loading was done 
since the load flow was not converging after the line 5-4 was taken out 

The system consists of 3 generators The generator at bus-1 was taken as the reference bus 
and tbpse at bus-2 and 3 were modelled by swing equations (4 26) and (4 27) All the base 
case loads have been modelled as constant impedance type. These constant impedance loads 
were clubbed with the Yf, ua of the system An additional load represented by the load model-II 
(equations (4 36) to (4 37) has been considered at bus-5 

This formulation results in a nonlinear system of ordinary differential equations consisting 
of six state variables viz the internal machine voltage angle S and speed deviation u) for the two 
generators (generator # 2 & 3) and voltage magnitude V and angle 6 at bus-5 Using AUT086 
software [34], the dynamics was studied on HP-9000/735 computer system, by varying l The 
power factor was assumed to be unity and the other parameters of the dynamic load D , a , b 
and k were taken to be 0 05 (sec), 0 0, 0 0 and 0 1 (sec), respectively 

The bifurcation diagram for variation in parameter l has been drawn as shown m Fig 4 22 
(load bus voltage magnitude V vs l) Some of the important points on the bifurcation diagram 
are as given below* 

• SNB at t = 1 195323 

• HB at t = 1 116463 and 1 192705 


• LP at i = 0 2351554 and 0 6105376 
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Figure 4 19 PF“(x) plot of 0 and V for low voltage u e p of 3-bus system 
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xl2] 

Figure 4 20 W+(x) plot of 6 and u for low voltage u e p of 3-bus system 



xlll 

Figure 4 21 W^(x) plot of 6 and V for low voltage u e p of 3-bus system 
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Voltage 



Figure 4 22 Variation of voltage magnitude V vs bifurcation parameter l for 9-bus system 

It may be seen from the bifurcation diagram m Fig. 4 22 that for l < 0 2351554 and l > 
0 6105376 there is only one type-1 u e p These two u e.ps are markedly different One is 
characterised by high voltage and the other by low voltage The low voltage u e p coalesces 
with s e p to form SNB 

The one-dimensional manifold of both the type-1 u e p corresponding to parameter ^ = 04 
has been studied Figs 4.23 to 4 25 show the plot of one dimensional half unstable manifold 
for high voltage u e p. computed in the reverse direction of the eigenvector corresponding to 
real positive eigenvalue It is seen from Figs 4 23 to 4 25 that this half manifold converged 
to the s e p indicating that this u e.p is on the boundary of the region of attraction of the 
s e p Figs 4 26 to 4 28 show the plot of the other half manifold for this u e p computed in the 
forward direction of the eigenvector It is clear from these figures that not only the synchronous 
machines loose angle instability but also the system voltages collapse 

The one dimensional half unstable manifold of the low voltage small angle u e p computed 
in the reverse direction of the value converged to the s e p (Figs 4 29 to 4 31), thus, confirming 
that this u e p is also on the boundary of the region of attraction of s e p The other half 
unstable manifold computed m the forward direction of the eigenvector as shown in Figs 4 32 
to 4 34, indicate that load voltage collapses and the synchronism amongst generators is also lost 
Therefore, as noted above for this system, the voltage and angle instability will occur together 
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4.5.2 Dynamic Bifurcation 

Dynamic bifurcations arise when the stability characteristic of the equilibrium points changes 
Hopf bifurcation is such an example of the (local) dynamic bifurcation This section presents the 
study of Hopf bifurcation on two sample systems, one consisting of two buses (a generator and a 
induction motor load bus) as shown m Fig A.l and the other a 10-bus system as shown m Fig 
F 1 The 10-bus system is derived from the 9-bus system of Appendix-E An induction motor 
was connected to bus-6 The studies were assisted by BIF0R2 program [11] on HP-9000/735 
computer system for the following cases 

Case— I: For different combinations of induction motor load torque coefficients (Ao,i?o, C 0 ), 
BIF0R2 program was used to find out critical value of generator AVR amplifier gam K a The 
bifurcation parameter A for this study was taken as K X J 2 which ensures positive value of gam 


Case— II: For the base value of AVR amplifier gam and different combination of load torque 
coefficients (Ao, B 0 ), BIF0R2 program was used to find out critical value of torque coefficient 
Co The bifurcation parameter A for this study was taken as C\j 2 . 

The results of the above two Cases for the two sample test systems are given below 

4. 5. 2.1 Results of 2-Bus System 

For both Case-I and Case-II studies on 2-bus system, the fourth order model of generator 
(equations (4 26) (4 27) (4 29) and (4 30)) together with AVR (equations (4 41) to (4 43)), speed 
governing control system (equations (4 44) to (4 50)) and load model-III (equations (4 38) to 
(4 40)) were used The hard limiters were represented as described in section 4 4 5 The reference 
voltage to AVR was selected as 1 1 pu and reference setting of governor was considered as 1 048 
pu The various data used in the system studies are given in Appendix-A 

For Case-I study, different values of all the load coefficients (Ao, Bo, Co) were fixed between 
0 1 to 0 8 The corresponding values of critical gain K a obtained from BIFOR2 and some of 
the other parameters are given in Table 4 1 For Case-II studies, the different values of load 
coefficients Ao and Bo were selected in the range of 0 00 to 0 81, the K a was fixed at 10 00 and 
critical values of Co and some of the other parameters obtained from BIFOR2 are listed in Table 
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Figure 4 29 W_(x) plot of 62 — x[l] and u >2 = x[2] for low voltage u e p of 9-bus system 
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Figure 4 30 W^(x) plot of S 3 = x[3] and oj 3 = x[4] for low voltage u e p of 9-bus system 
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Table 4 1: Case-I results of 2-bus system 


Load Coefficients 
Aq — Bo = Co 

Critical 
gain K a 

«(0) 

6(0) 

| 1 

M 

cu(0) 

xlO -2 

0 1 

15.9949 

-161 972 

-76 562 

0 2483 

0 1681 


02 

15.6009 

-040 397 

-18 514 




03 

14.8591 

-016 647 

-07 218 

0 2756 

0.1717 


04 

13.6547 

-007 857 

-03 130 

0 2999 

0 1747 

IB 

05 

11.9119 

-003 641 

BtlKTiEl 

0 3314 



06 

09.6294 

-001 391 

MilMrJ 

0 3689 

0 1825 

0 5969 

0 7 

06.8670 

-000 182 


0 4018 


lilSEEl 

08 

03.4227 

-000 416 



■iilflffft 



Table 4 2: Case-I I results of 2-bus system 


Load Coefficients 
Ao = Bo 

Critical 
value of Co 

“(0) 

5(0) 

eksm 

«'(0) 

xlO -2 

a,(0) 

xl0“ 2 

00 

HUBS 

-1 7651 

-0 6821 

0 8892 

0 1043 

0 6044 

0 1 


-1 7384 

-0 6580 

0 8671 


0 6047 

02 


Bum 

-0 6351 

Kraus! 

■wnwii 

0 6049 

05 

BiBOi 




mim:i 

0.6051 

07 



-0 5392 

0 4463 


KW 

08 

■eessehi 

-1 6242 

-0 5237 




0 81 

bohh 

-1 6237 

-0 5222 


0 0086 



4 2. It can be observed that the critical values of K a m Case-I and that of Co m Case-II studies 
which result into Hopf bifurcation lie m their normal operating range Further for periodic orbit 
to exist, /xa , (0)/a(0) should be strictly less than zero This implies that the periodic orbit will 
exist if /x = A — A 0 > 0, i e for Case-I study, if K a > K acrttlcal and for Case-II if C 0 > C 0criUcat 
It is seen that either when loading increases or exciter amplifier gam increases, a complex 
pair of eigenvalues moves towards imaginary axis At critical value of parameters the stationary 
point looses stability via Hopf bifurcation An eigenvalue sensitivity analysis using the techniques 
given m [72, 176], was performed to ascertain as to which of the states caused this instability 
Participation factors of state variables m the critical mode were computed Their magnitude 
for a typical case study (corresponding to Aq = Bq = 0.1 and Co = 1 5130) are given in Table 
4 3 The study reveals that the unstable modes correspond to E' q and E/d states 
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Table 4 3 Magnitude of participation factor 2-bus system 


'0 00006 (£j 


0 47715(E/«t)_ 
0 00040(P m ) 



Fl gure 4 35 2-bus system plot of l.rmt cycle projected on E‘, and E,a plane 

For val.dat.on of BIF0R2 results, nuxnencal s.mulat.on was done using Runge-Kutta-Merson 
method for solv.ng the system of ODE and a typical plot of hunt cycle between state E, 

Ejd is given m Fig 4 35 


4 5 2.2 Results of 10— Bus System 

The generators were modeUed by thud order model (er.uat.on (4 26) to (4 28)) along with AV* 
speed governing control loop and hard hunters. A dynanuc load mode, at bus was cons. 

ror nopi uuul respectively The reference 

1 02774 pu 1 08087 pu and 1 090037 pn for generators # 1, 2 and 3, respec y 
1 UZ ‘ ' 4 pu ’ 1 P , bmU1 k 1 03 and 1 05 pu for generator # 

power settings of speed governing system were taken 
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Table 4 4 Case-I results of 10-bus system 


Generator 

Critical 

«(0) 

6(0) 

«'(0) 

u/(0) 

“(0) 

Exciter No 

gam K a 

X 

1 — 1 
o 

1 

t— » 

xlO -1 

1 

o 

rH 

X 

1 

O 

r— I 

X 

xlO' 2 * 

1 

11 5588 

-0 5999 

-0 3424 

0 2203 

0 4068 

0 4466 

2 

10 1675 

-0 5480 

-0 3461 

0 0229 

102 30 

0 4492 

3 

10 4850 

-0 5392 

-0.3177 

0 7825 

03 372 

0 4489 


Table 4 5 Case-II results of 10-bus system 


Load Coefficients 
Ao = Bo 

Critical 
value of Co 

«( 0) 
xlO -1 

6(0) 

xlO -1 

«'(0) 

xlO -4 

u/(0) 

xl0~ 4 

"(0) 

xl(T 2 

0 01 

0 1836 

-0 5455 

-0 3121 

0 2980 

0 4537 

0 4474 

0 02 

0 1634 

-0 5482 

-0 3161 

0 3165 

0 4602 

0 4474 

0 03 

0 1431 

-0 5454 

-0 3061 

0 1626 

01287 

0 4474 

0 04 

0 1229 

-0 5528 

-0 3247 

0 3110 

0 4398 

0 4474 

0 05 

0 1027 

-0.5503 

-0 3265 

0 2790 

0 4436 

0 4474 

0 06 

0 0825 

-0 5501 

-0 3209 

0 2400 

0 3643 

0 4474 

0 07 

0 0623 

-0 5507 

-0 3232 

0 1271 

0 1130 

0 4474 

0 08 

0 0421 

-0 5517 

-0 3240 

0 1461 

0 1935 

0 4474 

0 09 

0 0218 

-0 5482 

-0 3163 

0 1153 

0 1753 

0 4474 

0 10 

0 0016 

-0.5510 

-0 3296 

0 0515 

0 1660 

0 4474 


1, 2 and 3, respectively Since static loads were represented as constant impedances, the Ward 
reduction method was used to reduce the system to a 4-bus system containing only generators’ 
and motor’s internal nodes. This reduction eliminates the algebraic variables and resulting 
are a set of nonlinear ordinary differential equations only One set of results was obtained for 

Case-I pertaining to the load coefficients (Ao = Bo = Co = 0 1) and the result of BIFOR2 
are presented m Table 4 4 For Case-II, seven sets of results were obtained for load coefficients 

(Ao, Bo) varying between 0 01 to 0 10 and results of BIFOR2 are presented in Table 4 5 

From these results it can be observed that the periodic orbits will exist for K a > K acr , tica i as 
computed by BIFOR2 for Case-I study and C 0 > Co crtttcal for Case-II study As the coefficients 
a(0) < 0 m all sets of studies, the bifurcation will be supercritical Participation analysis 
reveals that the dominant states that caused this instability correspond to E' q and Ejd The 
magnitude of participation factors of the critical mode for a typical case study (corresponding 
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Table 4 6 Magnitude of participation factor 10-bus system 


6'6l035(£J,) 

0 351 72^) 

0 11256(£;,) 

0 00037(6 2 ) 

0 00007 ( 63 ) 

0 00389(o> 2 ) 

0 00000(a;3) 

0 00317(01*) 

0 00994(5/(1,) 

0 36805( £,<*,) 

0 11315(£ /d ,) 

0 00010(4^) 

0 00071(V /2 ) 

0 00022 ( V/ 3 ) 

0 00029(14, ) 

0 09934(V a2 ) 

0 03052(I/ O 3 ) 

0 00008(A„) 

0 00008(X 51 ) 

0 00004(P Ml ) 

0 00008(A:i 2 ) 

0 00005(A 13 ) 

0 00020(X, 2 ) 

0.00012(^43) 

0 00067(^51 ) 

0 00039 (^ 53 ) 

0 00072( Pm 2 ) 

0 00045(/W 3 ) 

0 00002(6 m ) 

0.00060(F m ) 

0 00047(w r ) 







to A 0 = B 0 = 0 01 and C Q = 0 18211) are given m Table 4.6 

For validation of results, numerical simulation was performed using Runge-Kutta-Merson 
method for solving the system of ODE and a typical plot of limit cycle between E' q2 and Ej^ 
&: E' n and E/d 3 are shown m Figs 4 36 and 4 37, respectively 


4.6 Effect of Local Bifurcation on the Region of Stabil- 
ity in Parameter Space 

Further studies were done to observe the effect of static and dynamic bifurcations on the region 
of stability in parameter space on the 3-bus and 9-bus systems studied m section 4 5 1 Both real 
Pi and reactive loads Qi at the load node 2 and 5 for the 3 bus and 9 bus systems respectively, 
have been considered as bifurcation parameters and locus of saddle-node bifurcation and Hopf 
bifurcation has been drawn on P — Q plane using AUT086 [34] Fig 4 38 shows the locus of 
saddle-node and Hopf bifurcations for 3 bus system in P — Q plane and Fig. 4.39 for the 9 bus 
system It is seen that the effect of dynamic bifurcations such as Hopf bifurcation is to further 
reduce the region of stable operating parameters It is seen from Fig 4 38 that Hopf bifurcation 
occurs for 3-bus system when power factor of load is low while for 9-bus system (Fig. 4 39) it 
occurs when power factor is high 


4.7 Conclusion 

From the studies presented m this chapter the following conclusions are drawn 

(l) A power system may exhibit static and dynamic bifurcations even for variations of oper- 
ating parameters in their normal range 
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Real Power 



Figure 4 38 Loci of saddle-node and Hopf bifurcation for 3-bus system 


Reactive Power 



Figure 4 39 Loci of saddle-node and Hopf bifurcation for 9-bus system 
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( 11 ) For the same value of bifurcation parameter, the system may become unstable either 
purely due to angle instability or purely due to voltage instability or both may occur 
simultaneously depending upon disturbance 

(in) When a system undergoes Hopf bifurcation, the system looses linear stability. However, 
for small perturbations, the system becomes oscillatory It is shown that for certain choices 
of motor load torque coefficients Ao, Bo and Co or AVR gain K a values, the power system 
may undergo oscillatory type instability that is essentially bifurcation of a stationary point 
into a periodic orbit. 

(iv) The eigenvalue sensitivity reveals that the states responsible to cause Hopf bifurcation are 
E' q and Ejd when AVR, speed governor and motor dynamics are considered m the model 
along with the generator dynamics 

(v) The effect of Hopf bifurcation is to further reduce the stable operating range of parameter 

(vi) Hopf bifurcation may occur both at low or high power factor loadings 



Chapter 5 


Observance of Global Bifurcations 


5.1 Introduction 

Global bifurcations emerge out of periodic orbit and characterises a qualitative change m phase 
portrait not restricted to small neighbourhood of an equilibrium point When viewed m a con- 
tinuous time-setting, the bifurcation of limit cycles such as period doubling bifurcation (PDB), 
cyclic fold bifurcation (CFB) and bifurcation to torus (TRB) are the global bifurcation as they 
involve considerations over a large domain in phase space [161] Other global bifurcations may 
also arise in the study of system dynamics A periodic solution may collide with a saddle 
equilibrium m a saddle connection bifurcation A strange attractor may collide with a saddle 
equilibrium or periodic orbit and disappear in a boundary crisis Examples of some of these 
global bifurcations in electric power system models are given in ref [118,120, 154, 191] In math- 
ematical literatures [21,38,80], the global bifurcations refer to the homoclinic and hetrochnic 
bifurcations which are still an active area of research Kopell et al [17] and Salam [22,23] have 
studied the homochmc bifurcations in a model of power systems using Melnikov’s method 
PDB route to chaos and boundary crisis has been studied in ref [120, 122, 127, 154, 191] 
considering load and generator dynamics However, phenomena of voltage collapse is also greatly 
influenced by the on load tap changer (OLTC) and static var system (SVS) dynamics 

In this Chapter, extensive studies have been done to observe different types of dynamic 
bifurcations such as PDB, CFB, TRB, chaos and boundary crisis (blue sky bifurcation) on 
three sample systems namely 3-bus, 9-bus and 19-bus systems The effect of OLTC and SVS 
dynamics have been considered along with load and generator dynamics Existence of new 
dynamic bifurcations in power system models such as stable and unstable torus h chaos via 
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Figure 5 1 Period doubling bifurcation 

torus breakdown and voltage collapse due to some of these phenomenon has been explored 
To study the stability of the dynamic bifurcations and chaos behaviour, Floquet multipliers, 
Lyapunov exponents and dimension of the attractor have been computed 


5.2 Bifurcation of Periodic Orbit 

5.2.1 General 

The bifurcation of a periodic orbit (in continuous time setting) involves consideration over a 
large domain in phase space Under variation of the bifurcation parameter, a periodic orbit may 
undergo any of the following types of bifurcations 

5.2. 1.1 Period Doubling Bifurcation (PDB) 

In some of the systems, the stable limit cycle out of supercritical Hopf bifurcation may soon 
become unstable with another stable limit cycle of twice the period generated Emergence of a 
stable periodic orbit of twice the time period around another periodic orbit is known as period 
doubling bifurcation as shown in Fig 5 1 This acts as a gateway to chaos The original periodic 
orbit becomes unstable after bifurcation The PDB may also emerge out of unstable periodic 
orbit resulting into an unstable periodic orbit of twice the time period The original unstable 
orbit may or may not remain unstable Local analysis is generally not enough to predict the 
PDB 
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Figure 5 2 Cyclic fold bifurcation 
5 2 1.2 Cyclic Fold Bifurcation (CFB) 

This type of bifurcation arises when two periodic orbits coalesce Thus on one side of parameter 
at bifurcation point no periodic orbit exists and on the other side there exist two periodic 
orbits This is called fold bifurcation or limit point When a stable and an unstable periodic 
orbit coalesce, the resulting bifurcation is called cyclic fold bifurcation or saddle-node bifurcation 
of periodic orbits as shown in Fig 5 2 

5.2. 1.3 Torus Bifurcation (TRB) 

This type of bifurcation occurs at a point when a periodic orbit becomes quasiperiodic This 
generates frequencies of oscillations incommensurate or irrational to each other 

5.2.2 Stability of Periodic Orbit 

Consider an n th order autonomous parameterised vector field defined by the state equation 

x = f(x, A) (5 1) 

where / is a C T (r > 1) function on some open set in x R p (x,A) = (xo, A 0 ) at which 
/(x,A) = 0 is a fixed point of the vector field The equation (5 1) undergoes Hopf bifurcation 
(HB) if a complex pair of eigenvalues of the matrix 1 D x f(x, A)|( ro ,A 0 ) become purely imaginary 
A periodic orbit is, thus, bom It may be stable or unstable A periodic orbit is characterised 
by a point x* having a time period T such that the trajectory d> £ (x) emanating from x* returns 
back to x* in T time The system of equation (5 1) can be integrated for fixed A from 0 to T to 
yield a point x such that x(T) is equal to <j> ? (x*) 


lr The matrix D x f(x,X) = {§£} 
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In order for x* to lie on the periodic orbit and T to be time period 

x* — = F(x*, T) = 0 (5 2) 

The equation (5 2) is a system of n nonlinear equations in (n + 1) unknown variables At every 
iteration one component of x is fixed such that it lies on the periodic orbit, thereby reducing 
the number of unknowns by one The stability of a periodic orbit can be determined by the 
eigenvalues of the solution of variational equations, which can be calculated by appending the 
variational equations to the equations describing the dynamics 

x 

$ 

where $ t (x 0 ,t 0 ) = D xo <f> t (x 0 ,to) 

The above (n 2 +n) equations are integrated over the time period for initial conditions x(t 0 ) = 
x to and i>t 0 = / and the eigenvalues of the solution of variational equations are determined These 
eigenvalues are called Floquet multipliers A periodic solution is stable if all the eigenvalues lie 
within the unit circle except one that will be on the unit circle The stability of a periodic 
solution changes when some Floquet multiplier crosses the unit circle when A is varied If one 
multiplier which was earlier within the unit circle passes through the unit circle at —1 on the 
real axis, the original stable periodic orbit becomes unstable and a branch of periodic solutions 
with a two fold period branches off at this point With respect to the orientation of parameter 
variation the branching can be either supercritical or subcntical This bifurcation is usually 
called as period doubling bifurcation (PDB) If the original branch of the periodic solution is 
unstable, the new branch of the periodic solutions with two fold period will also be unstable 
Second important bifurcation takes place when a multiplier intersects the unit circle through 
+ 1 This case corresponds to a limit point (LP) when two periodic solutions coalesce If a 
stable periodic orbit coalesces with an unstable periodic orbit, it is called cychc fold bifurcation 
(CFB) 

Another bifurcation from a branch of periodic solution occurs when two complex conjugate 
eigenvalues intersect the unit circle The originally stable branch of periodic solution becomes 
unstable and a stable or unstable torus (TRB) may appear at the bifurcation point Torus is a 
quasipenodic motion However, if the periodic orbit is unstable, the torus will be unstable 


/(*) 
D x f(x)$ J 


(5 3) 
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5.3 Chaos 

5.3.1 Route to Chaos 

Chaos is seemingly a random behaviour in a deterministic system One fundamental character- 
istic of chaos is its extreme sensitivity to initial conditions Given two different initial conditions 
arbitrarily close to each other, the trajectories emanating from these two points diverge from 
each other with a wide separation. The transition to chaos may occur m one of the following 
four different ways 

5.3 1 1 Period Doubling (Sub-harmonic Cascade) Route to Chaos 

The period doubling route to chaos is the most widely known and studied route to chaos In 
this, a periodic orbit undergoes period doubling bifurcation m that the one original periodic 
orbit becomes unstable and a new periodic orbit of two fold period is born A complete penod 
doubling cascade is an infinite sequence of period doubling bifurcations after which the dynamics 
become chaotic 

5.3. 1.2 Quasi-periodic route to Chaos 

This route of chaos typically follows the sequences of bifurcations given below as the parameter 
is monotomcally varied 

Steady state =>■ Periodic Quasi-periodic (T 2 ) Quasi-periodic (T 3 ) =3* Chaos 

The main property of this route to chaos is that chaos is observed after a short sequence (finite) 
of bifurcation with continuous changes of the attractors 

5.3. 1.3 Intermittent Transition to Chaos 

Generally speaking, a signal (the value of a variable changing in time) is said to be intermittent 
if it is subject to infrequent variations of large amplitude On this route, a limit cycle becomes 
unstable as the parameter increases through a critical value. At the critical value a stable and 
an unstable limit cycle coalesce For very small increase of the parameter past its critical value, 
the same cycle persists but is occasionally and rarely interrupted by busts which are not small 
As the parameter is further increased, the busts occur more frequently but do not change much 
in magnitude, thus the dynamical system becomes chaotic 
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HOMOCLINIC ORBIT 

of infinite period 

Figure 5 3 Blue sky bifurcation 


5.3 1.4 Subcritical Instability 

On this route a familiar attractor, l e a fixed point, periodic or quasipenodic attractor, becomes 
unstable as the bifurcation parameter is slowly varied through a critical value The system jumps 
rapidly to a strange attractor which is not a continuous extension of the familiar attractor 

5.3.2 Blue Sky Bifurcation 

Thompson and Stewart [38] refer to a type of global bifurcation in stable space involving the 
discontinuous disappearance of a limit cycle as a blue sky catastrophe (Fig 13 5, page 271 of 
ref [38]) and reproduced as Fig 5 3 This is not a local bifurcation However, it does possess the 
feature of the saddle-node bifurcation, namely the disappearance of the attractor by collision 
with a saddle, in this case a saddle equilibrium point As seen from Fig 5 3, the blue sky 
bifurcation for a periodic orbit is the sudden disappearance of a limit cycle through collision 
with a saddle equilibrium point Prior to the critical parameter value A c at which the collision 
occurs, a saddle equilibrium coexists with the limit cycle At A c the limit cycle and a branch 
of both the stable and unstable manifolds of the saddle point coincide, forming a homochmc 
connection Past A c , the limit cycle no longer exists Collision with a saddle fixed point as 
a typical mechanism by which a limit cycle can abruptly vanish from state space This blue 
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sky bifurcation can take place m two forms, the disappearance of a stable limit cycle and the 
disappearance of an unstable limit cycle Moreover, it is a global, discontinuous and catastrophic 
bifurcation It also serves as a prototype of a blue sky bifurcation for a strange attractor or a 
boundary crisis of a strange attractor 

One type of global bifurcation involving a strange attractor is the sudden death of the 
attractor [120] Such a blue sky bifurcation occurs commonly for strange attractors of differential 
equations Like the blue sky bifurcation for a limit cycle, a chaotic blue sky bifurcation involves 
a collision with an object of saddle type and is analogous to the blue sky disappearance of a limit 
cycle as discussed above The chaotic blue sky bifurcation is also known by another term, namely 
the boundary crisis of a strange attractor The term crisis was introduced m [16,20] and applies 
to sudden qualitative changes m strange attractors with quasistatic changes m parameters A 
crisis involves the sudden destruction of a strange attractor through collision 

5.3.3 Lyapunov Exponents and Dimension 

Lyapunov exponents (LE) characterise the asymptotic orbital instability of dynamic systems 
As n th order dynamical system has n Lyapunov exponents Each exponent A(c, v)is a real number 
and is a measure of the divergence rate of trajectories m the direction of a vector v The positive, 
zero and negative Lyapunov exponents denote stretching, marginal and contracting directions 
in the phase space Some of the informations derived from Lyapunov exponents are as follows, 

(l) For conservative systems the sum of all Lyapunov exponents is zero 

(n) The attractor is a stable fixed point if all the LEs are negative 

(m) If all LEs are negative except one which is zero, system undergoes periodic motion More 
than one zero (other negative) indicate quasiperiodic motion or torus for multidimensional 
cyclic motion 

(iv) If one of the LEs is positive, one being zero and all other negative, the system undergoes 
chaotic oscillation Two positive, one zero and other negative LEs indicate the presence 
of Superchaos 

The Lyapunov exponents are defined as [59] 

A, = hm j In | m,(f) |, i = 1, , n (5 4) 
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where m,(t) are the eigenvalues of $ t (z 0 , t 0 ) 

The obvious way to find the Lyapunov exponents [59] is to integrate $ t (x 0 , <o) for T time 
and if T is large, the LE can be approximated as 


A, » ^ In | m,(t) |, i=l, ,n 


(5 5) 


However, there are two problems with this approach First, for a chaotic trajectory, A, > 0 and 
$ t (xo,io) 1S unbounded as t — > oo. Second, as t — * oo, each column of $ t (xo, io) tends to line 
up with the eigenvector associated with the largest eigenvalue m 1 (f) It follows that $<(x 0 , t 0 ) 
is ill conditioned and its eigenvalues may not be reliably computed 

A more robust algorithm [59] is based on the concept that almost every perturbation m the 
linearised system evolves on the average as Thus to find Ai , choose a time interval T > 0, 

an iteration count K : and any initial condition u^°\ normalised to 1 Calculate the solution to 
variational equation for t = kT normalising every T seconds as follows 


= 8x{kT,u( k -'\{k- 1)T) 

(5 6) 

r( fc ) 

u<*) = 

||x(‘)|[ 

(5 7) 

for k — 1 , , K 



where 6x(f, u, r) — $ t _ T (^ T (x 0 ),r)u is the solution at time t starting from initial condition u at 
time t Then, 


5x(M><°>,0) = ||x<'>|| . !|i“’||u (t| 


(5 8) 


and for K large enough 

A, » jt>Px(M><°',0)]| (5 9) 

® Alnjll|x w || (5 10) 

Kl k-l 

® Ax>n* (t| n (5ii) 

K1 fc=i 

Similarly, the area spanned by almost any n-linearly independent vectors evolves in the linearised 
system on the average as e< Al+ +A ")‘. Hence, (A x + . + A„) can be found by following the time 

evolution of almost any n-linearly independent vectors (x l5 ,x n ) However, as the evolution 
goes on, these vectors xi, ,i n ) tend to align in the same direction This effect can be avoided 
by orthonormahsing every T time using modified Gram-Schmidt orthonormalisation [59] At 
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each step in the iteration, the area is given by ||uji^|| where ||v|| is the orthonormalised 


vector If k is large 

A,+ . + A n = -rLf; In ||»W|| Ill'll (5 12) 

Kl k=l 

= + PrX>ll>4‘ ) ll (513) 

K1 k=l K1 k=l 

= Ai + + A n _! + — ^lnv^K (5 14) 

An = (5 15) 

K1 k= 1 


Thus by following the evolution, of volume elements with dimensions from 1 to ra, each of the 
Lyapunov exponents can be determined 

Let Xi > X 2 > A n be the Lyapunov exponents of an attractor of a continuous time 
dynamical system Let j be the largest integer such that X\ + + X 3 > 0 The Lyapunov 

dimension as defined by Kaplan et al [59] can be estimated 

« Ai T A, 

(516) 

If no such j exists, as is the case for a stable hyperbolic equilibrium point, Dl is defined to be 
0 For an attracting limit cycle, the generic situation is (Ai = 0 > A 2 > > A n ) Thus, the 

Lyapunov dimension of a generic attracting limit cycle is 1 Similarly, the Lyapunov dimension 
of a generic attracting K-penodic behaviour is K If the attractor is chaotic, Di is almost always 
a non-integer In a three-dimensional chaotic system with Lyapunov exponents (A + > 0 > A_, 
the Lyapunov dimension is given as 

(517) 

5.4 System Model 

For the studies in this chapter, the transmission lines have been modelled by their lumped 7r 
equivalent The generators and load models have been taken same as those described m section 
4 4 2 in Chapter 4 The models for OLTC and SVS utilized in the present work are presented 
below 
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5.4.1 OLTC Model 

The dynamic characteristics of tap changer and the operating practices vary from country to 
country [168] There is always a small voltage deadband within which tap changer is not 
activated 

In the present work, the dynamics of OLTC has been represented by a discrete time model 
which changes the tap setting n at each time step (k + 1) as follows [98] 

{ n(kAt) + An if V > V° + e 

n(kAt) if \V — V°\ < e (5 18) 

n(kAt) - An if V < V° - e 

with n mtn < n(k + 1) < n max 

where V = controlled secondary voltage 
V° = corresponding set point 
e = dead band 
An = OLTC step size 
At = time interval 
k — an integer time step 

5.4.2 SVS Model 

SVS has been modelled by a first order time delay block as shown in Fig 
descnbmg the dynamics of SVS can be written as 
d R 

T 3V3 — = I< 3V3 {V Te} - V) - B 

with B m xn — B ^ B max 

where B = shunt susceptance 
I< 3V3 = SVS gain 
Vref = reference voltage 
T 3V3 = SVS time constant 

The action of limiter has been realised by using tank function as discussed in Chapter 4 (equation 
(4 51)) 


5 4. The equation 

(5 19) 
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Figure 5 4 SVS Model 

5.5 Case Studies 

The studies were conducted on four sample systems viz a 3-bus system with SVS ( Appendix-C, 
a 4-bus system with OLTC (Appendix-D, a 9-bus system (Appendix-E) and 19-bus UPSEB 
system (Appendix-H) The systems m Appendices-C and D are derived from that in the ref [66] 
A transmission line between the two generators is added and through network reduction, the 
terminal node of the finite generator at bus-1 was eliminated The data for the reduced network 
retaining the internal node of the generator # 1 are given in Appendices-C &; D The studies 
were performed on HP9000/735 computer systems and assisted by AUT086 software [34] and 
INSITE software [59] 

5.5.1 3-Bus System with SVS 

The generator at bus-3 was taken as infinite bus and the generator at bus-1 was represented 
by swing equations (4 26) and (4 27) described in Chapter 4 The effect of SVS dynamics were 
studied using its models given in section 5 4 2 The load at bus-2 was assumed to be consisting 
of two parts, one being the static load of polynomial type and the other being an induction 
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PRRRMETER 


Figure 5 5 3-bus system with SYS bifurcation diagram 


motor load The real and reactive power components of static load were taken as follows 


P, = Pi + PkV"' 

Qs=Q\ + Q v V* 

where P x , Q x = constant real and reactive load demand 
Pv , Qv = coefficients of voltage dependent terms 
TJ, £ = exponents 


(5 20) 
(5 21) 


The model of induction motor considered in this study is the same as described by the equations 
(4 34) and (4 35) in Chapter 4 and also in ref [40] 

The system shown m Appendix-C was considered for the present study The bifurcation 
diagram has been drawn as shown in Fig 5 5 considering Q x as the parameter There are two 
Hopf bifurcations, subcritical HB at Q x = 14 4548 and supercritical HB at Q x = 14 4782 and a 
saddle-node bifurcation for <^i = 14 479395 Figs 5 6 and 5 7 show the movement of real part 
of complex eigenvalue and the largest real eigenvalue respectively After supercritical HB the 
critical mode corresponds to the real eigenvalue and causes the saddle-node bifuication 
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Figure 5 8 3-bus system with SVS participation factors of complex eigenvalue 

The participation factors of these two eigenvalues are plotted in Figs 5 8 and 5 9 for different 
values of parameter It is seen from these figures that in the Hopf bifurcation window, the 
dynamics is dominated by the generator states (<5,w) and for saddle-node bifurcation, the load 
bus voltage V dominates the dynamics 

Fig 5 10 demonstrates the presence of two chaotic regions, one between Q i = 14 44571 and 
14 45677 and the other between Qi = 14 45818 and 14 45930 For Q\ = 14 45592, the Lyapunov 
exponents are found to be (0 267, 0 000, —2 529, —63 380, —63 388) For Q x < 14 44521, the 
chaotic attractor disappears as a result of boundary crisis (possibly because of collision with 
saddle kind object like the saddle fixed point or the unstable limit cycle) Since for Q! < 
14 44571, the fixed point is the only stable attractor, the trajectories settle down to this stable 
fixed point, for Q i < 14 44571 (Fig 5 10) It is important to notice that the inclusion of 
SVS dynamics causes the flow to be overall contracting in the Hopf bifurcation window and 
no collapse of voltage or angular instability is observed However, for K sva = 20 and rj = 13, 
£ = 1 6 and T 3VS = 0 02 see the dynamics of the system even after consideung SVS dynamics 
is accompanied by voltage collapse due to boundary crisis Fig 5 11 shows the bifurcation 
diagram for this case, the system looses stability in the Hopf bifurcation window as a result of 
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Figure 5 9 3-bus system with SVS participation factors of real eigenvalue 
expanding flow The Lyapunov dimension of the attractor is 2 267 at Q x = 14 45592 

5.5.2 4-Bus System 

The model of generator at bus-1 and loads were taken same as m the study for 3-bus system 
with SVS The OLTC was modelled using equation (5 18) m section 5 4 1 The study on sys- 
tem shown in Appendix-D was carried out considering Q j (equation (5 20) as the bifurcation 
parameter while calculating the stationary points for the different values of Qi, three critical 
points were observed as shown in the bifurcation diagram of Fig 5 12 These include two Hopf 
bifurcation points, subcntical HB at Qi = 8 1650 and supercritical HB at 8 3199 and a saddle- 
node bifurcation at Qi = 8 3210168 At Qj = 7 8700 the OLTC hits its minimum limit (0 9 
p u ) and remains constant for further increase in Qi Fig 5 13 shows the movement of the real 
part of the complex eigenvalue and Fig 5 14, the movement of the largest real eigenvalue It is 
seen from the Figs 5 13 and 5 14 that in the stable region before the saddle-node bifurcation 
point the real eigenvalue is critical, while in the stable region before subcntical HB, the complex 
eigenvalue is critical 

The eigenvalue sensitivity analysis has been performed using the technique given in lef [72] 
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Figure 5 14 4-bus system with OLTC movement of real eigenvalue 

and the participation factors for complex eigenvalues and real eigenvalue are shown in Figs 5 15 
and 5 16 for dilferent values of Q\ It is interesting to note that the real eigenvalue is dominated 
by the load bus voltage V while the complex eigenvalue is mainly associated with (S,u) states 
Through numerical simulation, it is observed from Fig 5 17 that at subcntical HB a unstable 
limit cycle is born which surrounds the stable fixed point until Qi = 7 9163 The Floquet 
multipliers were computed and all of them were found to be within unit circle except one which 
was found to be outside the unit circle Thus, the periodic orbit is type-1 unstable It was 
observed during the numerical simulation, that the Floquet multiplier which was outside the 
unit circle enters the unit circle through +1, thus giving rise to a CFB at Q\ = 7 9163 The OFB 
arises when a stable and unstable limit cycle coalesce Thus it can also be called as saddle-node 
or fold bifurcation of the limit cycle Therefore, at Qi > 7 9163, a stable periodic orbit exists 
The stable periodic orbit undergoes a PDB at Qi = 7 9165 A Floquet multiplier crosses the 
unit circle through - 1 The original stable periodic orbit becomes unstable and remains so until 
Qi = 8 31057 At this point, a Floquet multiplier which was outside the unit circle enters the 
unit circle through -1 The periodic orbit gains stability at Qi = 8 31057 and remains so until 
supercritical HB at Qi = 8 3199 This periodic orbit undergoes further PDB leading to chaotic 
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Figure 5 15 4-bus system with OLTC participation factors of complex eigenvalue 



PRRRMETER 


Figure 5 16 4-bus system with OLTC participation factors of real eigenvalue 
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VOLTAGE 



Figure 5 17: 4-bus system with OLTC bifurcation diagram 

oscillations Starting from Q i = 8 3199 the bifurcation diagram using brute-force method (page 
210, ref [59]) has been shown in Fig 5 18 At Q x = 8 3067, the system looses stability and the 
simulation halts as a consequence of expanding flow i e divergence of the vector field becomes 
positive Figs 5 19 to 5 21 depict a typical simulation for Q\ = 8 30567 from the following 
initial conditions, (6 = 0 3149, u — 0 22016, 6 = — 0 00187, V — 0 87188) It is seen that the 
system voltage though oscillating remains fairly healthy for upto 27 6 seconds and then collapses 
as seen from Fig 5 19 The angle of load bus voltage shoots up at that point in time as shown 
in Fig 5 20 However, generator angle variation remains within acceptable range as shown in 
Fig 5 21 

The collapse of voltage takes place due to boundary crisis (blue sky bifurcation) at Q x = 
8 30567 In the boundary crisis, a strange attractor exists for parameter values upto the critical 
value, at which the collision takes place Subsequent to this value, the strange attractor no 
longer exists, but it leaves a signature of a transient chaotic motion The transient chaotic motion 
appears chaotic for a relatively long time (depending on the initial condition) and then suddenly 
experiences a sharp excursion either to another attractor or to infinity This excuision octuis 
through a tunnel in state space which necessarily follows the unstable manifold of the saddle 



CHAPTER 5 OBSERVANCE OF GLOBAL BIFURCATIONS 


115 



PRRRMETER 


Figure 5 18 4-bus system with OLTC chaotic dynamics 



TIME ( IN SECONDS ) 


Figure 5 19 4-bus system with OLTC time variation of V 
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Figure 5 22 4-bus system with OLTC time variation of V in chaotic region 

point or orbit with which the collision takes place The result is a discontinuous, catastrophic 
disappearance of the strange attractor 

Lyapunov exponents have been computed for Q i = 8 3067 LEs are found to be (0 331, 
0 000, -2 03, -33 615) One positive Lyapunov exponent confirms the presence of chaos Fig 
5 22 depicts the time variation of voltage in the chaotic region for parameter Q x = 8 3067 Fig 
5 23 shows the corresponding simulation projected m the phase onto load bus voltage and bus 
angle coordinates The Lyapunov dimension of the attractor was computed as 2 163 

Further it was observed that for for Q x > 7 916320 (i e after CFB) a stable periodic orbit is 
born which undergoes PDB at Q x = 7 916354 The stable periodic emerging as a result of PDB 
undergoes further PDB at Q x = 7 916370 Another PDB occurs at Q x = 7 916374 Since the 
bifurcations took place in very narrow range, the simulation could not be carried out to trace 
further bifurcations 

5.5.3 9-Bus System 

The generator and load model for this system was taken same as those considered in studies 
presented section 4 5 12 of the Chapter 4 The dynamics was studied by varying MVA load ( 
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Figure 5 23 4-bus system with OLTC plot of trajectory projected on V, 9 plane 

at bus-5 and power factor pf was assumed to be unity Two sets of studies were performed, one 
at the base case condition and the other being a contingency case considering the outage of the 
line between buses-4 and 5 

5.5 3.1 9— Bus System. Base Case 

The bifurcation diagram has been drawn in Fig 5 24 Some of the important bifurcations has 
been observed as follows 

• SNB at t = 3 598273 (point-3) 

• Subcntical HB at t = 3 577746 and 3 588495 (points-2 k 3) 

• TRB at e = 3 249934 and 1 037992 (points-6 & 13) 

• PDB at i = 3 135110,3 235849,3 042431 and 3 558070 (points-8, 9, 11 k 14) 

• LP of periodic orbits at £ = 3 070598,3 268290 and 3 022379 (pomts-7, 10 k 12) 

The numerical investigations reveal that unstable periodic orbits emerging at i = 3 5777 16 
and i = 3 588495 are of type-1 kind (l e one Floquet multiplier is outside unit aide) and at 
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Figure 5 24 9-bus system (base case) bifurcation diagram 

i = 3 249934 one periodic orbit bifurcates to unstable torus and the other one undergoes a PDB 
at l = 3 558070 For 3 249934 < i < 3 577746 there exists a periodic orbit of type-1 kind 
Since the boundary of region of stability of the stable equilibrium point is given by the union of 
stable m ani fold of the type-1 unstable equilibria and closed orbit [45], the region will shrink in 
presence of unstable periodic orbit having one Floquet multiplier outside unit circle 

The participation analysis was carried out to ascertain the participation of the states cor- 
responding to the critical eigenvalue causing the Hopf bifurcation The eigenvalues along with 
the participation factors at t = 3 577746 are shown in Table 5 1 It is found that the states 
associated with generator 2 are predominant in the critical mode Fig 5 2o shows a typical 
time variation of load voltage V 5 for l = 3 58 from initial condition (012 = ^3 = 0, 62 — — 0 2154, 
<5 3 = _o 26531, d 5 = -0 89383, V$ = 0 71238) This initial condition corresponds to stable 
operating point for t = 3 50 It is seen that on account of the subcritical Hopf bifurcation, the 
voltage at the load node will experience growing oscillations and collapses after about 135 16 
seconds 
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Table 5 1 9-bus system eigenvalues and state participation factors 


Eigenvalues 

Participation Factors 

^2, U>2 

$3, W3 


-1 2802 x 10 2 

8 19 x 10" 4 

5 26 x 10" 4 

5 82 x 10" 1 

4 19 x 10" 1 

-8 0958 

1 26 x 10" J 

4 99 x 10 -3 

5 03 x 10~ 3 

4 16 x lO' 1 

5 47 x 10' 1 

-3 5180 
±j 11823 

9 91 x 10" 2 

4 00 x lO" 1 

6 16 x 10' 4 

3 94 x 10" 4 

3 2445 x 10“ 7 
±j6 0707* 

3 95 x IQ" 1 

9 68 x 10~ 2 

1 29 x 10~ 2 

1 63 x 10- 2 


Voltage 
0 75 

0 7 

0 65 

06 

0 55 

05 

0 20 40 60 80 100 120 140 

Time (sec) 



Figure 5 25 9-bus system (base case) time vatiation of V 
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Figure 5 27 9-bus system (a contingency case), chaotic dynamics 

The unstable periodic orbit, for 1 076371 < £ < 1 116463, has one Floquet multiplier outside the 
unit circle Though there is a stable fixed point (operating point) for l 076371 < i < 1 116463 
the impact of this unstable periodic orbit is to decrease the stability region The stable periodic 
orbit emerging at i = 1 076371 undergoes a series of PDB between l — 1 115625 and l — 
1 171 130 leading to chaotic oscillations as shown in Fig 5 27 

For t = 1 14, Lyapunov exponents were found to be (0 135, 0 000, —0 032, —0 38, —6 657 
and —40 665) The presence of one positive LE confirmed the chaotic dynamics The Lyapunov 
dimension of the attractor is 3 284 Variation of the load voltage with time and its frequency 
spectrum are shown in Figs 5 28 and 5 29 The frequency spectrum in Fig 5 29 is continuous 
and broad band m nature 


5.5.4 19-Bus UPSEB System 

In this study, all the base case loads have been modelled as constant impedance type An 
additional load has been considered at bus-14 represented by load model-II described by the 
equations (4 36) and (4 37) in Chapter 4 The bifurcation diagram has been drawn as shown 
in Fig 5 30 considering t as the parameter and keeping the power factor unity The other 
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Figure 5 28 9-bus system (a contingency case) 


time variation of V 5 



freq Cfe), 40 c £-pcint FF7 


Figure 5 29 9-bus system (a contingency case) frequency spectrum 
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VOLTAGE 



Figure 5 30 19-bus system bifurcation diagram 

parameters of the dynamic load D , a , b and k were taken to be 0 05 (sec), 0 0, 0 0 and 0 1 (sec), 
respectively Some of the important points on the bifurcation diagram are given as below 

• Subcntical HB at I = 2 696254 

• Supercritical HB at I = 5 454979 

• SNB at l = 5 485668 

• CFB at £ = 2 519270 

• TRB at £ = 3 078993 and 3 911500 

The unstable periodic orbit for 2 519270 < l < 2 696254 is type-1 and will reduce the stability 
region of s e p At t = 2 519270, the Floquet multiplier which has been outside the unit circle 
enters the unit circle through +1, thus giving rise to a CFB The periodic orbit becomes stable 
and remains so until l — 3 078993 where a torus bifurcation occurs This torus breaks down into 
a chaotic attractor as I is further increased which disappears at I = 3 59421 due to boundaty 
crisis There is no attractor from I = 3 59421 to 3 911500 where a complex pair of Floquet 
multiplier enters the unit circle, thus giving rise to a stable periodic orbit 
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Figure 5.31 19-bus system chaotic dynamics 

The bifurcation diagram for the chaotic region is shown m Fig 5 31 It is observed that 
at places along the parameter there exist few periodic windows This is called mode-locking or 
phase-locking and is due to the fact that the frequencies of oscillations become rational, resulting 
in periodic behaviour The Lyapunov exponents and dimension were computed and are shown 
m Table 5 2 


Table 5 2 19-bus system Lyapunov exponents and dimension of chaotic attractor 


Parameter l 


Lyapunov ' 

Exponents 


Dimension 

3 1 

0 0306313 
-0 0536428 

0 00 

-0 0617678 

-0 001 
-183 165 

-10 0378236 
-201 391 

3 80984 

3 14646 

0 0300368 
-0 0565283 

0 0145376 
-0 126793 

-0 00 
-168 523 

-0 041033 
-200 27 

4 06265 

3 36451 

0 0237414 
-0 0473544 

0 0156513 
-0 223015 

0 00 

-165 617 

-0 0454897 
-199 915 

3 86597 
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5.6 Conclusion 

From the studies carried out on the three sample systems, in this chapter, to explore the existence 
of various types of global bifurcations and chaos, following observations can be made 

(i) Different types of dynamic bifurcations such as PDB, CFB, TRB, LP and chaos can exist 
m some of the power system networks for variation of the parameter even in their normal 
operating range 

(n) Some of the global bifurcations observed, possibly for the first time m power system model, 
include the bifurcation of periodic orbit to torus (TRB) and limit points of unstable to 
unstable periodic orbits 

(m) Chaos has been observed in all the three sample systems Chaos via PDB route has been 
observed in two of the systems and quasipenodic route m one of the systems Chaos via 
torus breakdown has been observed for the first time in a practical 19-bus Indian system 

(iv) The chaotic oscillations produce harmonics of all the frequencies with low frequencies (less 
than about 10 Hz) being dominant and a continuous broad band frequency spectrum is 
observed 

(v) Boundary crisis (blue sky bifurcation) and subcritical Hopf bifurcation cause voltage col- 
lapse in some of the systems 

(vi) By adjusting the SVS gain parameter the boundary crisis and hence the voltage collapse 
can be avoided 



Chapter 6 


Control of Dynamic Bifurcations and 
Chaos 


6.1 Introduction 

Presence of dynamic bifurcations and chaos has been shown to exist in several models of power 
system in Chapter 5 Though a chaotic oscillation is bounded, the system exhibits seemingly 
random behaviour and poses serious problems from operation view point Some of the implica- 
tions of these nonlinear phenomena already observed are the reduction in stability margin/region 
of attraction especially due to subcritical Hopf bifurcation, introduction of harmonics especially 
low frequency harmonics during chaotic oscillation and system voltage collapse due to subcritical 
Hopf bifurcation and boundary crisis 

Thus, these phenomena are unwanted in the system and need to be eliminated One of 
the easiest ways to remove these bifurcations in the system is by restricting the parameter to 
enter this region This may limit the transmission usable capacity resulting into reduced system 
loadability One of the efforts for stabilising unstable periodic orbit of a chaotic attractor is 
due to Ott et al [77], who applied the computed values of an externally manipulable system 
parameter This concept has been used by Bandyopadhyay et al [124] to stabilise a chaotic non 
isothermal CST reactor by controlling the bifurcation parameter or some other parameter to 
yield almost a periodic behaviour It has been demonstrated by Tan et al [154] that by varjing 
the damping, the dynamic bifurcation and chaos can be controlled Wang et al [157] employed 
a feedback law to control the stability and amplitude of the bifurcated solution using geneiatoi 
speed deviation signal 
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With developments in high power electronics, recently new devices, popularly called FACTS 
(Flexible AC Transmission System) devices, for control of voltage and power have been developed 
[108, 137, 138] The prominent among them are controllable series capacitors (CSC), static phase 
angle regulator (PAR) and static VAR compensators (SVC) 

In this Chapter, the studies have been conducted on the three sample systems consisting 
of 3-bus, 9-bus and 19-bus to explore the effectiveness of these new devices for controlling the 
dynamic bifurcations A first order model of these devices have been considered and the feedback 
control signal has been decided based on participation factor analysis 


6.2 Control of Bifurcation 

Bifurcation control deals with design of control inputs to modify the bifurcation characteristics 
of a parametensed system The control signal can take many forms, including open loop control 
or static or dynamic feedback The objective may be stabilisation and/or delay of a given 
bifurcation, reduction of the amplitude of bifurcated solutions, introduction of bifurcation at 
desired parameter values or a combination of these 

Consider a one-parameter family of nonlinear control systems 

x = f(x,u,n) (6 1) 

where x € R n is the state vector, u is the scalar control and ft 6 R 1 is the bifurcation parameter 
The vector function / is assumed to be smooth m x, u, and fi Let equation (6 1), with the 
control set to zero undergoes either a subcritical or supercritical HB from a fixed point x 0 (/i) at 
the critical parameter value /i = fi c 

Consider the case in which the HB is subcritical l e for n < /j. c , the Xo(fi) is stable and a 
unstable periodic orbit surrounds it Then, for n > /i c the system trajectory will diverge If 
the control u can be designed so as to render the HB to be supercritical, this situation will 
be remedied to a certain extent l e the same system trajectory will now converge to a small 
amplitude periodic solution Since this periodic solution is near the nominal equilibrium point, 
divergence of the system to a distant (and unacceptable) operating mode can thus be prevented 
Ref [123,161] present details of the laws 

In general, a static state feedback using the control law 

(6 2 ) 


u = u(x — x 0 (/l)) 
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designed with reference to the nominal equilibrium path xq ( fi) of equation (6 2) will not only 
affect stability but also the location and stability of other equilibria The mam disadvantage 
stems from the fact that system performance is often degraded by operating at an equilibrium 
which differs from the one at which the system is designed to operate 

For these reasons, the dynamic state feedback is recommended [123,161] This is achieved 
by incorporating high pass filters known as washout filters into the structure of the allowed 
controllers The use of washout filters is common in power system stabilisers A washout filter 
is stable high pass filter with transfer function 

= < 63 > 
Specifically, m equation (6 1), for each system feedback state variable x t7 i = 1,. n a 
washout filter is introduced The value of d is a positive parameter This helps m preservation 
of equilibrium for variation of parameter m its operating range However, m the present study, 
only the static feedback is considered Equilibrium preservation is guaranteed since the feedback 
signals used are the generator speed deviation and generator acceleration power which are zero 
at equilibrium 

There are many scenarios in which bifurcation can result into a chaotic invariant set Various 
routes to chaos have been discussed m Chapter 5 An important aspect in controlling chaos 
is that they may be suppressed by controlling a bifurcation in a given route to chaos The 
bifurcation control approach to this problem entails designing feedback control law which ensures 
a sufficient degree of stability for a primary bifurcation m such a scenario 

In order to control dynamic bifurcation and chaos in the power system networks, the FACTS 
devices represented by a first order delay block has been utilised The participation analysis 
has been carried out to ascertain the appropriate error signal to these devices The generator 
speed deviation and acceleration power has been considered as error feedback signals because 
their steady state value would be zero and for this reason the washout filters have not been 
considered in this study 


6.3 Model of FACTS Devices 


The FACTS devices (CSC, PAR and SVC) have been represented in this study by a first ordei 
time delay model Their output Y is governed by the following differential equation 



I< a (X r - X rn ) - Y 


(b 1) 
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3.n.d Y mjn ^ Y ^ Y ma x 

where T s is the time constant, K s the gam of the controller, X r is the reference value of the 
input quantity such as reference speed of the generator, or the mechanical power input to the 
generator X m is the measured value of the feedback signal such as generator actual speed or 
real power output of the generator The output Y is the series capacitance X csc for CSC, phase 
angle shift <j> p for PAR and the shunt sustenance b sh for SVC 

6.4 Case Study 

The studies have been performed on HP-9000/735 computer system using AUT086 software [34] 
The time constant T s for all the three FACTS devices was taken as 0 01 sec m each of the studies 
The effect of controller gam variation was observed on the dynamic bifurcations As shown m 
subsequent sections, the generator states (S,oj) were found to be dominant using participation 
analysis m the Hopf bifurcation window in each of the systems studied Hence, the generator 
speed deviation u> and the generator acceleration power P a (= Pm — P ge ) was tried out as error 
signals to the controller of the FACTS devices These signals were also selected due to the fact 
that they do not affect the steady state performance and contribute to the damping of power 
system [147, 156, 171, 180] Participation of state variables in the critical mode were computed 
using the technique described m ref [72] to decide the input signal to the FACTS devices 
The studies were performed on the three sample systems namely 3-bus, 9-bus (base case and 
a contingency case) and 19-bus systems Results considering the effect of each of the FACTS 
devices at a time are presented below 

6.4.1 3-Bus System 

The data and network topology for this system have been shown in Appendix-B The generator 
at bus-3 was taken as infinite bus and the generator at bus-1 was modelled by its swing equation 
as described in section 4 4 1 of Chapter 4 A dynamic induction motor load model due to 
Waive [40] and as described in section 4 4 2 of Chapter 4 was considered along with constant 
static P and Q load at bus-2 Without considering the effect of FACTS devices, the bifurcation 
diagram has been drawn in Fig 6 1 treating (static) reactive load power Qi as the bifurcation 
parameter Some of the important bifurcation points observed in Fig 6 1 are listed below 


• Subcritical HB at Qi = 6 272709 
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Figure 6 1 Bifurcation diagram for 3-bus system 

• Supercritical HB at Q\ = 6 542023 

• CFB at Qt = 6 181948 

• PDB at Q x = 6 215310 and 6 538772 

• SNB at Q x = 6 543211 

The locus of HB points under variation of mechanical input power to the generator P\j is 
shown in Fig 6 2 Fig 6 2 reveals that if Pm is reduced, the two HB points move closer and 
for Pm = 0 425, they coalesce and disappear This, however, results into under utilisation of 
generator and transmission capacity 

Six sets of studies were conducted to observe the effect of each controllable component, one 
at a time — three with the generator speed deviation as the error signal and the other three 
with generator acceleration power as the error signal 

The CSC and PAR were sited in the line connecting bus-1 and 2 and the SVC was located at 
bus-2 as shown in Fig 6 3 Participation analysis was carried out to ascertain states predominant 
in the critical mode and it was found that the generator # 2 states (S,uj) dominate the HB 
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Mech Power Input Pm 



Figure 6 2 Effect of Pm on Hopf bifurcations 


window Therefore, the error signals were derived from this generator The following limits 
were assumed for these devices These limits, however, will be encountered only in the time 
domain simulation studies involving disturbances 

X™ x = 0 1 p u <f>™ x = 0 5 rad = 1 0 p u 

X™ n = -0 1 p u 4>™ x = -0 5 rad 6^ n = -1 0 p u 

Figs 6 4 to 6 6 present the locus 1 of HB points considering CSC, PAR and SVC respectively with 
generator speed deviation Fig 6 4 shows the locus of CFB and HB points for gam variation of 
CSC When the gain of controller is increased from zero, the two HB points move closer Also 
the CFB moves closer to the HB point and for a gain value of 0 006, the CFB coalesces with the 
subcntical HB and disappears When the gain is further varied, the two HB points move closer 
and coalesce and disappear for a gain value of 0 0075 The strange attractor and HB window 
disappear and the dynamics become stable Similar effects of PAR and SVC were observed on 
the system dynamics From Figs 6 5 and 6 6, it is seen that HB window can be completely 
eliminated by PAR and SVC if the gam value of PAR is greater than 0 04725 and that of SVC 
is greater than 0 25 


1 It shows the position of the two Hopf bifurcations as the gain is varied 
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Gam of CSC 



Figure 6 7 3-bus system CSC gam vs Q\ with P a as error signal 

The studies were also conducted with generator acceleration power as error signal to these 
devices Figs 6 7 to 6 9 show the locus of HB points It can be seen from these these figures 
that if the gam of CSC or PAR is greater than 0 72 and 0 38 respectively or that of SVC is less 
than —0 7 the HB window and the strange attractor can be eliminated 

6.4.2 9-Bus System 

6.4 2.1 Base Case 

The data and network topology for this system have been shown in Appendix-E The generator 
at bus-1 was taken as the infinite bus and the generators at bus-2 and 3 were represented by their 
swing equations (4 26) and 4 27) given in Chapter 4 The stator transients were ignored All the 
base case loads were represented by their constant impedance and in addition a dynamic load 
model as considered in the study presented in section 4 5 1 2 of Chapter 4 was assumed to be 
present at bus-5 The dynamics were studied by varying loading parameter £ and power factor 
p f was assumed to be unity Some of the important bifurcations observed from the bifurcation 
diagram in Fig 5 24 of Chapter 5 are listed below 


• SNB at l = 3 598273 
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Gain of PAR 



Figure 6 8 3-bus system PAR gam vs Q i with P a as error signal 


Gain of SVC 
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Figure 6 9 3-bus system SVC gain vs Q\ with P a as error signal 
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• Subcntical HB at l — 3 577746 and 3 588495 

• TRB at l = 3 249934 and 3 037992 

• PDB at t = 3 135110,3 235849,3 042431 and 3 558070 

• LP of periodic orbits at i = 3 070598,3 268290 and 3.022379 

From participation analysis carried out in the studies presented in section 5 5 3 1 in Chapter 5 
(Table 5 1), it is observed that the states predominantly responsible for causing Hopf bifurcations 
were the states associated with generator #2 It is seen from Table 5 1 in section 5 5 3 1 of 
Chapter 5 that critical mode is extremely sensitive to (S 2 , u> 2 ) states (l e Generator # 2) Hence, 
the input quantity to the controller was derived from Generator # 2 Accordingly the CSC and 
PAR were placed in line between buses-7 &: 5 and SVC at load bus-5 The following limits were 
assumed for these devices 

= 0 1 p u <t>™ ax = 0 5 rad b™£ x = 1 0 p u 

X™ n = -0 1 p u <j>™ ax = -0 5 rad b^ n = -1 0 p u 

Fig 6 10 shows the locus of Hopf bifurcation for variation of the controller gam with u 2 as 
input signal to CSC It is found that for K C3C > 0 0125 with u 2 signal the Hopf bifurcation come 
closer, coalesce and disappear A similar study was made with P a signal and it was found (Fig 
6 11) that for K csc > 0 1, it is possible to eliminate the Hopf bifurcations completely. With the 
application of u> 2 signal to PAR, it turned out (Fig 6 12) that for K paT < —0.125, the Hopf 
bifurcation phenomena could be eliminated while the application of P a signal to PAR did not 
help as is evident from Fig 6 13 

The study was also performed with P a signal to SVC at bus-5 and it was found that Hopf 

bifurcations could be eliminated (Fig 6 14) if K svc > 1 25 However, with u> 2 signal, it was 

found that the two Hopf bifurcations are persistent and SVC could not eliminate them as shown 
in Fig 6 15 
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Figure 6 10 9-bus system (base case) CSC gam vs £ with u? 2 as error signal 



Gam of CSC 



Parameter £ 


Figure 6 11 9-bus system (base case) CSC gam vs £ with P aj as error signal 
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Gain of J 



Figure 6 14 9-bus 
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Fig 6 16 shows a typical time variation of load voltage V 5 for £ = 3 58 from initial condition 
(w 2 = w 3 = 0, S 2 = -0 2154, S 3 = -0 26531, 6 S = -0 89383, V s = 0 71238) This initial 
condition corresponds to stable operating point for £ — 3 50 The simulation has been performed 
for two cases, one without considering the effect of CSC and the other with CSC present in the 
line with a gam (/f CJC ) setting of 0.15 In this case, P a was used as input signal to CSC It was 
found that the application of CSC could avert the system voltage collapse 

6. 4. 2.2 A Contingency Case 

In this study the line 5-4 was taken out and the existing static load of (125 + j50) MVA was 
reduced to ( 100+j30) MVA This change in loading was done as the load flow for the contingency 
case was not converging 

Without considering the effect of controllable components, the bifurcation diagram for •vari- 
ation in parameter £ has been drawn as shown m Fig. 5 26 m Chapter 5 Some of the important 
points on the bifurcation diagram are as below 

• SNB at £ = 1 195323 

• Subcntical HB at £ = 1 116463 

• Supercritical HB at £ = 1 192705 

• CFB at £ = 1 076371 

• PDB at £ = 1 115625 and 1 171130 

• LP of fixed points at £ = 0 2351554 and 0 6105376 

From the participation analysis, it was found that the states predominantly responsible for 
causing Hopf bifurcations in this case also were the states associated with the generator # 2 
The CSC and PAR were placed in line 7-5 and SVC at bus-5 The working limit of these 
components are assumed to be the same as given for the base case study 

The speed deviation and acceleration power of generator $ 2 have been taken as the input 
signal to the controller of FACTS devices The locus of HB points are plotted in Figs 6 17 to 
6 22 It is found that by suitable choice of gains for these devices, the strange attractor and the 
HB window can be eliminated by either of the two signals 
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Voltage 



Figure 6 16. 9-bus system (base case) time variation of V 5 


Gam of CSC 



Figure 6 17 9-bus system (a contingency case) CSC gain vs i with u; 2 as error signal 
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Gam of CSC 



Figure 6 18 9- bus system (a contingency case) CSC gam vs I with P a2 as error signal 
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Figure 6 19 9-bus system (a contingency case) PAR gam vs t with u ; 2 as error signal 
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Figure 6 20 9-bus system (a contingency case) PAR gam vs £ with P a2 as error signal 
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Figure 6 21 9-bus system (a contingency case) SVC gain vs £ with uii as error signal 
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6.4.3 19-Bus UPSEB System 

This system is described m Appendix-H The generator at bus-1 was treated as infinite bus and 
those at buses-2, 3 and 4 were represented by swing equations (4 26) and (4 27) as described in 
Chapter 4 All the base case loads were assumed to be constant impedance type In addition, a 
dynamic load model as described by equations (4 36) and (4 37) in Chapter 4 was assumed to 
be present at bus- 14 The nominal MVA demand l was taken as bifurcation parameter and the 
power factor was assumed to be unity The other parameters of the dynamic load was taken to 
be the same as in section 5 5 4 of Chapter 5 The bifurcation diagram as shown in Fig 5 30 of 
Chapter 5 was drawn considering l as the bifurcation parameter, and the important bifurcation 
points as given in section 5.30 of Chapter 5 are 

• Subcntical HB at I = 2 696254 

• Supercritical HB at I = 5 454979 

• SNB at l = 5 485668 

• CFB at I = 2 519270 

• TRB at i = 3 078993 and 3 911500 

The CSC and PAR were sited in the line between buses-5 and 7 and SVC at bus 14 The 
following limits were assumed for these devices 

X™ x = 0 002 p u <f>™ ax = 0 5 rad = 1 0 p u 

X™ n = -0 002 p u 4>p' n = -0 5 rad b^ n = -1 0 p u 

The participation analysis was carried out and the results are presented in Table 6 1 for i = 
3 0793 It is observed that the states predominantly responsible for causing the HB is associated 
with generators # 2 and 3 The CSC and PAR were placed on the line between buses-5 and 
7 and SVC at bus- 14 The error signals have been taken as the sum of the generator speed 
deviations of generator # 2 and 3 and the sum of acceleration powers of the generator # 2 and 
3 

Fig 6 23 shows the locus of Hopf bifurcation for variation of the controller gam with gen- 
erator speed deviation as error signal to CSC It is found that for K C3C < —0 0012 the Hopf 
bifurcation come closer, coalesce and disappear Control of dynamic bifurcations using CSC 
was tried with generator acceleration power as the error signal to the controller with its location 
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1 08 11 1 12 1 14 1 16 1 18 12 

Parameter I 

Figure 6 22 9-bus system (a contingency case) SVC gam vs I with P a2 as error signal 


Table 6 1 19- bus system eigenvalues and state participation factors 


Eigenvalues 

Participation Factors 

<?2, UJ2 

S3, U>3 

S4, U4 

#14> Vl4 

-2 3641 x 10+ 02 
±jl 3122 x 10 +O2 

6 3782 x 10~ 05 

6 3784 x 10~ 05 

2.9551 x 10~ 05 

2 9552 x 10’° 5 

3 2744 x 10 -05 

3 2748 x lO' 05 

5 0361 x 10“ 01 

5 0372 x lO" 01 

-1 6098 x 10~ 02 
±jl 5169 x 10 +O1 

5 1885 x 10-° 2 

5 1885 x 10"° 2 

5 2224 x lO -03 

5 2224 x 10~ 03 

4 4289 x lO" 01 

4 4289 x lO -01 

2 3651 x 10"° 5 

3 3128 x lO" 06 

-7 2877 x 10" 03 
± 3 1 2656 x 10 +m 

2 0340 x 10~ O1 

2 0340 x 10~ 01 

2 8736 x lO" 01 

2 8736 x 10~ 01 

9 2448 x lO' 03 

9 2447 x lO' 03 

2 1481 x lO" 05 

7 3717 x lO -06 

1 5435 x 10-° 3 
±j7 3118 

2 4467 x lO" 01 

2 4467 x 10~ 01 

2 0740 x 10“ 01 

2 0740 x 10“ 01 

4 7844 x 10~ 6i 

4 7843 x lO" 02 

1 5211 x 10" w 

2 2559 x lO" 04 
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Gam of CSC 



Figure 6 23 19-bus system. CSC gam vs £ with sum of and uj 3 as error signal 

considered at one of the lines between buses-5 and 7, 5 and 8 and 14 and 15 The study revealed 
that the HB window could not be eliminated m either of the cases 

Figs 6 24 and 6 25 show the locus of HB points using PAR with generator speed deviation 
and acceleration power as error signals It is observed that the HB window can be eliminated for 
K V ar < —0 00175 using speed deviation as the error signal and K par < —0 016 for acceleration 
power as the error signal 

Similarly Figs 6 26 and 6 27 reveal that SVC HB window can also be eliminated using SVC 
with either the speed deviation as error signal ( K svc > 0 0225 or the acceleration power as error 
signal (K 3vc > 0 125 


6.5 Conclusion 

From the studies performed in this chapter, the following conclusions may be drawn 

(l) On varying the controller gain, the subcntical HB becomes supercritical after a certain 
gain value with CFB disappearing Further variation of controller gain value may cause 
the other dynamic bifurcations to disappear yielding a stable steady state opeiating point, 
in the earlier Hopf bifurcation window 











CHAPTER 6 CONTROL OF DYNAMIC BIFURCATIONS AND CHAOS 


150 


(n) The CSC with generator speed deviation as error signal has been able to eliminate dynamic 
bifurcation and chaos in all the three base case system studies and also m the contingency 
case 

(in) The choice of appropriate input signal to the controller of FACTS devices is also very 
important as evident from the studies on 9-bus system (base case conditions) PAR could 
not eliminate the dynamic bifurcation and chaos with generator acceleration power (P a ) as 
error signal and the SVC could not eliminate the dynamic bifurcations with speed deviation 
as error signal However, PAR and SVC could eliminate the dynamic bifurcations when 
the choice of input signals was reversed Similarly, for 19 bus system, CSC did not work 
with P a as error signal 

(iv) The FACTS devices such as CSC, PAR and SVC can be effectively used for controlling 
the occurrence of dynamic bifurcation and chaos by proper choice of error signal and gam 
of the controller 



Chapter 7 


Conclusions 


7.1 General 

With growing interconnection along with economic and environmental pressures, the possible 
threat of voltage instability is becoming increasingly pronounced in pow'er system networks 
Researchers and utilities, all over the world, are devoting efforts to understand, analyse and 
develop strategies to cope up with the menace of voltage instability /collapse This thesis has 
attempted to make an iota of contribution m this direction. 

The thesis has addressed to both static and dynamic aspects of voltage stability Investiga- 
tions considering static aspects include determining point of voltage instability in load parameter 
space for both AC and combined AC-DC systems (Chapter-2) and maximising the voltage sta- 
bility margin through rescheduling of real and reactive power outputs of sources (Chapter-3) 
Dynamical studies have been conducted to explore the occurrence of various local bifurcations 
(Chapter-4), global bifurcations and chaos (Chapter-5) considering dynamical models of gen- 
erators, composite loads, induction motors, SVS, OLTC etc and also the control of dynamic 
bifurcations (Chapter-6) 

The aim of this chapter is to summarise the main findings of this thesis and provide sugges- 
tions for further research work Some of the main findings are given below 


7.2 Summary of Important Findings 

Chapter 2 has applied the concept of minimum singular value and condition number of load flow 
Jacobian for prediction of static voltage stability of integrated AC as well as AC-DC systems 
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From the works presented in this chapter, the following conclusions are drawn 

(1) The concept of minimum singular value and condition number of power flow Jacobian can 
be effectively used to predict the static voltage instability of integrated AC-DC networks 
as well as the AC systems having polynomial type voltage dependent loads 

(2) At static voltage instability point corresponding to system maximum loadability, the con- 
dition number approaches to infinity and the minimum singular value approaches to zero 
With increase in loading, an abrupt change in the minimum singular value and condition 
number of the load flow Jacobian is observed whenever a generator hits its reactive power 
limit or control mode changes in the DC link 

(3) Considering the static model, it has been found that if the load is varied at only one bus 
in an integrated AC-DC system, the collapse of voltage takes at a bus where the load is 
being increased However, with simultaneous changes m the loadings at all the buses, all 
the bus voltages approach the nose point of the P-V curves simultaneously 

(4) Consideration of reactive power output limits of synchronous machines reduces the loading 
at which the system becomes voltage unstable Thus, the maximum system loadability is 
reduced, whereas it is increased by considering the loads to be voltage dependent 

Chapter 3 compares the effect of four different generation rescheduling schemes in increasing the 
voltage stability margin in reactive parameter space The four schemes of generation reschedul- 
ing has been formulated in view of minimising total fuel cost of generation (OBJ-I), minimising 
total system transmission loss (OBJ-II), maximising minimum singular value of the load flow Ja- 
cobian (OBJ-III) and minimisation of the slack reactive power injection (new scheme OBJ-IV) 
The main findings are as follows 

(1) The generation rescheduling using any of the four objectives increases the voltage stability 
margin m stressed condition of system The stability margin can be measured in terms of 
Q-margin or the minimum singular value of the power flow Jacobian 

(2) At light loading conditions, none of the schemes is able to increase the Q-margm The 
minimum singular value, however, is slightly improved with OBJ-III and IV In such cases, 
the correlation between minimum singular value of the power flow Jacobian and Q-margin 
is lost 
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(3) Amongst the four different objectives tried out, the new objective (OBJ-IV) hats consis- 
tently resulted into better results than other objectives The results obtained are very 
close to the OBJ-III 

(4) The new objective of minimising the slack bus reactive injection is computationally less 
expensive and can be effectively used for voltage stability margin enhancement 

Chapter 4 has investigated certain aspects of local bifurcations and its implications on power 
system voltage stability Considering the detailed model of generators along with AVR, speed- 
governor control loop and induction motor loads Investigations reveal that* 

(1) Power system models may exhibit static and dynamic bifurcations even for parameter 
values in their normal operating range. 

(2) For the same value of bifurcation parameter, the system may become unstable either 
purely due to angle instability or purely due to the voltage instability or both may occur 
simultaneously, depending on the disturbance in the system. 

(3) A system becomes unstable due to Hopf bifurcation and exhibits oscillatory behaviour for 
small perturbations It has been shown that for certain choices of induction motor load 
torque coefficients or AVR gam values, the power system may undergo oscillatory type 
instability that is essentially bifurcation of a fixed point into a periodic orbit 

(4) Hopf bifurcation may occur both at low or high power factor loadings The effect of Hopf 
bifurcation is to further reduce the stable operating range of parameters 

Chapter 5 has investigated various types of global bifurcations and its implications on voltage 
stability The main findings are as follows 

(1) Different types of dynamic bifurcations such as PDB, CFB, TRB, LP and chaos can exist 
in the power system networks for variation of some of the parameters even in their normal 
range 

(2) Some of the global bifurcations observed, possibly for the first time, in the power system 
networks include the bifurcation of periodic orbit to torus (TRB) and limit points of 
unstable to unstable periodic orbits 
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(3) Chaos has been observed m all the three sample systems, in two of the systems via PDB 
route. In one of the systems, quasipenodic route to chaos has been observed for the 
first time in power system model The chaotic oscillations produce harmonics of all the 
frequencies and a continuous broad band frequency spectrum is observed 

(4) Boundary crisis (blue sky bifurcation) and subcritical Hopf bifurcation cause voltage col- 
lapse in some of the systems By adjusting the SVS gam parameter, the boundary crisis 
and hence the resulting voltage collapse can be avoided 

Chapter 6 has dealt with elimination of dynamic bifurcation and chaos using FACTS devices 
such as controllable series capacitors (CSC), static phase angle regulator (PAR), static VAR 
compensators (SVC) The important conclusions drawn from the studies presented m this chap- 
ter are as follows 

(1) On varying the controller gain, the subcritical Hopf bifurcation becomes supercritical 
after a certain gain value with cyclic fold bifurcation disappearing Further variation of 
controller gain value causes other dynamic bifurcations to disappear yielding stable steady 
state operating point, m the earlier Hopf bifurcation window 

(2) The choice of appropriate feedback control signal is very important In one of the systems 
studied, PAR could not eliminate the dynamic bifurcation and chaos with generator ac- 
celeration power ( P a ) as the error signal and the SVC could not eliminate the bifurcation 
with generator speed deviation as the error signal However, they could eliminate the HB 
window when the choice of input signals was reversed Similarly in another system, CSC 
did not work with P a as error signal 

(3) The CSC with generator speed deviation as error signal has been able to eliminate dynamic 
bifurcation and chaos in all the sample systems studied 

(4) The FACTS devices such as CSC, PAR, SVC can be effectively used for controlling the 
occurrence of dynamic bifurcations and chaos by proper choice of the error signal and gain 
of the controller 


7.3 Scope for Further Research 

As a consequence of investigations carried out in this thesis, the following aspects aie being 
suggested as future research work to be carried out 
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(1) The static voltage stability margin in the present work and also by other researchers has 
been computed with respect to saddle node bifurcation m multidimensional parameter 
space. The studies may be extended to estimate the stability margin with respect to the 
Hopf bifurcation m multidimensional parameter space. 

(2) More rigorous steady state model of synchronous machine may be considered to include 
the effect of machine saturation and under frequency operation, on the heating of the 
machine This will affect the static voltage stability margin 

(3) In Chapter 4, it was observed that the dominant states at Hopf bifurcation considering 
excitation system model are the E' q and Ejj. Studies to eliminate dynamic bifurcation can 
be extended for such cases by utilising these variables variables as feedback control signals 
and incorporating washout filters to ensure preservation of steady state equilibrium The 
effect of power system stabilisers on control of dynamic bifurcation can also be studied 

(4) In the present work, the location of FACTS devices were selected heunstically in the 
sample systems for control of dynamic bifurcations However, a more systematic approach 
need to be evolved for optimal siting of these devices 

(5) Due to extensive computing requirements to study the global bifurcations and chaos, the 
research activities in this thesis and also in other works are confined to smaller systems only 
It will be worthwhile developing techniques and models for study of nonlinear dynamics of 
large size systems This may require exploring new methods to obtain network equivalents 
suitable for the voltage stability analysis 
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Appendix A 


Data for 2 Bus Test System (at 100 
MVA base) 

This system is shown m Fig A 1 The relevant data are provided in the following tables 


Table A 1 Generator Data 



MBiTirireoi 


H = 7 76 sec 

t— 

oo 

CD 

o 

II 

* 

£3SH2E£1 

mEQQZDi 

T qo = 0 02 sec 

R a = 0 000909 



Table A 2 Exciter Data 


K a = 10 0 

T a = 0 06 sec 

gWHilifehfcl 

E39^3QE2S3S 

Tj = 1 00 sec 

K e = -0 08 

T e = 0 405 sec 



v; mm = -3 5 


Table A 3 Governor Data 


T x = 124 47 

T 2 = 8 590 

Unj^igSEl 

^ J 

II 

o 

o 

o 

o 

T s = 0 740 


F = -2 000 

P ma * = 2 670 
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Table A 4 Induction Motor Data 



Table A 5 Transmission Line Data 



Figure A 1 ‘2-bus system 





Appendix B 

DATA FOR 3 BUS TEST SYSTEM 

The 3 bus system is shown in Fig B 1 The relevant data for this study system is given below 


Table B 1 Generator Data 


Pm 

19 

D 

Voltage of Inf bus-1 and Gen bus-3 = 10 


t 

| 



0 05 

1 


Table B 2 Load Data 


Pi =00 

Qx = 0 0 

* 

H 

O 

o> 

K pw = 0 4 

'-O 

o 

II 

ii 

00 

Ol 

Qo = 1 3 

K qw = -0 03 

K qv =-2 8 







Table B 3 Line Data 


Line no 

From bus 

BiCTIBIM 

Admittance 

UBUS/ 2 

i 

1 

m 

20 Z— 85 deg 

0 

2* 

2 

3 

10 Z — 85 deg 



* Bus 2 has compensation of 3 0 p u 
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Fig 


3-bus system 



Appendix C 


DATA FOR 3 BUS TEST SYSTEM 
WITH SVS 

The 3 bus system is shown in Fig Cl The relevant data for this study system is given below 


Table C 1 Generator Data 


Pm 

M g 

D 

Voltage of infinite bus-1 and generator bus-3 = 10 





1 0 

05 

0 05 ] 



Table C 2 Load Data 


Pi = 0 0 

© 

II 

o 

o 

II 

<s 

Qv = 10 

CD 

II 

e = ie 

P o = 0 3 

K pu = 0 6 

o 

II 

p 

*! 

T = 85 

Q o = 13 

K q „ = -0 05 

K qv = -0 8 

I< qV2 =23 






Table C 3 SVS Data 


o 

r — 4 

ii 

K sv , = 50 0 

T 3va = 0 01 

B mtn = -15 0 

Rnax = 15 0 


Table C 4 Line Data 


Line no 

From bus 

To bus 

Admittance 

VBUS / 2 

1 

i 

2 

202-82 deg 

0 

2 

2 

3 

102—85 deg 

0 

3 

3 

1 

2 82—82 deg 

0 

★ 

Bus 2 has compensal 

non of 10 0 p u 
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Appendix D 


DATA FOR 4 BUS TEST SYSTEM 

The 4-bus system is shown in Fig D 1 The relevant data for this study sj stem is given below 


Table D T Generator Data 


Pm =10 


M a = 0 5 


D = 0 05 


Voltage of infinite bus-1 and generator bus-3 = 10 


Table D 2 Load Data 


Mil 


<0 

II 

O 

o 

Qv = 1 0 

T] =2 1 

^ = 14 

Mil 

K puj = 05 

K pv = 0 3 


Qo = lZ 

K qw = -0 03 

K qv = -0 8 

K qVi =23 






Table D 3 OLTC Data 


An = 0 025 

At = 10 sec 

e = 0015 

n mt „ = 09 

^mar 11 

O 

f-H 

II 

O 


Table D 4 Line Data 


Line no 

From bus 

To bus 

Admittance 

Y sh , 2 

i 

1 

4 

20Z-82 deg 

0 

2 

4 

3 

10 Z — 85 deg 

0 

3 

3 

1 

2 8 Z— 80 deg 

0 

4* 

4 

2 

100 Z0 deg 

0 

*1 

3us-4 has compensation of 12 0 p u 
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Figure D 1 4-bus system with OLTC 



Appendix E 

Data for 9 Bus Test System (at 100 
MVA base) 

The 9 bus system is shown in Fig E The system data is taken from ref [7] The relevant data 
are provided in following tables 


Table E 1. Line Data 


From 

bus 

To 

bus 

Series impedence 

Shunt susceptance 
B(pu) 

R(p u ) 

X(p u ) 

1 

4 

0 0000 

0 0576 

0 0000 

4 

6 

0 0170 

0 0920 

01580 

6 

9 

0 0390 

0 1700 

0 3580 

9 

3 

0 0000 

0 0586 

0 0000 

9 

8 

00119 

0 1008 

0 2090 

8 

7 

0 0085 

0 0720 

0 1490 

7 

2 

0 0000 

0 0625 

0 0000 

7 

5 

0 0320 

0 1610 

0 3060 

5 

4 

0 0100 

0 0850 

01760 
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Table E 2 Machine Data 


Parameter 

Machine- 1 

Machme-2 

Machine-3 

^(p U ) 

0 1460 

0 8958 

1 3125 

*5(p u ) 

0 0608 

0 1198 

0 1813 

^ 0 (sec ) 

8 96 

6 00 

5 89 

X q (p u ) 

0 0969 

HBUgg] 

1 2578 

X'(p « ) 

0 0608 

urn 

0 1813 

wmam 

0 310 

0 535 

0 600 


23 64 

6 40 

3 Of 

| D(p u ) 

1 0 

1 0 

1 0 


Table E 3 Bus Data 


Bus 

No 

Bus 

Type 

Generation(p u ) | 

Load(p u ) 

Voltage 

Magnitude 

Real 

Reactive 

Real 

Reactive 

1 

Swing 

— 

— 

0 00 

■KH9 

1 0400 

2 

PV 

1 63 

— 

0 00 

wami 

1 0253 


PV 

0 85 


0 00 


1 0253 


PQ 

0 00 

0 00 

pHi| 

0 00 

— 

mm 

PQ 

WM 

0 00 

KE23 

0 50 

— 

mi 



0 00 

nn 

0 30 

— 



BiltTil 

0 00 


0 00 

— 

8 



mmm 

WM 

■■ill 

■■■■BH 

ma 



wmm 

mm 

IfPKIll 

— 


























































APPENDIX E DATA FOR 9 BUS TEST SYSTEM (AT 100 MVA BASE) 


Load C 



Appendix F 


Data for 10 Bus Test System (at 100 
MVA base) 

This system is derived from 9-bus system given in Appendix-E by connecting an induction motor 
load in addition to the existing load to bus 6 (Fig F.l) The data for excitation system and 
speed governing loop were taken from ref [7] in accordance with the capacity of the machines 
The relevant data are provided m following tables 


Table F.l Excitation System Data 



Exciter-1 

Exciter-2 

Exciter-3 

K a 

10 00 

10 00 

10 00 

K f 

0 0648 

0 0910 

0 1080 

K e 

-0 0200 

-0 0505 

-0 0601 



0 20 

0 20 

EBB 

HBiEfli 

0 35 

0 35 

T e (sec ) 

0 1000 



A ex 

MiMM 

■iIiliM 


Bex 

Iflif 

1 3733 

■liMM 

pm 

■EMM 

2 00 

2 00 

msm 


-2 00 

-2 00 
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Table F 2 Speed Governor Data 


Parameter 

Governor- 1 

Governor-2 


Ti( sec ) 

30 0 

0 083 

0 083 

T 2 ( sec ) 

35 

0 00 

0 000 

T 3 (sec ) 

0 520 

0.200 

0 200 

T 4 ( sec ) 

00 

0 050 

0 050 

T 5 (sec ) 

0 415 

8 000 

5 000 

R'g ( P u ) 

0 050 

0.050 

0 050 

F 

-2 000 

0.271 

0 280 

P 

1 max 

2 50 

1 75 

1 32 


Table F 3 Induction Motor Data 



R 1 

H 

X 2 

Ri 


3 5065 

0 03057 

1 0 

3 5065 1 

0 03596 ! 

3 399 





























Appendix G 


Data for IEEE-14 Bus Test System (at 
100 MVA base) 


The IEEE-14 bus system is shown in Fig G 1 The system data is taken from ref [39] and buses 
renumbered The relevant data are provided m following tables 


Table G 1 Generator Bus Data 


Bus 

No 

Scheduled real 
power generation 

P G { MW) 

Specified Voltage 
magnitude 
Vspec (p u ) 

Load 

Real 

(MW) 

Reactive 

(MVAR) 

1 (slack) 

— 

1 060 

00 00 

00 00 

2 

040 0 

1 045 

21 70 

IBSEI 

3 

020 0 


nail 


4 

— 


P’TfTi 


5 

— 


00 00 

HEEl 


Table G 2 Generator Data 


Gen 

no 

Real generation 
limit 

React generation 
limit 

Cost characteristics* 


Maximum 

(MW) 

Minimum 

(MW) 

Maximum 

(MVAR) 

Minimum 

(MVAR) 

a, 

($/MW 2 -hr) 

6, 

($/MW-hr) 

c, 

(S/hr) 


200 0 


100 0 

-45 0 

1 0 

2 45 

105 00 

EHI 

100 0 


050 0 

-40 0 

1 0 

3 51 

044 40 

EH 

100 0 


024 0 

-06 0 

1 0 

3 89 

040 60 

E 

— 

— 

040 0 

000 0 

— 

— 

— 


— 

— 

024 0 

-06 0 

— 

— 

— 


* Fuel cost of i th gen unit F t = (\a t P Gt + b t P Gi + c,)$/hr 
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Table G 3 Load Bus Data 


Bus 

No 

Load 

Shunt compensation 
m(MAVR) 

External shunt 
susceptance(p u ) 

Real(MW) 

Reactive(MVAR) 

6 

00 0 

00 0 

000 0 

0 00 

7 

29 5 

16 6 

-19 0 

0 19 

8 

07 6 

01 6 

000 0 

0 00 

9 

47 8 

-3 9 

000 0 

0 00 

10 

09 0 

05 8 

000 0 

0 00 

11 

03 5 

01 8 

000 0 

000 

12 

06 1 

01 6 

oooo 

0 00 j 

13 

13 5 

05 8 

000 0 

0 00 

14 

14 9 

05 0 

oooo 

000 1 


Table G 4 Transformer Data 


Transf no 

From 

To 

Series reactance 

Tap settings 

(Line no ) 

bus 

bus 

X(p.u) 

(p u )(base case) 

i 

8 

3 

0 25202 

0 962 

2 

9 

6 

0 20912 

0 978 

3 

9 

7 

0 55618 

0 969 


Table G 5 Line Data 


Line 

From 

To 

Senes impedence j 

Shunt susceptance 

no 

bus 

bus 

R(pu) 

X( P u) 

B(p u ) 

4 

1 

8 

0 05403 

0 22304 

0 0246 

5 

2 

8 

0 05695 

0 17388 

0 0170 

6 

4 

9 

0 06701 

6 17103 

0 0173 

7 

9 

8 

0 01335 

0 04211 

0 0064 

8 

1 

2 

0 01938 

0 05917 

0 0264 

9 

2 

4 

0 04699 

0 19797 

0 0219 

10 

6 

5 

0 00000 

0 17615 

0 0000 

11 

2 

9 

0 05811 

0 17632 

0 0187 

12 

6 

7 1 

0 00000 

0 11001 

0 0000 

13 

7 

10 

0 03181 

0 08450 

0 0000 

14 

3 

11 

0 09498 

0 19890 

0 0000 

15 

3 

12 

0 12291 

0 25581 

0 0000 

16 

3 

13 

0 06615 

0 13027 

0 oooo 

17 

7 

14 

0 12711 

0 27038 

0 0000 

18 

10 

11 

0 08205 

0 19207 

0 0000 

18 

12 

13 

0 22092 

0 19988 

0 0000 

20 

13 

. 14 

0 17093 

0 34802 

0 0000 









Appendix H 


Data for 19 Bus Test System (at 100 
MVA base) 


The UPSEB 19-bus system is shown in Fig H 1 The system data is taken from UP STATE 
ELECTRICITY BOARD STUDY CELL, IIT KANPUR and buses renumbered The relevant 
data are provided in following tables 


Table H 1 Generator Bus Data 


Bus 

No 

Scheduled real 
power generation 
Rg(MW) 

Specified Voltage 
magnitude 
Vspec (p u ) 

Load 

React Generation 

Real 

(MW) 

Reactive 

(MVAR) 

Maximum 

(MVAR) 

Minimum 

(MVAR) 

1 (slack) 

— 

1 010 

0 00 

0 00 

— 

— 

2 

1600 0 

1 058 

0 00 

0 00 

960 0 

-480 0 

3 

0900 0 

1 056 

0 00 

0 00 

540 0 

-050 0 

4 

0550 0 

1 057 

0 00 

0 00 

330 0 

-050 0 


Table H 2 Transformer Data 


Transf no 

From 

To 

Series resistance 

Series reactance 

Tap settings 

(Line no ) 

bus 

bus 

R(p u ) 

X(p u) 

(p u )(base case) 

1 

5 

2 

00016 

01991 

1000 

2 

6 

1 

.00073 

01460" 

1 000 

3 

7 

3 

00030 

03499 

1000 

4 

8 

4 

00049 

04743 

1 000 
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Table H 3 Load Bus Data 


Bus 

No 

Load 

External shunt 
susceptance(p u ) 

Real(MW) 

Reactive( M VAR) 

5 

00 0 

00 0 

-2 02269 

6 

191 6 

16 0 

-0 907029 

7 

1000 0 

00 0 

-0 571429 

n 

00 0 

00 0 

-0 453515 

9 

135 0 

78 0 

-0 453515 

10 

236.0 

45 0 

-0 453515 

11 

360 0 

11 0 

-1 02495 

12 

337 0 

59.0 

-0 453515 

13 

90 0 


00 

14 

520 0 

bhhhesei 

-0 453515 

15 

387 0 

64 0 

-0 907029 

16 

288 0 

00 

-2 6304 

17 


0 0 

-0 453515 

18 


69 0 

00 

19 


00 

00 1 


Table H 4 Line Data 


Line 


Bl 

Senes impedence 

Shunt susceptance 

no 

mi 

iS9i 

R(pu) 1 

X(p u ) | 

B(p u ) 

5 

5 

Bi 

iiMiSfl 


04056 

6 

7 

16 

■j 


1 21680 

7 

16 

Bl 

00927 

09429 

1 23293 

8 

16 

Bl 

00833 

08478 

1 10855 

wm 

Bl 

B! 

00031 

00310 

■■EisiifliBi 


6 

m 

■ran™ 

00517 

■hbi 


8 


KB 

04880 

63614 


10 

9 

00254 

02584 

33798 

13 

6 

10 

00468 

04770 


14 

10 

11 

00294 

02997 

ihhiI 

15 

mm 

5 

00823 

08375 

1 09503 

16 

mm 

12 

■m 

06617 

86521 

17 

mm 

12 


— H 1 

||BKlM9Bi| 

18 

13 

14 

ilTOM 

— 

\mmmmam] 

19 

14 

12 

1— i 

02646 

imihhsbhI 

20 

15 

14 

IBBl 

08169 

1 06808 

21 

6 

15 

i— a 

07990 

1 04738 

22 

15 

16 

00015 

HU 

02704 

23 

16 

17 

00559 

mmm 

74354 

24 

17 

IB 

00121 

01109 

72815 

25 

14 

HB 

00051 

00517 

06760 
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Table H 5 Machine Data 


Parameter 

Generator # 2 

Generator #3 

Generator #4 

M. 

0 266 

0 178 

0 084 

X'i 

0 014 

0 023 

0 028" 

D 

1 0" 

1 0 

1 0 


LKOW SLNT 



Figure H 1 19-bus UPSEB system 



















Appendix I 


Data for IEEE- 30 Bus Test System (at 
100 MVA base) 

The IEEE-30 bus system is shown m Fig II The system data is taken from ref [39] and buses 
renumbered The relevant data are provided m following tables 


Table I 1 Generator Bus Data 


Bus 

No 

Scheduled real 
power generation 
P g (MW) 

Specified Voltage 
magnitude 
Vspec (p u ) 

Load 

Real 

(MW) 

Reactive 

(MVAR) 

1 (slack) 

— 

1 060 

00 00 

00 00 

2 

040 0 

1 045 

21 70 

12 70 

3 

020 0 

1 010 

00 00 

■Mif 

4 

— 

1 082 

30 00 

BiEB 

5 

— 

1 010 

09 42 

■BDEH 

6 

— 

1 071 

00 00 

00 00 


Table I 2 Transformer Data 


Transf no 
(Line no ) 

From 

bus 

To 

bus 

Senes reactance 
X(p u) 

Tap settings 
(p u )(base case) 

1 

13 

|07 


0 978 

2 

mm 


HKESSHH 

0 969 

3 

mm 


n 

0 962 

4 

mm 


0 3960 

0 968 
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Table I 3 Generator Data 



Real generation 
limit 


Maximum Minimum 


(MW) 


300 0 


150 0 


150 0 


(MW) 


050 0 


020 0 


020 0 



React generation 
limit 

Maximum 

(MVAR) 

Minimum 

(MVAR) 

150 0 j 

-100 0 

050 0 1 

-040 0 

040 0 

- 010 " 0 

024 0 

-006 0 

040 0 

-040 6 

024 0 

-006 0 


Cost characteristics 



d t 

($/MW 2 -hr) 


1 0 


1 0 


1 0 


* Fuel cost of i th gen unit 

Table I 4 Load Bus Data 
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Table I 5 Line Data 

























































































































(a) ONE LINE DIAGRAM 

(b) 3 - WINDING TRANSFORMERS 
EQUIVALENTS 



SHUNT REACTOR / 
CAPACITOR 

GENERATOR 


SYNCHRONOUS 

CONDENSER 

TRANSFORMER 


Figure I 1 IEEE-30 bus system 




Appendix J 


Data for UPSEB-156 Bus Test System 
(at 100 MVA base) 


The UPSEB 156-bus system is shown in Fig J.l The system data is taken from UP STATE 
ELECTRICITY BOARD STUDY CELL, IIT KANPUR and buses renumbered The rele\ant 
data are provided m following tables 


Table J 1 Generator Bus Data 


Bus 

Scheduled real 

Specified Voltage 

Load 

React. Generation 

No 

power generation 

magnitude 

Real 

Reactive 

Maximum 

Minimum 


P a { MW) 

Vspec (p.u ) 

(MW) 

(MVAR) 

(MVAR) 

(MVAR) 

1 (slack) 

— 

1030 

00 00 

00 00 

— 

— 

2 

030 0 

wmBSEEMM 

00 00 

00 00 

20.00 

00 00 

3 

360 0 

■j 

00 00 

00 00 


00 00 

4 

180 0 

1 050 

KnilinM 

00 00 

mzMm 

■H 

5 

099 0 

1030 

Buna’ 

00 00 

44 00 

1 

6 

100 0 

wKmERMM 

00 00 

lifiliM 

60 00 

■ 

7 

250 0 

BHEES9HH 

00 00 

wSmSM 

■Hf 

00 00 

8 

230 0 

HHES3HHI 

■IIIIIIIB 

00 00 


00 00 

9 

140 0 

1 050 

iE&ESI 



00 00 

10 

180 0 

1 030 

00 00 

— 

Bia 

00 00 

11 

060 0 

ihheeshhi 

Binlritl 


19 00 

iHKB 

12 

180 0 

II 



mmm 

— 

13 


II 

IKItgifni 

WMM 


mmm 

14 


1030 

n 

mmm 

50 00 

mmm 

15 

020 0 

1 030 

1— 

IHH 

10 00 

00 00 


Contd 
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Table J.l (Contd ) . . 


Bus 

No. 

Scheduled real 
power generation 

Po( MW) 

Specified Voltage 
magnitude 
Vspec (p u ) 

Load 

React Generation 

Real 

(MW) 

Reactive 

(MVAR) 

Maximum 

(MVAR) 

Minimum 

(MVAR) 

16 

018 0 

1020 

00 00 

00 00 

04 00 

00 00 

17 

032 0 

1 020 

00 00 

00 00 

08 00 

00 00 

18 

020 0 

1 020 

00 00 

00 00 

09 00 

00 00 

19 

1600 0 

1 030 

00 00 

00 00 

144 00 

00 00 

20 

900 0 

1 030 

00 00 

00 00 

80 00 

00 00 

21 

150 0 

1 030 

00 00 

00 00 

84 00 

00 00 

22 

454 0 

1 000 

00 00 

00 00 

35.00 

00 00 

23 

130 0 

1 020 

00 00 

00 00 

56 00 

00 00 

24 

040 0 

1.000 

56 00 

09 00 

31 00 

00 00 

25 

147 0 

1 000 

00 00 

00 00 

00 00 

00 00 ] 

26 

1000 0 

1 000 

00 00 

00 00 

00 00 

00 00 


Table J 2 Load Bus Data 


Bus 

Load 

Shunt compensation 

External shunt 

No 

Real (MW) 

Reactive(MVAR) 

m(MAVR) 

susceptance(p u ) 

26 

735 0 

00 0 

000 0 

0 00 

27 

00 0 

00 0 

000 0 

0 00 

28 

05 8 

02 0 

00 0 

0 00 

29 

11 2 

07 5 

000 0 

0 00 

30 

00 0 

-0 0 

-19 0 

019 

31 

07 6 

01 6 

000 0 

0 00 

32 

22 8 

10 9 

000 0 

0 00 

33 

00 0 

00 0 

000 0 

0 00 

34 

06 2 

01 6 

000 0 

0 00 

35 

08 2 

02 5 

000 0 

0 00 

36 

03 5 

01 8 

000 0 

0 00 

37 

09 0 

05 8 

000 0 

0 00 

38 

03 2 

00 9 

000 0 

0 00 

39 

09 5 

03 4 

000 0 

0 00 

40 

02 2 

00 7 

000 0 

0 00 

41 

17 5 

11 2 

000 0 

0 00 

42 

00 0 

00 0 

000 0 

0 00 

43 

03 2 

01 6 

000 0 

0 00 

44 

08 7 

06 7 

-04 3 

0 043 

45 

o"o6 

00 0 

000 0 

0 00 

46 

03 5 

02 3 

000 0 

0 00 


Contd 
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Table J.2 (Contd ) . . 


Bus 

Load 

Shunt compensation 

External shunt 

No 

Real(MW) 

React i ve ( M VA R) 

m(MAVR) 

susceptance(p u ) 

47 

02 4 

01.2 

000 0 

0 00 

48 

00 0 

00.0 

000 .0 

0 00 

49 

02 4 

00 9 

000 0 

0 00 

50 

■HI 

01 9 

000 0 

0 00 

51 

17 5 

11 2 

000 0 

0 00 

52 

00 0 

00 0 

000 0 

0 00 

53 

03 2 

01 6 

000 0 

0 00 

54 

08 7 

06 7 

-04 3 

0 043 

55 

HEEQHI 

U 

000 0 

000 

■a 

■KSSH 


000 0 

0 00 

mi 

02 4 

01 2 

000 0 

0 00 

58 

00 0 

mh&hhii 

000 0 


59 

BHI 

■I 

000 0 

0 00 

60 

1 

01.9 

000 0 

0 00 

si 

07 6 

01.6 


0 00 

m 

22 8 

10 9 


0 00 

63 

00.0 

00.0 

000 0 

0 00 

64 

06 2 

■BDH1 

000 0 

0 00 

65 

08 2 


hhesssbhii 

000 

66 

03 5 

01 8 


0 00 

67 

09 0 

05.8 

■ 

0 00 

68 

03 2 

00.9 

000 0 

0 00 

69 

09 5 

03.4 

■■■EOiEHiHi 

0 00 

70 

02 2 

00 7 


0 00 

71 

17 5 

11 2 


0 00 

72 

00 0 

00 0 


0 00 

73 

03 2 

ihhbohi 

ihhsuzihhb 

0 00 

74 

mmm 


-04 3 

0 043 

nm 


00 0 

000 0 

0 00 

lm 

IHUH 

02 3 

■■■■■ 

0 00 

77 

02 4 

01 2 

IHHH 

0 00 

78 

1—TilM 

00 0 

000 0 


79 

mmam 

00 9 

IHHH 

IHH 

80 

10 6 

01 9 

jlHBOlMHHi 


81 

17 5 

11 2 

IH0H 


82 



iwmmmm 

0 00 

83 

IHHI 

IHHI 

oooo 

IHHI 


Contd 
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Table J 2 (Contd ) . . 


Bus 

Load 

Shunt compensation 

External shunt 

No. “ 

Real(MW) 

Reactive(MVAR) 

in(MAVR) 

susceptance(p.u ) 

84 

08 7 

06 7 

-04.3 

0.043 

85 

00 0 

00 0 

000 0 

0 00 

86 

03 5 

023 

000 0 

0 00 

87 

02 4 

01 2 

000 0 

0 00 

88 

00 0 

■HI 

000 0 

0.00 

89 

02 4 

00 9 

000 0 

0 00 

90 

10 6 

01 9 

000 0 

0 00 

91 

07 6 

01 6 

000 0 

0 00 

92 

22 8 

10 9 

000 0 

0 00 

93 

PBliflH 

00 0 

000 0 

0 00 

94 

m 

■HE99HM 

000 0 

0 00 

95 

08 2 

n 

000 0 

0.00 

96 

03 5 

01 8 

000 0 

0 00 

97 

09 0 

05 8 

000 0 

0.00 

98 

03 2 

00 9 

000 0 

0 00 

99 

09 5 



0 00 

100 


■j 

000 0 

0 00 

101 

17 5 

11 2 

000 0 

0 00 

|K m 

00 0 


000 0 

0.00 

In 

03 2 


000 0 

0 00 

104 

iHBSBHI 

hhbsqhh 

-04 3 

0 043 

105 



000 0 

0 00 

106 

03 5 

■mkiiohmi 

000 0 

0.00 

107 

02 4 


000 0 

0 00 

108 

■KnOH 


■■SHI 

0 00 

IllSltsI 


wmmmam 


0 00 

lug 

10 6 

01 9 


0 00 

111 

07 6 

01 6 

000 0 

0 00 

112 

22 8 

(■■Hi 

000 0 

0 00 

113 

1021 

IBHBH 

000 0 

0 00 

114 

06 2 

01 6 

000 0 

0 00 

115 

08 2 

02 5 

000 0 

0 00 

116 

iwmmm 

01 8 

000 0 

0 00 

117 

!■ 

05 8 

000 0 

0 00 

118 

032 

00" 9 

ihhkhhb 

0 00 

119 

09 5 

03 4 


0 00 

120 

02 2 

00 7 

000 0 

0 00 


Contd 
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Table J 2 (Contd ) 


Bus 

No 

Load 

Shunt compensation 
m(MAVR) 

External shunt 
susceptance(p.u ) 

Real(MW) 

Reactive(MVAR) 

121 

17 5 

11 2 

000 0 

0.00 

122 

00 0 

00 0 

000.0 

0.00 

123 

03 2 

01 6 

000 0 

0 00 

124 

08 7 

06 7 

-04.3 

0.043 

125 

00 0 

00 0 

000.0 

000 

126 

03 5 

02 3 

000 0 

0 00 

127 

02 4 

01 2 

000 0 

0 00 

128 

00 0 

00 0 

000 0 

0 00 

129 

02 4 

00 9 

000 0 

0 00 

130 

10 6 

01 9 

000 0 

0 00 

131 

07 6 

01 6 

000 0 

0 00 

132 

22 8 

10 9 

000 0 

0 00 

133 

00 0 

00 0 

000 0 

0 00 

134 

06 2 

01 6 

000 0 

0.00 

135 

08 2 

02 5 

000 0 

0 00 

136 

03 5 

01 8 

000 0 

0 00 

137 

09 0 

05 8 

000 0 

0.00 

138 

03 2 

00 9 

000 0 

0.00 

139 

09 5 

03 4 

000 0 

0.00 

140 

02 2 

00 7 

000 0 

0 00 

141 

17 5 

11 2 

000 0 

0.00 

142 

00 0 

00 0 

000 0 

0 00 

143 

03 2 

01 6 

000 0 

0 00 

144 

08 7 

06 7 

-04 3 

0 043 

145 

00 0 

00 0 

000 0 

0 00 

146 

03 5 

02 3 

000 0 

0.00 

147 

02 4 

01 2 

000 0 

0 00 

148 

00 0 

00 0 

000 0 

000 

149 

02 4 

00 9 

000 0 

0.00 

150 

10 6 

01 9 

000 0 

0 00 

151 

07 6 

01 6 

000 0 

0 00 

152 

22 8 

10 9 

000 0 

0 00 

153 

00 0 

00 0 

000 0 

0 00 

154 

06 2 

01 6 

000 0 

0 00 

155 

08 2 

02 5 

000 0 

0 00 

156 

03 5 

01 8 

000 0 

0 00 
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Table J 3 Transformer Data 


Transf no 

From 

To 

Series resistance 

Series reactance 

Tap settings 

(Line no ) 

bus 

bus 

R(pu) 

X(p.u ) 

(p.u )(base case) 

1 

42 

54 

.00170 

04330 

1 000 

2 

42 

57 

.03540 

18180 

1 000 

3 

42 

3 

00123 

02469 

1000 

4 

93 j 

16 

01666 

33333 

1 000 1 

5 

86 

10 

00000 

02917 

1 000 1 

6 

50 

42 

00065 

02604 

1 000 | 

7 

50 

1 

00073 

01460 

1.000 

8 

90 

88 

00133 

04660 

1.000 

9 

29 

71 

00143 

05000 

1.000 

10 

30 

29 

00065 

02604 

1 000 ® 

11 

28 

27 

00110 

03770 

1.000 

12 

32 

31 

00095 

03330 

1.000 

13 

97 

46 

00065 

02604 

1.000 

14 

33 

32 

00065 

.02600 

1.000 

15 

35 

34 

00090 

03200 

1 000 

16 

36 

35 

00065 

02604 

1 000 

17 

38 

37 

00286 

.10000 

1 000 

18 

40 

39 

00143 

05000 

1 000 

19 

46 

45 

00153 

.05250 

1.000 

20 

56 

6 

00306 

.06135 

1 000 1 

21 

56 

55 

00136 

04760 

1 000 

22 

60 

59 

00143 

05000 

1 000 

23 

63 

62 

00143 

05000 

1.000 

24 

64 

63 

00043 

01736 

1 000 

25 

66 

65 

00133 

04644 

1 000 

26 

69 

68 

00143 

05000 

1000 

27 

73 

72 

00150 

05160 

1 000 

28 

76 

9 

00344 

06875 

1 000 

29 

78 

77 

00150 

05160 

1 000 

30 

80 

79 

00148 

05200 

1 000 

31 

82 

81 

00143 

05000 

1 000 

32 

89 

12 

00235 

04710 

1000 

33 

95 

17 

01040 

20800 

1 000 

34 

85 

18 

01766 

35313 

1 000 

35 

99 

2 

00000 

19034 

1 000 

36 

91 

13 

00514 

10285 

1000 

37 

87 

11 

00549 

10978 

1 000 1 


Contd 
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Table J 3 (Contd ) 


?ransf no 

From 

To 

Series resistance 

Series reactance 

Tap settings 

(Line no ) 

bus 

bus 

R(pu) 

X(p u.) 

(p.u )(base case) 

38 

51 

4 

00000 

.04860 

1.000 

39 

92 

14 

00049 

01943 

1000 

40 

112 

108 

00065 

02604 

1.000 

41 

74 

73 

00130 

05208 

1 000 

42 

94 

96 

00286 

10000 

1.000 

43 

58 

44 

00095 

03333 

1 000 

44 

57 

7 

00153 

03600 

1 000 

45 

54 

5 

00556 

11110 

1.000 

46 

31 

15 

01250 

25000 

1 000 j 

47 

32 

23 

00243 

04860 

1.000 1 

48 

108 

107 

00143 

05000 

1 000 i 

49 

129 

105 

00143 

05000 

1 000 1 

50 

126 

104 

00286 

10000 

1 000 1 

51 

124 

8 

00770 

02720 

i ooo j 

52 

103 

102 

00286 

10000 

14)00 

53 

109 

19 

00016 

00591 

1 000 | 

54 

113 

20 

00030 

01199 

1000 

55 

123 

j 41 

00286 

10000 

1000 

56 

61 

21 

00000 

02841 

1 000 

57 

127 

22 

00000 

02273 

1000 

58 

128 

127 

00056 

02222 

0 950 

59 

155 

140 

00286 

10000 

1 000 

60 

26 

154 

00147 

05830 

1000 

61 

133 

153 

00286 

10000 

1 000 

62 

51 

152 

00286 

10000 

1000 

63 

124 

131 

00286 

10000 

1 000 1 

64 

67 

156 

00143 

05000 

1000 

65 

118 

130 

00150 

05160 

1 000 I 

66 

119 

125 

00300 

10320 

1 000 

67 

43 

122 

00150 

05160 

1000 

68 

84 

121 

00286” 

10000 

1 000 1 

69 

52 

151 

00160 

05600 

i ooo i 

70 

116 

114 

00286 

10000 

1 ooo 

71 

111 

53 

00150 

05160 

1 ooo 1 

72 

120 

47 

00300 

10320 

1 ooo 

J 
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Table J 4: Line Data 


Line 

From 

To 

Series impedence 

Shunt susceptance 

no 

bus 

bus 

R.(p.u.) 

X(p-u.) 

B(p.u.) 

73 

42 

90 

00810 

03880 

14030 

74 

42 

28 

00993 

04746 

38643 

75 

50 

30 

00468 

04770 

62450 

76 

50 

33 

00785 

07990 

1 04738 

77 

33 

64 

00806 

08169 

1 06808 

78 

54 

57 

01582 

03777 

.01735 

"79 

54 

75 

01303 

03102 

01429 

80 

57 

75 

02606 

06205 

02857 

81 

75 

88 

03583 

08532 

03929 

82 

29 

103 

01830 

09270 

07390 

83 

55 

59 

04188 

10518 

01087 

84 

59 

62 

05584 

14024 

01449 

85 

30 

36 

00294 

02997 

39206 

86 

94 

28 

01093 

05221 

18892 

87 

94 

118 

00662 

03164 

11451 

88 

32 

35 

00505 

02416 

08730 

89 

32 

52 

01291 

06171 

22327 

90 1 

34 

37 

08000 

20100 

02080 

91 

35 

38 

01600 

08100 ! 

06440 

92 

37 

39 

08560 

21500 

02220 

93 

38 

40 

01550 

07940 

06300 

94 

39 

41 

03820 

09090 

04180 

95 

41 

45 

04142 

10490 

04260 

96 

45 

49 

05680 

14260 

01470 

97 

45 

77 

03257 

07756 

03350 

98 

45 

156 

08000 

20100 

02080 

99 

156 

62 

04002 

10051 

01039 

100 

46 

78 

01660 

08430 

06720 

101 

46 

61 

01270 

06410 

05220 

102 

48 

49 

02280 

05710 

02370 

103 

48 

77 

02690 

06420 

02960 

104 

51 

52 

01770 

08510 

30260 

105 

52 

56 

01060 

05060" 

18320 

106 

56 

61 

00580 

02900 

02370 


Contd 
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Table J 4 (Contd ) . 


Line 1 

from 

To : 

Senes mipedence 

Shunt susceptance 

no 

bus 

>us 

El(p u ) : 

C(p.u.) 

B(p.u) 

107 

56 

60 

00370 

01780 

06440 

108 

60 

63 

00490 

02370 

08590 

109 

62 

65 

01120 T 

02800 

01160 

110 

63 

66 

00750 1 

03840 1 

03110 

111 

43 

63 

00679 T 

03412 

02782 

112 

43 

80 

00666 j 

03390] 

02672 

113 

63 

80 

01440 j 

07310 j 

05850 

114 

66 

67 

00670 

03390 

02670 

115 

66 

69 

00583 

02956 

02346 

116 

68 

79 

04560 

11450 

01180 

117 

69 

73 

01410 

07180 

05700 

118 

70 

77 

04080 

09640 

04440 

119 

70 

81 

01440 

03430 

01580 

120 

70 

24 

01260 

03160 

01300 

121 

24 

76 

00791 

01991 

00820 

122 

24 

72 

02230 

05620 

00580 

123 

72 

83 

01860 

04680 

01930 

124 

73 

91 

.01440 

.07250 

.05900 

125 

73 

i 87 

00705 

03590 

H 11400 

126 

52 

84 

00990 

05090 

04010 

127 

58 

84 

00780 1 

03980 

03140 

128 

76 

77 

08100 

20330 

02100 

129 

80 

82 

01160 

05830 

04750 

130 

82 

91 

00690 

03500 

11260 

131 

89 

91 

00050 

00253 

00805 

132 

93 

95 

00460 

01170 

00120 

133 

83 

95 

03910 

09840 

01010 

134 

83 

85 

01950 

04920 

02020 

135 

85 

95 

01870 

04700 

1 00490 

136 

9f 

93 

01031 

11874 

^ 00193 

137 

92 

112 

! 00479 

04880 

i 63614 

138 

90 

10E 

! 01732 

08784 

06973 

139 

30 

115 

> "00254 

02584 

r 33798 

140 

50 

92 

00051 

00511 

' 06760 

141 

29 

12' 

1 00580 

.0294C 

) 09420 

142 

71 

13 

l .06978 

1 1752! 

) 01817 


Contd. 
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Table J.4 (Contd.) . 


Line 

From 

To 

Senes impedence 

Shunt susceptance 

no 

bus 

bus 

R(pu) 

X(p.u.) 

B(p u.) 

143 

124 

126 

00830 

04240 

03340 

144 

92 

109 

00031 

00310 

04056 

145 

72 

76 

01390 

03510 

00360 

146 

117 

109 

00927 

09429 

1 23293 

147 

117 

109 

00833 

08478 

1 10855 

148 

117 

110 

00559 

05686 

.74354 

149 

36 

97 

00650 

06617 

86521 

150 

64 

97 

00260 

02646 

34610 

151 

36 

109 

00823 

08375 

1 09503 

152 

149 

55 

04560 

11453 

01183 

153 

149 

44 

04560 

11450 

01180 

154 

96 

130 

03723 

.08863 

.04082 

155 

130 

31 

00698 

.01662 

.00765 

156 

126 

129 

00930 

04750 

03740 

157 

129 

124 

01330 

.06680 

.05420 

158 

74 

64 

00370 

03762 

48870 

159 

74 

97 

00325 

.03307 

.43264 

150 

101 

102 

05030 

12620 

.01300 

161 

88 

98 

04372 

.08520 

03352 

162 

98 

106 

24980 

.04714 

01915 

163 

88 

98 

08744 

.17040 

.01676 

164 

98 

106 

04996 

09428 

00957 

165 

105 

106 

04930 

11750 

05410 

166 

104 

105 

05960 

14960 

01550 

167 

101 

104 

06140 

15460 

01590 

168 

71 

101 

05490 

13790 

01430 

169 

41 

100 

08470 

21320 

02190 

170 

49 

100 

05360 

13560 

01400 

171 

141 

99 

01163 

02770 

01276 

172 

86 

94 

00437 

02552 

09399 

173 

106 

107 

02050 

04880 

02250 

174 

63 

111 

00248 

01186 

04294 

175 

109 

113 

00031 

00310 

04056 

176 

113 

117 

00918 

09306 

1 21680 

177 

88 

115 

02311 

04361 

01772 

178 

29 

116 

. 01325 

06667 

05416 


Contd 
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Table J 4 (Contd.) 


Line 

From 

To 

Senes impedence 

Shunt susceptance 

no. 

bus 

bus 

R(p u ) 

X(p u ) 

B(p u.) 

179 

115 

153 

05956 

14959 

.01546 

180 

33 

117 

00015 

.00155 

02704 

181 

32 

118 

00124 

00593 

.02147 

182 

125 

71 

05119 

12856 

05332 

183 

86 

119 

00336 

01963 

.01808 

184 

86 

120 

01344 

07852 

.07230 

185 

35 

119 

01344 

07852 

07230 

186 

35 

120 

00336 

01963 

01808 

187 

46 

123 

01598 

08108 

06436 

188 

61 

67 

.01315 

06696 

05278 

189 

61 

60 

00499 

02397 

08523 

190 

58 

127 

.00998 

04794 

17046 

191 

110 

128 

00121 

01109 

72815 

192 

64 

26 

00051 

00517 

06760 

193 

108 

133 

01212 

06100 

04956 

194 

136 

134 

10889 

27347 

02826 

195 

136 

31 

09307 

23370 

02415 

196 

134 

135 

06140 

15460 

01590 

197 

34 

47 

02327 

05840 

00600 

198 

47 

138 

02327 

05840 

00600 

199 

34 

137 

09307 

22670 

02340 

200 

139 

34 

04188 

10518 

01087 

201 

140 

39 

09307 

22670 

02340 

202 

140 

150 

02280 

05725 

02360 

203 

141 

41 

02792 

07012 

02898 

204 

142 

37 

04188 

10518 

01087 

205 

143 

100 

04467 

11219 

01159 

206 

100 

49 

05380 

12840 

01480 

207 

100 

41 

08560 

21500 

02220 

208 

144 

151 

04560 

11450 

01100 

209 

145 

151 

04188 

10518 

01087 

210 

146 

151 

05680 

14260 

01470 

211 

147 

149 

03510 

08440 

00910" 

212 

147 

44 

04188 

10518 

01087 

213 

148 

45 

02880 

06860 

00790 

214 

150 

45 

02094 

05254 

02174 

215 

152 

134 

01115 

02810 

.01429 


Contd 
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Table J 4 (Contd ) 


Line 

From 

To 

Senes impedence j 

Shunt susceptance 

no 

bus 

bus 

R(pu ) 

X(p u ) 

B(p u ) 

216 

155 

123 

00878 

04430 

03580 

217 

132 

59 

04002 

.10051 

.01039 

218 

132 

156 

04188 

10518 

01087 

219 

122 

79 

03723 

09349 

00966 

220 

151 

121 

05960 

1496 

0155 

221 

71 

114 

08000 

20100 

0208 

222 1 

114 

27 

01440 

03430 

.01580 

223 

62 

53 

01440 

0343 

0158 

224 

25 

111 

00083 

00396 

01431 

225 

25 

63 

00695 

03500 

02843 
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Figure J 1 156-bus Indian system 
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